SH 
11 
A2 
ros 
FISH 


U.S. Department 
of Commerce 


Volume 122 
Number 1-2 
January-April 2024 


Fishery 
ulleti 


| = 
Ota £7 oo 


y= a “aay 

Mt mn Val 

CE 

Naf \ 
R= yu : 


d 
Be s CEE SWS ASS | 
ees Sits ears S8se 
Coc Soc 
e os eS Seve SS 
a —— \ tg DN Ss 
ONT | ne nal A 


¥ cme 
ES -=< (SF 
ON | A 


U.S. Department 
of Commerce 


Gina M. Raimondo 
Secretary 


National Oceanic 
and Atmospheric 
Administration 


Richard W. Spinrad 


Administrator 


National Marine 
Fisheries Service 


Janet Coit 
Assistant Administrator 
for Fisheries 


Fishery Bulletin (ISSN 0090-0656) 
is published quarterly by the Scientific 
Publications Office, National Marine 
Fisheries Service, NOAA, 7600 Sand 
Point Way NE, Seattle, WA 98115-0070. 

Although the contents of this publica- 
tion have not been copyrighted and may 
be reprinted entirely, reference to source 
is appreciated. 

The Secretary of Commerce has deter- 
mined that the publication of this peri- 
odical is necessary according to law for 
the transaction of public business of this 


Scientific Editor 

José |. Castro 

National Marine Fisheries Service, NOAA 
263 13th Avenue South 

Saint Petersburg, Florida 33701 


Managing Editor 

Kathryn Dennis 

Scientific Publications Office 

National Marine Fisheries Service, NOAA 
1845 Wasp Boulevard, Building 176 
Honolulu, Hawaii 96818 


Associate Editor 

Cara Mayo 

Scientific Publications Office 

National Marine Fisheries Service, NOAA 
7600 Sand Point Way NE 

Seattle, Washington 98115-0070 


Editorial Board 


Henry L. Bart Jr. 
Katherine E. Bemis 
Matthew D. Campbell 
William B. Driggers III 
Gretchen L. Grammer 
Eric R. Hoffmayer 
Richard Langton 

J. Fernando Marquez-Farias 
Richard S. McBride 
Richard L. Merrick 
Richard D. Methot 
Lisa J. Natanson 

Mark S. Peterson 
André E. Punt 

Joseph M. Quattro 
John F. Walter III 


Tulane University, New Orleans, Louisiana 

National Systematics Laboratory, Washington, DC 

National Marine Fisheries Service, Pascagoula, Mississippi 

National Marine Fisheries Service, Pascagoula, Mississippi 

South Australian Research and Development Institute, Adelaide, Australia 
National Marine Fisheries Service, Pascagoula Mississippi 

National Marine Fisheries Service, Walpole, Maine (retired) 
Universidad Autonoma de Sinaloa, Mazatlan, México 

National Marine Fisheries Service, Woods Hole, Massachusetts 
National Marine Fisheries Service, Falmouth, Massachusetts (retired) 
National Marine Fisheries Service, Seattle, Washington 

National Marine Fisheries Service, Narragansett, Rhode Island (retired) 
University of Southern Mississippi, Ocean Springs, Mississippi 
University of Washington, Seattle, Washington 

University of South Carolina, Columbia, South Carolina 

National Marine Fisheries Service, Miami, Florida 


Fishery Bulletin website: https://spo.nmfs.noaa.gov/fb.htm 


Department. Use of funds for printing of 
this periodical has been approved by the 
Director of the Office of Management and 
Budget. 

For Sale by the Superintendent of 
Documents, U.S. Government Publishing 
Office, Washington, DC 20401. Subscrip- 
tion price per year: $32.00 domestic and 
$44.80 foreign. Cost per single issue: 
$19.00 domestic and $26.60 foreign. See 
back for order form. 


Fishery Bulletin carries original research reports on investigations in fishery science, engineering, and 
economics. It began as the Bulletin of the United States Fish Commission in 1881; it became the Bulletin 
of the Bureau of Fisheries in 1904 and the Fishery Bulletin of the U.S. Fish and Wildlife Service in 1941. 
Separates were issued as documents through volume 46; the last document was no. 1103. Beginning with 
volume 47 in 1931 and continuing through volume 62 in 1963, each separate appeared as a numbered 
bulletin. A new system began in 1963, with volume 63, in which papers are bound together in a single issue. 
Beginning with volume 70, number 1, January 1972, Fishery Bulletin became a periodical, issued quarterly. 
In this form, it is available by subscription from the Superintendent of Documents, U.S. Government 
Publishing Office, Washington, DC 20401. It is also available free in limited numbers to libraries, 
research institutions, and state and federal agencies and in exchange for other scientific publications. 


U.S. Department 


of Commerce 
Seattle, Washington 


Volume 122 
Number 1-2 
January-April 2024 


The National Marine Fisheries 
Service (NMFS) does not approve, 
recommend, or endorse any proprie- 
tary product or proprietary material 
mentioned in this publication. No 
reference shall be made to NMFS, 
or to this publication furnished by 
NMFS, in any advertising or sales 
promotion that would indicate or 
imply that NMFS approves, rec- 
ommends, or endorses any propri- 
etary product or proprietary mate- 
rial mentioned herein or that has 
as its purpose an intent to cause 
directly or indirectly the advertised 
product to be used or purchased be- 
cause of this NMFS publication. 


The NMFS Scientific Publications 
Office is not responsible for the con- 
tents of the articles. 


Fishery 
Bulletin 


Contents 
Articles 
1-12 Jeff Kneebone, Connor F. White, Caroline Collatos, and 


Nicholas M. Whitney 


High survivorship of sandbar sharks (Carcharhinus plumbeus) following catch 
and release in a growing land-based fishery in Massachusetts 


13-25 Valéria M. Lemos, Ana C. G. Mai, Ana C. Martins, José V. Martins, 
Sandra Andrade, Henrique Cabral, and Joao P. Vieira 


Variation in otolith microchemistry for Lebranche mullet (Mugi/ liza) in 
southern Brazil and its potential use in identifying their nursery grounds 


26-43 Matthew D. Damiano, Mandy Karnauskas, Wessley Merten, and Jie Cao 


Spatiotemporal dynamics of dolphinfish (Coryphaena hippurus) in the western 
Atlantic Ocean 


44-57 Benjamin R. LaFreniere, Briony Donahue, Jillian E. Price, Alicia Cruz-Uribe, 
Nate Miller, Benjamin T. Manard, Richard McBride, and John A. Mohan 


Chemical clocks: using otolith geochemistry to enhance estimation of age and 
growth of white hake (Urophycis tenuis) 


59-62 Guidelines for authors 


JUL 24 2024 
LIBRARIES 


— 


National Marine 
Fisheries Service 


NOAA 


Abstract—Directed land-based rec- 
reational catch-and-release fishing 
for sandbar sharks (Carcharhinus 
plumbeus) is growing in popularity in 
Massachusetts. Working with 21 vol- 
unteer fishermen of varying experience 
levels, we observed and documented 
the fishing gear, tackle, and techniques 
used to catch and release 67 sand- 
bar sharks. The postrelease fate of 
each shark was monitored by using 
an acceleration data logger (ADL) tag 
embedded in a custom float package 
that was secured to the first dorsal fin 
with a galvanic timed release. All 67 
packages were recovered after detach- 
ment following monitoring periods of 
0.15—-9.98 d. Examination of the depth, 
tailbeat period (TBP), pitch, and roll 
time series from 65 ADLs that recorded 
data revealed high survivorship for 
tagged sharks, with all of them alive 
at the time of tag detachment. Behav- 
ioral recovery was estimated to have 
occurred at an average of 6.36 h after 
release on the basis of trends in TBP 
for 54 sandbar sharks with at least 10 h 
of postrelease acceleration data. These 
results indicate that sandbar sharks 
are remarkably resilient to catch and 
release in the land-based shark fishery 
in Massachusetts. 
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Recreational shark fishing is a popular 
activity with significant socioeconomic 
value to coastal communities along the 
East Coast of the United States and 
in the Gulf of Mexico region (Hutt and 
Silva, 2019; Guay et al., 2021). Tradi- 
tionally, a large portion of the recre- 
ational shark fishing effort has occurred 
at sea aboard federally permitted char- 
ter boats and private fishing vessels. 
However, a recent rise in the popularity 
of shore- or land-based shark angling 
has resulted in increased effort within 
state waters (senior author, personal 
observ.; Guay et al., 2021). Because of 
how recreational catch and effort data 
are collected by the National Marine 
Fisheries Service Marine Recreational 
Intercept Program and Large Pelagics 
Survey, as well as the lack of a federal 
permit requirement for land-based 
shark fishing, data on the magnitude 
of catch or effort across this fishery are 
very limited (Weber et al., 2020; Guay 
et al., 2021). Nonetheless, it is well- 
documented that several shark species 
targeted or incidentally captured by 


land-based anglers are experiencing 
population declines and are prohibited 
from retention under state and federal 
laws (Kilfoil et al., 2017; Gibson et al., 
2019). One such species is the over- 
fished sandbar shark (Carcharhinus 
plumbeus) (SEDAR, 2017). 

Directed land-based recreational fish- 
ing for sandbar sharks has occurred in 
Massachusetts for decades (Skomal, 
2007). However, evidence indicates 
that land-based fishing effort for sand- 
bar sharks and shark catch rates have 
increased dramatically in Massachu- 
setts over the last decade (Skomal’). 
For example, conversations with recre- 
ational land-based anglers have revealed 
that catches of 5-10 sandbar sharks in 
a single trip are now common (senior 
author, personal observ.), and individual 
land-based anglers have convention- 
ally tagged over 100 sandbar sharks in 


1 Skomal, G. 2023. Personal commun. Mass. 
Div. Mar. Fish., 8836 S Rodney French 
Blvd., New Bedford, MA 02744. 


a single season (Zewinski”). Under current Massachusetts 
state law, circle hooks must be used when targeting sharks 
with natural baits, and all sandbar sharks must be released 
(Regulation .. . 2023). However, anglers in Massachusetts 
are not required to possess an additional permit to engage 
in land-based shark angling nor undergo training on catch- 
and-release best practices prior to engaging in the fishery; in 
contrast, both a permit and training are required in Florida 
(Sharks ... 2019). 

In few studies have the effects of land-based capture and 
handling on shark mortality rates been examined. Weber 
et al. (2020) reported a postrelease mortality (PRM) rate 
of 17% for blacktip sharks (C. limbatus) caught in South 
Carolina, with mortality occurring because of physical 
trauma incurred during capture and foul hooking. Simi- 
larly, Kilfoil et al. (2017) reported low (6%) PRM of sand 
tigers (Carcharias taurus) caught in Delaware, with all 
mortality attributed to internal hooking. Most recently, 
Binstock et al. (2023) reported variable, species-specific 
PRM rates for bull (Carcharhinus Leucas) (7%), blacktip 
(50%), and tiger (Galeocerdo cuvier) (0%) sharks caught in 
Texas, with evidence that increasing water temperatures 
resulted in higher mortality rates for blacktip sharks. 
Mortality has also been suspected to have occurred in 
land-based fishing operations because of the removal of 
sharks from the water (i.e., air exposure) for extended 
periods during handling (Shiffman et al., 2017). However, 
to date, no direct relationship between time out of water 
and mortality has been reported. Because PRM rates for 
sharks are highly specific to species (e.g., Whitney et al., 
2021) and geographic location (e.g., Mohan et al., 2020), 
there is a continued need to estimate these rates across all 
species, fisheries, fishing areas, and environmental condi- 
tions. Accordingly, the goals of this study were to assess 
the PRM rate of sandbar sharks caught in the land-based 
shark fishery in Massachusetts and to identify factors 
that adversely affect animal condition, postrelease behav- 
ior, and survival. We also briefly assessed the effect of ani- 
mal behavior on tag retention. 


Materials and methods 


Activities of this project were conducted under the New 
England Aquarium Animal Care and Use Committee pro- 
tocol 2019-06. 


Study area and design 


All project activities were conducted in cooperation 
with volunteer recreational land-based shark fishermen 
on Cape Cod and Nantucket in Massachusetts. Scien- 
tific researchers accompanied and observed fishermen 
during all trips, but anglers were free to choose the loca- 
tion of angling and received no input from scientific staff 


2 Zewinski, K. 2021. Personal commun. Popul. Ecosyst. Monit. 
Anal. Div., Northeast Fish. Sci Cent., Natl. Mar. Fish. Serv., 
NOAA, 28 Tarzwell Dr., Narragansett, RI 02882. 
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regarding the gear (e.g., rod and reel type or size), tackle 
(e.g., hook type or size), or techniques used to capture and 
handle sandbdar sharks. Each participant’s experience 
with land-based shark fishing was determined at the time 
of each trip on the basis of the following criteria: 


Novice: less than 1 year of experience with land-based 
shark fishing or a total of <5 sharks captured and 
handled from land in their angling career; 

Intermediate: more than 1 year of experience with 
land-based shark fishing or a total of 5-10 sharks 
captured and handled from land in their angling 
career; and 

Expert: more than 2 years of experience with land- 
based shark fishing or a total of >10 sharks captured 
and handled from land in their angling career. 


When a volunteer angler captured a sandbar shark, 
the following data were collected: rod and reel type (e.g., 
conventional versus spinning) and size (e.g., model or 
pound class); line type and strength; leader type, length, 
and strength; hook type (circle or J) and size; and bait 
type. Fight time (i.e., the time from when the shark was 
hooked to the time it was landed) was recorded to the 
nearest second. Handling location, or position of the shark 
during landing, unhooking, and tagging, was recorded as 1 
of 3 locations: shark left in water (i.e., the majority of the 
body, including the head and gills, remained in the water), 
shark brought to the water line (i.e., only a portion of 
the animal remained in the water; this location includes 
instances when wave action briefly exposed the head and 
gills to air), or shark removed from water (i.e., the entire 
shark was out of the water for the entirety of the handling 
event). Air and water temperature at the time of capture 
were also noted. 

Biological assessment and tagging activities were con- 
ducted where the shark was landed, and efforts were made 
to note the sex, measure the fork length (FL) to the near- 
est centimeter, and apply a tag package that included an 
acceleration data logger (ADL) (methods described in the 
next section) while fishermen were removing or attempt- 
ing to remove the hook, so as not to artificially extend 
handling time. Hooking location was noted (jaw, inside 
mouth [esophagus or pharynx], gills, gut [hook not visible 
or accessible], or body [foul hooked]) along with the dis- 
position of the hook (removed or retained). The presence, 
location, and severity of bleeding (none, light, moderate, or 
heavy), overt signs of physical trauma (e.g., abrasions and 
lacerations), and the extent to which fishermen physically 
restrained the shark (e.g., holding the body of the shark 
in place during unhooking) were also noted. Handling 
time (i.e., the time from landing to release, including time 
needed for tagging) was recorded to the nearest second, 
and total capture and handling duration was noted as the 
sum of fight and handling times. Following tagging, fish- 
ermen were responsible for releasing the shark and left 
to determine if revival was necessary. Revival time was 
noted to the nearest second, if applicable. Finally, each 
shark was assigned a release condition based in part on 
the ordinal scale reported in Weber et al. (2020) (Table 1). 
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Table 1 


Release condition categories assigned to sandbar sharks 
(Carcharhinus plumbeus) captured off Cape Cod and Nan- 
tucket in Massachusetts from June through September in 
2019-2021. Conditions are based in part on those used by 
Weber et al. (2020). 


Physical condition and release behavior 


Category 


No physical trauma other than hook wound, 
no bleeding, swam away rapidly 

Good No physical trauma other than hook 
wound, no bleeding, swam away slowly or 
appeared disoriented 

Fair Overt signs of physical trauma (abrasions 
and minor cuts), minor bleeding, swam 
away slowly or appeared disoriented 

Poor Overt signs of physical trauma (abrasions, 

minor cuts, and lacerations), moderate 

or severe bleeding, swam away slowly or 

appeared disoriented 


Excellent 


Postrelease monitoring 


To monitor the postrelease behavior and fate of sandbar 
sharks, custom tag packages that included an ADL and a 
float were created following the design outlined in Whit- 
more et al. (2016) and Whitney et al. (2021). Each pack- 
age consisted of an ADL (Axy 4°. TechnoSmArt, Rome, 
Italy) that recorded triaxial acceleration at 25 Hz and 
depth and temperature every second, a SPOT-6 satellite 
transmitter (Wildlife Computers Inc., Redmond, WA), and 
a VHF transmitter (Advanced Telemetry Systems Inc., 
Isanti, MN). Each of these components was embedded in 
a custom tag float (107 x 68 x 50 mm in size and 151 g in 
weight in air). The package was designed to be <2% of the 
body weight of a tagged sandbar shark, meeting the rec- 
ommended threshold intended to minimize the likelihood 
of tag-induced impairment to movement (Use of Fishes in 
Research Committee, 2014). 

During tagging of each shark, 2 holes were drilled in 
the first dorsal fin, and a float package was secured on 
the left side of the fin (relative to the cranial—caudal axis) 
with cable ties, vinyl Peterson discs (Floy Tag & Mfg. 
Inc., Seattle, WA), and a galvanic timed release (GTR) 
(International Fishing Devices Inc., Northland, New 
Zealand). The GTRs dissolved in seawater and were sized 
in a manner to keep the ADL packages secured to sharks 
for a predetermined duration of approximately 1-4 d, 
and once that time was reached, the tags would release 
from the shark and float to the surface (Whitmore et al., 
2016). The satellite transmitter would report its location 
through the Argos satellite network once it had reached 
the surface, and land- or sea-based excursions were made 


3 Mention of trade names or commercial companies is for identi- 
fication purposes only and does not imply endorsement by the 
National Marine Fisheries Service, NOAA. 


to recover the tag by using a multichannel VHF receiver 
(Advanced Telemetry Systems) following the methods 
described by Lear and Whitney (2016). 


Data analysis 


All analyses were performed by using the statistical soft- 
ware R (vers. 4.1.2; R Core Team, 2021), and statistical 
significance was accepted at P<0.05. 


Capture and handling characteristics A series of qualita- 
tive and quantitative assessments were used to exam- 
ine basic relationships between physical and biological 
aspects of the capture and handling events that may 
relate to release condition and PRM of sandbar sharks. 
Because of the opportunistic nature of sampling in this 
study, sample sizes were heavily unbalanced among many 
capture- and handling-related variables (e.g., line class, 
hooking location, and handling location) and autocorrela- 
tion structures violated model assumptions about sample 
independence (e.g., some anglers used only a single type or 
size of hook or line strength). Accordingly, box plots were 
created to qualitatively demonstrate general trends in the 
relationship between categorical variables related to cap- 
ture and handling, with individual categories containing 
<5 observations. 

Individual generalized linear models (GLMs) were 
constructed to quantitatively examine the effect of con- 
tinuous and categorical variables with a sufficient num- 
ber of observations on fight time and handling time. In 
fight time GLMs, data from a single anomalous cap- 
ture event with an artificially inflated fight time were 
excluded from analysis (for a description of this event, 
see the “Results” section). Handling time GLMs were 
constructed with FL and restraint (binomial) as fixed 
effects. Data from 2 handling events were excluded 
from analysis because difficulties in the application of 
the ADL package artificially extended handling time 
beyond what was necessary for fishermen to prepare the 
shark for release. All GLMs were fit with data that had 
a Gaussian distribution in the R package Ime4 (vers. 1.1- 
30; Bates et al., 2015). Model fit was assessed by visually 
examining residual plots. 


Postrelease fate and mortality and survival rates After recov- 
ery of the ADL package, acceleration, depth, and tem- 
perature data were downloaded from the ADL, and the 
time of detachment of the package was noted as the time 
when the tag achieved an upright position and floated 
to the surface. To correct for variation in tag placement 
on the fin and attachment angle, acceleration data were 
rotated to be in the same orientation across all deploy- 
ments. Rotation angles were determined by ensuring that 
the average acceleration of the x- and y-axis was close 
to 0; the x-axis correlated with the vertical velocity (i.e., 
ascending or descending in the water column), and the 
y-axis contained the largest dynamic acceleration. This 
approach ensured that the x-axis corresponded to the 
cranial—caudal axis of the animal (this axis pitches up 
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and down as the animal ascends and descends), the y-axis 
corresponded to the lateral axis of the animal (which 
experiences dynamic acceleration during tailbeats), and 
the z-axis corresponded to the dorsal—ventral axis of the 
animal (Jorgensen et al., 2015). 

Static acceleration (accel) was calculated by using a 3-s 
(75-point) running mean of the raw data, and dynamic 
acceleration was calculated as the remainder when sub- 
tracting the dynamic acceleration from the raw accel- 
eration (Shepard et al., 2008). Body pitch and roll were 
calculated from the static acceleration, and the overall 
dynamic body acceleration (ODBA) (Wilson et al., 2006) 
was calculated from the dynamic acceleration: 


Roll = tan” Ca Z seep (1) 
Y. 


Pitch = tan* (-X, aces] Sin Roll (2) 


ccel> 


+ Laece) COS Roll); and 


ODBA = abs(Xnnca) + ADS(Y accep) + ADS(Zaccey)- (3) 


A continuous wavelet transformation of the y-axis accel- 
eration data was used to calculate the power of period- 
icities over a period of 0.2—-4.0 s by using the R package 
biwavelet, vers. 0.20.21 (Gouhier et al., 2021). The period- 
icity with the greatest power each second was ascribed as 
the individual’s tailbeat period (TBP). 

Postrelease fate was determined through visual exam- 
ination of fine-scale vertical movement and acceleration 
profiles recorded by the ADL (Whitney et al., 2016, 2017, 
2021). The sandbar shark is an obligate ram-ventilating 
species, and the fate of each shark at the time of tag 
detachment was assessed by using the following criteria: 


Alive: consistent tailbeat patterns and changes in 
depth (i.e., vertical movements) and body position 
(i.e., pitch and roll) indicative of active swimming; 

Dead: prolonged cessation of tailbeats in concert with a 
constant depth (indicative of rest on the seafloor); and 

Predation: shift in the pitch and roll, an alteration of 
TBP and dive patterns, and potentially a lack of ver- 
tical thermal variation, which can indicate presence 
of the tag in a predator’s stomach. 


Total PRM was calculated as the proportion of all mor- 
tality events observed across monitored individuals. We 
calculated 95% confidence intervals (CIs) by using the Wil- 
son method (Brown et al., 2001) in the R package binom, 
vers. 1.1-1 (Dorai-Raj, 2014). 


Postrelease recovery period Behavioral recovery periods 
were evaluated for surviving sandbar sharks to address 
the sublethal effect of capture and handling on swimming 
or predator avoidance capabilities. Results from previous 
studies indicate that TBP is a suitable metric for assess- 
ing postrelease behavioral changes (e.g., Whitney et al., 
2017, 2021; Grainger et al., 2022). Therefore, average 
TBP was calculated every 10 min from the time of release 
to 48 h after release for all sharks. With this approach, it 
has been found that behavioral recovery often occurs less 


than 10 h after release in multiple species (e.g., Whitney 
et al., 2017, 2021). Therefore, recovery time was not cal- 
culated for any sharks for which less than 10 h of postre- ° 
lease acceleration data were available. Patterns in TBP 
were modeled by using nonlinear mixed effects models as 
an asymptotic relationship executed in the lme4 package 
(Bates et al., 2015): 


TBP = Ag, +(Ro Avs) ( - el (4) 


where A,,,, = the fully recovered TBP; 
Ry = the TBP immediately after release; 
Irc = the log rate constant that describes the rate 
at which recovery occurs; and 


T = the time after release. 


To account for the repeated measures of TBP for each 
individual over 48 h, we included random effects so that all 
3 model parameters (A,,,,,, Ro, and /rc) varied by individ- 
ual. Recovery period was defined as the time after release 
when TBP reached 80% of the difference between Ry and 
Agym (Whitney et al., 2016, 2017, 2021). By incorporating 
individual as a random effect, time to recovery could be 
calculated for each shark. 

Fixed-effect GLMs were fit to time-to-recovery data to 
assess the influence of FL, fight time, handling time, and 
total time by using the Ime4 package (Bates et al., 2015). 
Because of low and unequal samples sizes evident in some 
conditions nested in categorical variables (e.g., recovery 
time was only estimated for 2 sharks in poor condition and 
for 3 sharks that were left in the water), the following vari- 
ables were excluded from models: hook disposition, han- 
dling location, and fish condition. Instead, box plots were 
created to qualitatively describe relationships between 
these variables and recovery period. 


Tag detachment The results of the preliminary analy- 
sis of pitch, roll, and ODBA data revealed that numer- 
ous tagged sandbar sharks exhibited rolling behavior 
in which the body of the shark temporarily rolled >75° 
to the left (along the axis of forward motion, to the side 
with the ADL package) or to the right while the animal 
was actively swimming. In addition to these observations, 
multiple recovered float packages had scratches and 
abrasions, indicating that the rolling behavior might have 
resulted in contact of the tag with the seafloor. To investi- 
gate the extent to which rolling was linked to premature 
tag detachment, the depth, pitch, and roll of each individ- 
ual were visually inspected for the 50 s surrounding the 
time of tag detachment. Tag detachment was identified as 
a positive pitch, indicating that the tag nose was pointed 
to the seafloor, while the depth was decreasing (i.e., the 
tag was floating to the surface). Rolling behavior was 
identified if individuals rolled to either the left or right 
side (>75°) while pitch values indicate that they were 
horizontal (<25° and >—25°) (Suppl. Figure) (online only). 
Additionally, this rolling behavior had to occur when the 
shark was presumed to be at the seafloor (i.e., within 1 m 
of the maximum depth the shark experienced in the 5 min 
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before and after the roll) and <2 s prior to the estimated 
time of tag detachment. Each individual’s tag detachment 
was categorized as associated with a left roll (<—75°) or a 
right roll (>75°) or was not associated with a roll. 

To determine if rolls influenced tag retention, we ana- 
lyzed the estimated error in tag retention across the 3 tag 
detachment categories. Release error was calculated as 
the duration the tag remained on the animal divided by 
the expected duration based on GTR size. Releases were 
classified as premature if the tag detached more than 25% 
(i.e., <-25%) earlier than expected. Release error was then 
compared across the 3 detachment categories to deter- 
mine if rolling behavior immediately before a tag detach- 
ment influenced tag retention. Pairwise comparisons were 
performed by using Tukey’s honestly significant difference 
(HSD) test. Data from tags that were retained on sharks 
for greater than 100% longer than the estimated GTR 
release time were excluded from analysis and were con- 
sidered outliers. 


Results 


A total of 67 sandbar sharks (47 females and 20 males), 
measuring 113-187 cm FL (mean: 140 cm FL [standard 
deviation (SD) 17]), were captured by 21 volunteer recre- 
ational fishermen (9 novice, 7 intermediate, and 5 expert 
anglers) and were tagged with ADL packages from June 
through September in 2019 (sample size [n]=14), 2020 
(n=8), and 2021 (n=45) (Suppl. Table) (online only). Females 
and males ranged from 113 to 187 cm FL (mean: 144 cm 
FL [SD 18]) and from 113 to 152 cm FL (mean: 132 cm 
FL [SD 12]), respectively. Individual anglers were observed 
capturing 1—20 sharks (mean: 3 sharks [SD 4]), with expert 
anglers accounting for the largest number of sharks cap- 
tured (n=39). Novice and intermediate anglers captured 17 
and 11 sharks, respectively. 

To capture sandbar sharks, fishermen used a range of 
bait and tackle, including 13-36 kg (30-80 lb) class con- 
ventional (n=1) and spinning (n=66) rod and reel setups 
spooled with 13-36 kg braided main fishing line and steel 
single-strand wire (n=33) or multistrand cable (n=34) 
leaders with lengths of 0.6—1.2 m (~2—4 ft) and rated at 
weight test of 36-181 kg (80-400 lb). Sharks were pri- 
marily captured with 8/0—18/0 circle hooks (n=58, 86.6% 
of capture events), but 9 sharks (13.4% of events) were 
captured on 8/0—10/0 J-hooks (Table 2). To capture sharks, 
volunteers used a variety of natural baits, including cut 
sections of Atlantic menhaden (Brevoortia tyrannus), 
bluefish (Pomatomus saltatrix), and little tunny (Euthyn- 
nus alletteratus) and whole live or dead American eels 
(Anguilla rostrata). 

Sandbar sharks were primarily hooked in the jaw (n=59, 
88.1%) but were also hooked inside the mouth (n=2), in 
the gills (n=4), and on the body (pectoral fin, n=2). All 
hooking locations were observed with both circle and 
J-hooks; however, a larger percentage of sharks caught 
on J-hooks were hooked inside the mouth, in the gills, or 
on the body (Fig. 1). Anglers removed hooks from 92.5% 


Table 2 


Breakdown of data collected for the 67 sandbar sharks 
(Carcharhinus plumbeus) for which land-based cap- 
ture and handling events were observed during June— 
September in 2019, 2020, and 2021 in Massachusetts. 


Observations 


Variable Number Percentage 
Sex 
Females 70.1 
Males 29.9 
Hook type 
Circle 
J 
Hooking location 
Jaw 
Mouth 
Gills 
Body 
Hook disposition 
Removed 
Retained 
Handling location 
Brought out of water 
Brought to water line 
Left in water 
Shark restrained 
Yes 
No 
Fish condition at release 
Excellent 
Good 
Fair 
Poor 
Bleeding 
None 
Light 
Moderate 
Heavy 


of captured sharks (n=62); the leader was cut and the 
hook retained in 5 sharks caught on circle hooks, 3 sharks 
hooked in the jaw and 2 sharks hooked in the gills. The 
majority of sharks had no bleeding or had minor bleed- 
ing localized around the hook wound. Two sharks hooked 
in the gills experienced heavy bleeding. In one of these 
animals, the J-hook was removed by using a dehooking 
device, and in the other, the leader was cut, leaving the 
circle hook embedded in the gills. Fight times ranged from 
2.57 to 32.40 min (mean: 6.90 min [SD 4.75]), and 98.5% 
of fight times were less than 20 min. The longest fight 
time (32.40 min) was observed when a fisherman, fishing 
alone, hooked 2 sandbar sharks at the same time and did 
not begin retrieving the second animal until the first was 
released. Excluding this anomalous capture event, there 
was a significant positive relationship between shark 
FL and fight time (coefficient of multiple determination 
[R?]=0.306, P<0.001) (Fig. 2A). 
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Figure 1 


Relative percentage of hook locations observed by hook 
type for 67 sandbar sharks (Carcharhinus plumbeus) cap- 
tured off Cape Cod and Nantucket in Massachusetts from 
June through September in 2019-2021. Fifty-eight sharks 
were caught on circle hooks, and 9 sharks were caught on 
J-hooks. Use of circle hooks yielded a higher percentage 
of capture events in which sharks were hooked in the jaw, 
and use of J-hooks had a higher probability of hooking 
sharks on the body, in the gills, or inside the mouth. 


Once brought into the surf zone, sharks in the majority 
of cases were dragged to the water line (n=39, 58.2%) or 
completely out of the water (n=23, 34.3%) by fishermen 
holding them by the tail. Anglers left the shark in the 
water in 5 instances (7.5% of observed capture events). 
Forty-seven sharks (70.1%) were restrained during han- 
dling. Handling times ranged from 1.50 to 9.25 min (mean: 
3.82 min [SD 1.27]), yielding total capture and handling 
event times of 5.25—-34.33 min (mean: 10.72 min [SD 4.96]). 
There was a significant, positive relationship between FL 
and handling time (R?=0.131, P<0.001) (Fig. 2B), but han- 
dling time was not significantly affected by the restraint of 
sharks (P=0.152). Qualitatively, sharks that were handled 
in the water or that were hooked in the gills generally 
experienced longer handling times (Fig. 3). At the time of 
release, all sharks were alive and assigned the following 
release conditions: excellent (n=16), good (n=44), fair (n=5), 
and poor (n=2). During handling, sharks experienced air 
temperatures ranging from 15.9°C to 28.9°C (mean: 23.4°C 
[SD 2.9]) and were captured in water temperatures rang- 
ing from 21.2°C to 27.7°C (mean: 24.1°C [SD 1.7]). 


Tag recovery and performance 


All 67 ADL packages deployed on sandbar sharks detached 
and were recovered. Forty-one tags were recovered at 
sea, and 26 tags were recovered after washing ashore on 
Cape Cod or Nantucket. The ADLs remained attached to 


sandbar sharks for periods of 0.15—-9.98 d (mean: 2.67 d 
[SD 2.08]) or 3.6—239.5 h (mean: 64.0 h [SD 49.9]). Two 
ADLs malfunctioned and did not record any data during 
the deployment period, and 3 ADLs stopped recording 
prior to detaching. Altogether, data sufficient for determin- 
ing postrelease fate were available from 65 ADLs. A total 
of 27 tags (40.3%) released prematurely, earlier than their 
expected time based on the size of the GTR used. 


Postrelease survival and recovery 


Results from the examination of the depth, TBP, pitch, and 
roll time series from the 65 ADLs reveal high survivorship 
for tagged sharks. All sharks were alive at the time of tag 
detachment, yielding a PRM rate of 0% (95% CI: 0.0-5.6%) 
or survival rate of 100% (95% CI: 94.4—100%). 

Time to recovery was estimated for 54 of 56 sandbar 
sharks with at least 10 h of postrelease acceleration data. 
Recovery time was not estimated for 2 sharks whose TBP 
failed to reach their recovery value (i.e., 80% of the fully 
recovered TBP) and whose change in TBP was in the low- 
est 25th percentile during the time from release to 48 h 
after release. Of note, both of these sharks consistently had 
a TBP that was similar to the “recovered” values of the 54 
sharks for which recovery time was estimated. Behavioral 
recovery in TBP for the 54 sharks was evident from 1.33 
to 17.18 h after release (mean: 6.36 h [SD 3.52]) (Fig. 4). 
Results from the GLMs reveal no significant relationships 
between recovery time and fight time, handling time, total 
time, or FL (P>0.11). However, findings from qualitative 
analyses with box plots provide some evidence that sharks 
hooked in the gills and classified in fair or poor condition 
took longer to recover (Fig. 5). 


Tag detachment 


Tag detachment in association with rolling behavior 
was examined for 60 individuals whose ADLs recorded 
data at the time of release. In 32 instances (53.3%), tags 
appeared to detach without any clear indication of roll- 
ing behavior. These tags detached close (mean: —9.2% 
[SD 36.1]) to the expected GTR duration, and only 25.0% 
(n=8) of these tags released prematurely (Fig. 6). Rolling 
was observed in 28 instances (46.7%) of tag detachment, 
including 11 instances (39.3%) when sharks rolled to the 
right (non-tag side) and 17 cases (60.7%) when sharks 
rolled to the left side (tag side). Only 4 tag detachments 
associated with rolls to the right (36.4%) occurred prema- 
turely, a number that is not significantly different than 
the number of premature detachments not associated 
with rolling behavior (n=6) (Tukey’s HSD test: P=0.99). 
However, 88.2% (n=15) of tag detachments associated 
with rolls to the left side occurred prematurely, signifi- 
cantly earlier than detachments associated with rolls to 
the right (Tukey’s HSD test: P<0.001) and with detach- 
ments that were not associated with rolling behavior 
(Tukey’s HSD test: P<0.001). In total, 55.6% (n=15) of all 
premature releases (n=27) occurred when sharks rolled 
to their left side. 
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Figure 2 


Relationship of (A) fight time (sample size [n]=66) and (B) handling time (n=65) with fork length 
(FL) for sandbar sharks (Carcharhinus plumbeus) captured in Massachusetts during 2019-2021. 
The gray line represents the linear relationship, and the gray shaded area indicates the standard 
error around the predicted relationship. In panel A, data for a single 32-min event of the capture 
of a 132-cm-FL female shark was excluded (and not plotted) because the volunteer angler hooked 
2 sharks simultaneously and did not begin reeling this shark in until the first one was released. 
In panel B, data for 2 handling events were excluded because of difficulties in applying the accel- 
eration data logger package that artificially extended handling time. 


Discussion 


The lack of mortality (100% survival) observed in sand- 
bar sharks following land-based catch-and-release activ- 
ities adds to the growing body of evidence that this 
species is resilient and capable of surviving both the 
physical and physiological effects of capture and han- 
dling (e.g., Marshall et al., 2015; Whitney et al., 2021). 
Postrelease mortality of sharks caught on recreational 
hook-and-line fishing gear has been linked to physi- 
cal trauma from internal hooking in the gut or gills 
(Kneebone et al., 2013; French et al., 2015; Kilfoil et al., 
2017), foul hooking on the body (Sepulveda et al., 2015; 
Weber et al., 2020), extended fight times due to foul 
hooking (Heberer et al., 2010), and prolonged handling 
times (Mohan et al., 2020; Knotek et al., 2022). However, 
in our study, no mortality was identified in 2 sandbar 
sharks that were hooked in the gills and were bleeding 
severely at the time of release, with one foul-hooked 
individual having been fought for 32.40 min (i.e., nearly 
5 times longer than the mean fight time observed across 
all sharks). In addition, no mortality was observed in 
23 individuals that were completely removed from the 
water for 2.20-6.67 min during handling. These observa- 
tions and the finding that the majority (90%) of sharks 
monitored in our study were released in excellent or 


good condition indicate that sandbar sharks are likely to 
survive land-based catch-and-release fishing in Massa- 
chusetts, as long as best practices are used consistently 
throughout the fishery. 

Despite the resilience of this species, it should be under- 
stood that the true PRM rate of sandbar sharks caught 
and released by using land-based techniques in Massa- 
chusetts is most certainly greater than zero. Indeed, our 
research team necropsied 2 sandbar sharks that washed 
ashore 24-36 h after they were caught and released by 
land-based fishermen, as confirmed by the presence of an 
M-tag of the National Marine Fisheries Service Coopera- 
tive Shark Tagging Program in each carcass and informa- 
tion reported (by the capturing anglers) to this program 
(Zewinski”). Although the capture-and-handling event was 
not directly observed for either shark, anglers mentioned 
in their testimonials that 1 shark was hooked “deep” and 
was bleeding heavily during hook removal and that the 
other shark seemed physically exhausted and would not 
swim away despite being revived for several minutes. The 
specific cause of these mortalities cannot be determined, 
but their occurrence confirms that mortality can, and does, 
occur in the land-based shark fishery. For this reason, we 
recommend that a PRM rate of 5.6% (the upper limit of 
the 95% CI) be used when estimating mortality for stock 
assessment purposes. 
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Figure 3 


Box plots of the relationship of (A) hook location and (B) handling location with han- 
dling time for sandbar sharks (Carcharhinus plumbeus) captured in Massachusetts in 
2019-2021. The numerals below the plots are the sample sizes for the types of hook and 
handling locations. The thick line within each box indicates the median, the upper and 
lower parts of each box represent the first and third quartiles (the 25th and 75th per- 
centiles), the whiskers extending above and below each box correspond to 1.5 times the 
interquartile range, and the dots represent values outside this range. 


Behavioral recovery from land-based catch and release 
occurred relatively rapidly in sandbar sharks (1.33—17.18 h), 
with 63% of individuals in our study recovering within 
6.36 h after release (i.e., the mean recovery observed over 
all individuals). This mean recovery time is comparable to, 
but also faster than, what has been estimated for blac- 
knose (C. acronotus) (11.7 h, Knotek et al., 2022), blacktip 
(10.5 h, Whitney et al., 2017; 29.6 h, Binstock et al., 2023), 
and bull (8.0 h, Binstock et al., 2023) sharks following 
catch and release in recreational fishing events and for 
sandbar sharks following capture in bottom longline fish- 
ing (11.7 h, Whitney et al., 2021). Although low and 
unequal sample sizes precluded statistical analysis in our 
study, sharks that incurred physical injury from being 
hooked in the gills or that were foul hooked and in poorer 
condition at release had longer recovery times than unin- 
jured sharks (Fig. 5). A similar relationship has also been 
found for tiger sharks following bottom longline capture 
(Whitney et al., 2021). The ability of sandbar sharks to 
rapidly recover from land-based catch and release is likely 
a major factor contributing to the high survivorship 
observed in this fishery. 

Postrelease monitoring periods achieved in our study 
were relatively short (mean: 2.5 d, 60 h) but were suffi- 
cient for monitoring the period during which the majority 
of PRM typically occurs in sharks (Ellis et al., 2017). For 
example, Whitney et al. (2021) documented that all (n=61) 


PRM events for several carcharhinid sharks, including 
sandbar sharks, occurred within 12 h, and Weber et al. 
(2020) reported that all PRM events for blacktip sharks 
caught in a land-based shark fishery occurred within 6 h 
of release, despite their monitoring of animals for weeks 
after release. These trends, combined with the direct 
observation of the mortalities of 2 sandbar sharks within 
24—36 h of release (see the previous mention of necropsies 
in this “Discussion” section), provide strong evidence that 
the durations of ADL package deployments achieved in 
this study were sufficient for the detection of the major- 
ity of mortality events. This notion is also supported by 
reports from several other studies on recreationally caught 
sharks documenting that most or all PRM events occur 
within 24 h (e.g., Heberer et al., 2010; French et al., 2015; 
Whitney et al., 2017; Mohan et al., 2020). Nonetheless, it is 
possible that mortality occurred after monitoring ceased, 
particularly for the few sharks that were hooked in the 
gills or experienced physical trauma (e.g., Kneebone et al., 
2013). Regardless, undocumented, delayed mortalities of 
sharks in poor condition were unlikely to have yielded 
total PRM rates exceeding the upper value of the 95% CI 
(i.e., 5.6%). 

The experience level and attitudes of volunteer anglers 
may have contributed to the high survivorship of sandbar 
sharks monitored during this study. Efforts were made to 
observe anglers of all experience levels, but there was no 


Kneebone et al.: High survivorship of Carcharhinus plumbeus in a land-based fishery 9 


Tailbeat period (s) 
= he) 
on oO 


= 
© 


18 


30 


Postrelease period (h) 


Figure 4 


Estimated periods for postrelease behavioral recovery of 54 sandbar sharks (Carchar- 
hinus plumbeus) for which at least 10 h of acceleration data were available. Sharks 
were captured in Massachusetts from June through September in 2019-2021. Recovery 
was modeled as a function of the tailbeat period (TBP) in seconds. The solid dark gray 
lines represent nonlinear mixed model estimates of TBP, and the vertical dashed gray 
lines indicate when behavioral recovery occurred (i.e., when TBP reached 80% of its 
asymptote). The solid red line depicts the nonlinear mixed model estimate for all sharks 
combined, and the vertical dashed red line represents the predicted recovery time over 
all sharks. The gray dots represent the average 10-min TBP calculated for each shark; 
darker shades of gray indicate overlapping values. 


discernable effect of angler experience on fight time, han- 
dling time, or animal release condition. Furthermore, vol- 
unteer fishermen of all experience levels used similar 
techniques and fishing gear to capture and handle sand- 
bar sharks, and many anglers expressed their intent to 
minimize fight and handling times however possible. Sim- 
ilar acknowledgement of the need to reduce negative 
effects of land-based recreational shark fishing has been 
observed in surveys of land-based shark fishermen in 
Florida (Guay et al., 2021) and Texas (Gibson Banks et al., 
2023). Nevertheless, multiple anglers used J-hooks to cap- 
ture sandbar sharks, a violation of current Massachusetts 
state law. Candid conversations with these anglers 
revealed that some of them were unaware of the regula- 
tion requiring circle hooks and that others assumed the 
hooks they purchased were circle hooks or never transi- 
tioned to circle hooks because of their familiarity with the 
use of J-hooks but not circle hooks. Given that J-hooks 
were more likely to hook sandbar sharks in places other 
than the jaw, improved outreach is needed to better inform 
land-based shark fishermen of fishing regulations and the 
benefits of using circle hooks to reduce physical injury. 


Premature detachment of an ADL package was fre- 
quently encountered during this study (40.3% of packages 
were released more than 25% prematurely). The cause of 
premature release could not be determined in every 
instance; however, the correlation between shark rolling 
behavior and tag detachment provides evidence that some 
of the packages released because the sharks rubbed 
against the bottom. Identifying positional changes of ani- 
mals from acceleration-only data during highly dynamic 
movements can be difficult (Noda et al., 2014), and it was 
particularly challenging to identify the exact instant when 
a tag detached from an animal (i.e., when acceleration 
data transitioned from indicating shark movements to 
indicating the release of a tag). Together, these factors 
introduced the potential for incorrect categorization of 
rolling behavior during tag release. However, the observa- 
tion that premature releases were more common when 
sharks rolled to the left side supports our hypothesis of 
behaviorally driven detachment because rolling to the left 
would bring the tag package toward the seafloor. In the 
context of postrelease survival, the observance of rolling 
behavior in 19 of 27 instances (70.4%) of premature tag 
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Figure 5 

Box plots of the relationship between release condition and 
behavioral recovery period for sandbar sharks (Carcharhi- 
nus plumbeus) captured in Massachusetts during 2019- 
2021. The numerals below the plots are the sample sizes 
for the release conditions. The thick line within each box 
indicates the median, the upper and lower parts of each 
box represent the first and third quartiles (the 25th and 
75th percentiles), the whiskers extending above and below 
each box correspond to 1.5 times the interquartile range, 
and the dots represent values outside this range. 


release is evidence of animal vigor and indicates that the 
individuals in these instances were likely to survive fol- 
lowing premature detachment of the ADL package. 


Conclusions 


The results of this study provide strong evidence that 
land-based catch-and-release fishing for sandbar sharks 
in Massachusetts is a sustainable practice when state 
regulations and best-practice guidelines are followed. 
Although total catch of sandbar sharks of the land-based 
shark fishery in Massachusetts, as well as across all U.S. 
states, is not well quantified, our estimated PRM rate of 
5.6% indicates that incidental mortality in this fishery is 
very low and should minimally affect the ongoing recovery 
of the sandbar shark population in the Atlantic Ocean and 
Gulf of Mexico. Regardless, outreach efforts, particularly 
those targeting new entrants into the land-based fishery, 
should promote the awareness and use of best practices 
and compliance with existing state shark fishing regula- 
tions (e.g., requirements to use circle hooks and to leave 
animals in the water). Future research should be con- 
ducted to estimate PRM of sandbar sharks in land-based 
fisheries occurring in other states or regions where envi- 
ronmental conditions (e.g., warmer air and water tem- 
peratures) and handling practices may be different. 
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Figure 6 


Box plots of the timing of the release of tag packages in rela- 
tion to observations of rolling behavior at the time of tag 
release for sandbar sharks (Carcharhinus plumbeus) cap- 
tured in Massachusetts during 2019-2021. Release error 
represents the deviation from the expected time of release 
(based on the size of the galvanic timed release used), with 
positive values indicating delayed release and negative val- 
ues indicating premature release. The asterisk (*) denotes a 
significant difference in the occurrence of premature release 
following rolls to the left side (tag side). The numerals below 
the plots are the sample sizes for the rolling behaviors. The 
thick line within each box indicates the median, the upper 
and lower parts of each box represent the first and third 
quartiles (the 25th and 75th percentiles), and the whiskers 
extending above and below each box correspond to 1.5 times 
the interquartile range. 


Resumen 


La captura y liberacién de tiburones aleta de cartén 
(Carcharhinus plumbeus) en la pesca recreativa de orilla es 
cada vez mas popular en Massachusetts. En colaboraci6n 
con 21 pescadores voluntarios con distintos niveles de expe- 
riencia, observamos y documentamos los artes de pesca, los 
cebos y las técnicas empleadas para capturar y liberar 67 
tiburones aleta de cartén. El destino de cada tibur6én liber- 
ado se monitore6 utilizando una marca de registro de datos 
de aceleracién (ADL) incrustada en un paquete flotador per- 
sonalizado que se fij6 a la primera aleta dorsal con una lib- 
eracion galvanica temporal. Los 67 paquetes se recuperaron 
tras su desprendimiento después de periodos de seguimiento 
de 0.15-9.98 d. El examen de las series temporales de pro- 
fundidad, periodo de latido de la cola (TBP), cabeceo y bal- 
anceo de 65 ADL que registraron datos revel6é una elevada 
supervivencia de los tiburones marcados, todos ellos vivos 
al momento del desprendimiento de la marca. Se estim6 
que la recuperaci6n del comportamiento se habia producido 
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de 6.36 h en promedio, después de la liberacién, con base 
en las tendencias del TBP de 54 tiburones aleta de cartén 
con al menos 10 h de datos de aceleracién tras la liberaci6n. 
Estos resultados indican que los tiburones aleta de cartén 
son notablemente resistentes a la captura y liberaci6n en la 
pesqueria de orilla de tiburones de Massachusetts. 
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Abstract—The Lebranche mullet 
(Mugil liza) is an economically import- 
ant estuarine species found along the 
coast of Brazil. The southern population 
of this species ranges from the coast 
of the state of SAo Paulo (23°S) to the 
coast of Argentina (36°S). It migrates 
annually among estuaries in Argentina, 
Uruguay, and southern Brazil to repro- 
duce in their marine spawning grounds 
(~26°S). We evaluated whether varia- 
tions in the otolith chemical composition 
for fish acquired at nursery grounds can 
be used to make distinctions among 
4 nursery areas of the southern popula- 
tion. Analysis of otolith microchemistry 
included concentrations of 9 elements: 
lithium (Li), magnesium (Mg), copper, 
zinc, strontium (Sr), cadmium, barium, 
lanthanum, and lead. When using ran- 
dom forest classification algorithms, the 
maximum accuracy of 96% in assign- 
ment of nursery habitat for Lebranche 
mullet between 2 groups (fish caught at 
nurseries in Brazil and fish migrating 
from waters of Uruguay and Argentina) 
was achieved with the combination of Li, 
Mg, and Sr. Our results indicate that the 
elemental composition of otoliths can be 
an important tool for establishing con- 
nectivity between nursery areas used by 
Lebranche mullet. We discuss the impli- 
cations of this result for the structure 
of the population and the management 
of the mullet fishery in southern Brazil 
relative to the limitations of the meth- 
ods we employed. 
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The Lebranche mullet (Mugil liza) is 
one of the most economically import- 
ant marine species in southern Brazil 
(Sant’Ana et al., 2017; dos Santos et al., 
2018). This species ranges from the 
Gulf of Mexico to Argentina (Menezes 
et al., 2015), and the southern popu- 
lation, discriminated on the basis of 
population genetic analyses, is distrib- 
uted from the coast of the state of Sao 
Paulo in Brazil (23°S) to Buenos Aires 
Province in Argentina (36°S) (Mai et al., 
2014; Lemos et al., 2017). The states of 
Rio Grande do Sul and Santa Catarina 
alone account for approximately 
90-95% of the total catch in Brazil 
(MPA and MMA, 2015). In this area in 
Brazil, the average annual catch of com- 
mercial fisheries was 7150 metric tons 
(t), with a maximum catch of 13,600 t in 
2007 (MPA and MMA, 2015; Sant’Ana 
et al., 2017). To the south, in Uruguay, 
the average annual catch was 200 t for 
the period 2000—2007 (Lorenzo, 2010). 
In Argentina, the Lebranche mullet is a 
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target of asmall-scale fishery, mainly off 
the northern coast of the Buenos Aires 
Province, where it is a popular food 
resource and a favorite of sport fisher- 
men (Gonzalez Castro et al., 2009). 

The Lebranche mullet is a marine 
migrant estuarine species (Mai et al., 
2018, 2019). The southern population, 
to reproduce, migrates annually from 
estuaries and lagoons along the coasts 
of Argentina, Uruguay, and southern 
Brazil to offshore waters off the north- 
ern coast of Santa Catarina and the 
southern coast of the state of Parana 
in Brazil (~26°S) (Vieira and Scalabrin, 
1991; Gonzalez Castro et al., 2009; 
Lemos et al., 2014, 2016). This migra- 
tion has been reported to start gradually 
and to follow the temporal movement of 
the range of sea-surface temperatures 
of 19—21°C in a south—north direction 
from April to July (Lemos et al., 2014, 
2016). According to Lemos et al. (2014), 
the broad geographical range of spawn- 


ing locations along the coast of Santa 


Catarina indicates that spawning of Lebranche mullet 
is more closely related to oceanographic conditions than 
to specific locations. These oceanographic conditions are 
salinities of 32-35 and sea-surface temperatures of ~20 °C, 
found at depths of ~50 m in waters on the continental 
shelf (Lemos et al., 2014, 2016). The larvae and early juve- 
niles of this species are neustonic and are carried, mainly 
through wind-driven currents, from spawning grounds 
toward the coast, where they concentrate in surf zones 
prior to recruitment to estuarine environments (Vieira, 
1991; Lemos et al., 2021). 

Many estuaries along the southern coast of Brazil 
and the coasts of Uruguay and Argentina are used by 
Lebranche mullet as nursery grounds (Ramos and Vieira, 
2001; Loebmann and Vieira, 2005; Gonzalez Castro 
et al., 2009; del Favero and Dias, 2015; Callic6 Fortunato 
et al., 2017; Mai et al., 2018), and the Patos Lagoon estu- 
ary is considered the largest nursery area for the south- 
ern population of Lebranche mullet (Lemos et al., 2021). 
The Lebranche mullet has a flexible life history with a 
high degree of behavioral plasticity and wide salinity 
tolerance as juveniles. According to Mai et al. (2019), at 
least 3 patterns of coastal use and recruitment have been 
observed for juveniles. The majority (~90%) of juvenile 
Lebranche mullet migrate toward estuaries before com- 
pleting their first year of life (64%) or after (25%) their 
first year of life (Mai et al., 2019). However, approxi- 
mately 10% of the population does not always use estu- 
aries as nurseries and can complete its entire life cycle in 
the marine environment (Callicé Fortunato et al., 2017; 
Mai et al., 2019). Furthermore, some adults return to an 
estuary after spawning in the ocean, but others never 
return to an estuary after spawning and remain in the 
marine environment (Mai et al., 2019). 

Schroeder et al. (2023), combining otolith shape analy- 
ses and meristic and morphometric analyses of Lebranche 
mullet captured in waters of Brazil, discriminated 2 sub- 
populations. From the same study, results of analyses of 
the deposition of multiple elements in otolith cores from 
Lebranche mullet collected in Brazil and Argentina indi- 
cate the existence of more than one nursery area for this 
species (Schroeder et al., 2023). The identification of nurs- 
ery origin of Lebranche mullet is extremely important and 
may help guide strategies to ensure the efficient manage- 
ment and sustainability of fisheries that target this spe- 
cies (Lotze et al., 2006). 

The use of the chemical composition of otoliths has 
allowed stock identification (Thresher, 1999; Campana 
et al., 2000), descriptions of migration patterns and habitat 
utilization (Secor et al., 1995; Kennedy et al., 2002; Elsdon 
and Gillanders, 2003; Gillanders, 2005), descriptions of the 
life histories of different commercial fish species (Rooker 
et al., 2001; Artetxe-Arrate et al., 2019; Moll et al., 2019; 
Bae and Kim, 2020; Avigliano et al., 2021), and identifica- 
tion of nursery habitats (Gillanders and Kingsford, 2000; 
Vasconcelos et al., 2008; Reis-Santos et al., 2013; Tournois 
et al., 2017). Properties of otoliths allow them to form 
growth increments and retain chemical elements, asso- 
ciated with their metabolically inert capacity, and these 
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signals of growth and chemistry provide a permanent 
record of the environmental conditions (such as the tem- 
peratures and chemical composition of the water) to which 
a fish was exposed (Campana and Neilson, 1985; Campana, 
1999; Panfili et al., 2002). For this reason, researchers have 
long appreciated the importance of otoliths in studies on 
the biology of teleost fish (Campana and Neilson, 1985). 

In this study, we tested the applicability of the elemental 
fingerprint in otoliths of Lebranche mullet to identification 
of nursery grounds utilized by the southern stock. For this 
purpose, the microchemistry of otoliths from Lebranche 
mullet was compared between 4 sampling sites, all nursery 
areas, off the southern coast of Brazil. The concentrations 
of the following chemical elements in otoliths were exam- 
ined: lithium (Li), magnesium (Mg), copper (Cu), zinc (Zn), 
strontium (Sr), cadmium (Cd), barium (Ba), lanthanum 
(La), and lead (Pb). 


Material and methods 
Study area and fish collection 


Adult Lebranche mullet (n=83) were collected directly 
from landings of the artisanal and commercial fisher- 
ies (gill-net fisheries) that operate in coastal waters and 
estuarine areas along the coasts of the southern states of 
Brazil (Fig. 1). Lebranche mullet use as nurseries several 
estuaries and lagoons that are part of the complex hydro- 
graphic system in and along the coastal plain in south- 
ern Brazil (Ramos and Vieira, 2001; Lemos et al., 2014). 
Therefore, we chose to collect specimens caught in 3 estu- 
arine systems used as nursery grounds by juveniles on the 
southern coast of Brazil: Laguna (LAG, 28°S), Tramandai 
(TRA, 29°S), and Patos Lagoon estuary (PLE, 32°S). We 
compared those specimens with individuals in prespawn- 
ing migration that were coming from nurseries in Uru- 
guay and Argentina; they were captured at the border 
between Brazil and Uruguay, near the city of Chui (CHU, 
33°S), by purse seine at the time of northern migration 
(Table 1). Individuals were collected between May and 
October 2011. All fish were transported on ice to the labo- 
ratory and preserved frozen until dissection. All specimens 
were measured for total length (TL) in millimeters and 
total weight (TW) in grams, and the sagittal otoliths were 
extracted, cleaned, dried, and stored for further analysis. 


Otolith preparation and sample analysis 


The method we used for obtaining otolith microchemical 
data follows that of Mai et al. (2018, 2019). Otoliths, one 
from each of the 83 collected specimens, were embedded 
in polyester resin, and transverse sections were cut by 
using an IsoMet Low Speed Saw’ with a diamond blade 
(Buehler, Lake Bluff, IL). Prior to analysis through laser 


1 Mention of trade names or commercial companies is for identi- 
fication purposes only and does not imply endorsement by the 
National Marine Fisheries Service, NOAA. 
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The LA-ICP-MS_ analysis was 
done with a Thermo Scientific iCAP 
@ ICP-MS (Thermo Fisher Scientific 
Inc., Waltham, MA) coupled to a 213-nm 
New Wave neodymium-doped yttrium- 
aluminium-garnet laser (Electro Scien- 
tific Industries Inc., Beaverton, OR) at 
the Instituto de Geociéncias da Univer- 
sidade de Sao Paulo. Using this analy- 
sis, we quantified the concentrations of 
the following isotopes in the otoliths of 
Lebranche mullet: ‘Li, 7°Mg, °Cu, °°Zn, 
8o, AO ip a, 1? Vacand °° Pb, A eri- 
ter with a width of approximately 30 um 
was generated in otoliths with a pulse 
frequency of the laser of 10 Hz, a scan 
speed of 30 um/s, and energy outputs of 
10 and 12 J/cm? per pulse. The ablated 
material was placed through a Teflon- 
coated tube into the ICP-MS by using a 
combination of helium (*He) and argon 
(*°Ar) as a carrier gas. 

Before each otolith ablation, the back- 
ground signal intensity was measured 
with the ICP-MS for 60 s. Counts per 
second for Sr and calcium (Ca) were 
subtracted from the background level, 
and the ratios of the concentrations of 
elements to those of Ca were then calcu- 
lated for all otoliths analyzed. The max- 
imum analyte intensities and minimum 
interferences in the performance of 
the ICP-MS were determined by using 
oxides and double charged ions. The 
equipment was calibrated and tuned 
with the certified reference material 
MACS-3 (U. 8S. Geological Survey): 5.7% 
86Sr and 3.2% '°’Ba. To correct for mass bias and machine 
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Figure 1 


Map of southern Brazil showing the 4 sampling sites where Lebranche mullet 
(Mugil liza) were caught in May—October 2011. Laguna, Tramandai, and Patos 
Lagoon estuaries are estuarine systems used by juvenile Lebranche mullet as 
nursery grounds, and the fourth site is near Chui, a city at the border between 
Brazil and Uruguay. 


ablation inductively coupled plasma mass spectrometry 


(LA-ICP-MS), the sections (0.56 mm) were mounted on 
glass slides (with cyanoacrylate glue), polished (with 
silicon-carbide sanding paper no. 8000), sonicated for 
1 min, rinsed 3 times with ultrapure distilled water 
(Milli-Q, Merck-Millipore), and dried in a laminar flow 
cabinet (Mai et al., 2018, 2019). 


drift, after every 10 otoliths, a certified glass reference 
standard, NIST 612 (National Institute of Standards and 
Technology), was run. 

Otolith transects were ablated from the core to the 
edge to capture the entire life history of the collected fish. 
Knowing that Lebranche mullet spawn at sea (Lemos 


Table 1 


Number and size of Lebranche mullet (Mugil liza) caught at each sampling site on the coast 
of southern Brazil from May through October 2011. n=sample size; SD=standard deviation. 


Total length (mm) Total weight (g) 


Site (latitude) Range Mean (SD) Range Mean (SD) 
377-1538 
271-1260 

506-919 


1145-2322 


749 (300) 
597 (212) 
698 (106) 
1771 (413) 


353-558 
319-540 
360—430 
490-590 


Laguna (28°S) 

Tramandai (29°S) 

Patos Lagoon estuary (32°S) 
Chui (33°S) 


435 (59) 
403 (49) 
402 (20) 
545 (34) 
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et al., 2014, 2017; Mai et al., 2019), the sampling zone 
around the otolith core was chosen manually to represent 
the otolith section formed while the fish was in a nursery 
area, with the ratios of concentrations of Sr and Ba to 
concentrations of Ca in otoliths assumed to be charac- 
teristic of the water chemistry of an estuarine zone (Mai 
et al., 2018, 2019). According to Mai et al. (2019), values 
of 0.02—0.13 mol/mol for ratios of Ba concentration to Ca 
concentration after the core zone indicate that a species 
was living in a nursery area that was an estuarine zone. 
For characterization of the nursery origins of specimens, 
at least 20 chemistry readings (by laser ablation spot) in 
the nursery zone of each otolith were made. 

Data reduction, including background corrections, con- 
version of mass count data to concentrations (in parts per 
million), and limits of detection, was done for each sam- 
ple by using Glitter software (Australia Research Coun- 
cil National Key Centre for Geochemical Evolution and 
Metallogeny of Continents, Macquarie University, Sydney, 
Australia). In further analysis, we considered the concen- 
trations of 9 elements (Li, Mg, Cu, Zn, Sr, Cd, Ba, La, and 
Pb) that were all well above the detection limits calcu- 
lated by using Glitter (0.020, 0.080, 0.025, 0.057, 0.155, 
0.004, 0.007, 0.014, and 0.015 ppm, respectively). Concen- 
trations of trace elements (originally in parts per million) 
were transformed to ratios of concentrations of elements 
to concentrations of Ca, expressed as units of millimoles 
per mole or micromoles per mole (Thorrold et al., 1998; 
Forrester and Swearer, 2002; Swearer et al., 20038; Dorval 
et al., 2005; Brown, 2006). 


Statistical analyses 


We used a nonmetric multidimensional scaling plot to 
visualize the elemental composition of otoliths and poten- 
tial chemical differences among sampling sites from a 
dissimilarity matrix of calculating the Euclidean dis- 
tances. For each nursery area, we analyzed the element 
concentration distributions by using box plots. We visual- 
ized otolith chemical signatures for each nursery area by 
using radar plots, based on the mean proportion of each 
element concentration among the 83 otoliths analyzed. 
Concentrations of all elements were checked for normality 
and homoscedasticity prior to statistical analyses. Para- 
metric assumptions were not met (according to results of 
Shapiro—Wilk and Bartlett tests), even after log,, trans- 
formation; therefore, we used nonparametric statistics. 
To detect whether otolith elemental fingerprints differed 
among nursery areas, we employed the Kruskal—Wallis 
test with Dunn’s post hoc test. Spearman correlation 
analysis was used to examine the correlation between the 
ratios of element concentrations in the otoliths. 

We used the elemental fingerprinting index (EFI) cre- 
ated by Moll et al. (2019) as a measure of the similarity 
among individuals when considering numerous chemical 
signatures in their otoliths. The EFI ranges from 0 to 1, 
with a value of 0 indicating that the elemental composi- 
tions in the otoliths of the 2 compared individuals are most 
different and a value of 1 indicating the highest similarity 


in the elemental compositions in the otoliths of the 2 com- 
pared individuals. We conducted comparisons of element 
concentrations in otoliths for all individuals within one 
sampling area and comparisons for all individuals among 
the different sampling areas. The EFI values were tested 
for significant differences within and between the 4 sam- 
pling areas by using the Welch test and Games—Howell 
post hoc test. | 

To identify the optimal combination of elements for hab- 
itat discrimination, we used a machine learning method 
based on the random forest algorithm developed by 
Breiman (2001). Machine learning methods, such as ran- 
dom forest, have been reported to have greater classifica- 
tion efficiency for nursery identification based on otolith 
microchemistry, especially when data sets include multi- 
ple specific otolith signatures or several chemical elements 
(Mercier et al., 2011). Following this method, 5000 classi- 
fication trees were built, considering all possible combina- 
tions of the 9 elements measured. The analysis proceeded 
in 2 steps. First, a random forest classifier was built by 
using the chemical signatures in otoliths for 75% of indi- 
viduals that were randomly selected from specimens of 
known nursery origin. Second, the otolith signatures from 
the remaining individuals (25%) were used in the random 
forest classifier to predict origin and to test the prediction 
accuracy by measuring the proportion of individuals cor- 
rectly assigned to their known origin (CHU, LAG, TRA, or 
PLE). Additionally, the algorithm was also executed with 
the origin site defined as binary, classified as either the 
coasts of Uruguay and Argentina (CHU) or estuaries of 
Brazil (LAG, TRA, and PLE). Details of the method and 
application are provided in Mercier et al. (2011, 2012). 

All statistical analyses were done in statistical software 
R (vers. 4.2.2; R Core Team, 2022) with the packages dplyr 
(vers. 1.0.10; Wickham et al., 2022), class (vers. 7.3-22; 
Venables and Ripley, 2002), FactoMineR (vers. 2.9; Lé 
et al., 2008), FSA (vers. 0.9.5; Ogle et al., 2023), ggplot2 
(vers. 3.4.4; Wickham, 2016), gtools (vers. 3.9.3; Warnes 
et al., 2022), MASS (vers. 7.3-58.1; Venables and Ripley, 
2002), mda (vers. 0.5-3; Hastie et al., 2022), randomFor- 
est (vers. 4.7-1.1; Liaw and Wiener, 2002), and stats (vers. 
4.2.2: R Core Team, 2022). 


Results 
Otolith microchemistry 


Mean concentrations of elements in examined otoliths 
varied greatly across the 4 sampling sites, CHU, PLE, 
TRA, and LAG (Table 2, Fig. 2). The ratios of concen- 
trations of 5 elements, Li, Cu, La, Pb, and Sr, to concen- 
trations of Ca in the otoliths from specimens differed 
significantly (Kruskal—Wallis pairwise test: P<0.05) 
between 2 or 3 sites. The ratios of the concentrations of Li 
to Ca in otoliths were significantly higher for fish caught 
at CHU than for those caught at the other 3 sampling 
sites. Concentration ratios of La to Ca were lower in oto- 
liths from fish caught at CHU than in those from fish 
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Table 2 


Mean ratios of the chemical concentrations of 9 elements to concentrations of calcium in the nursery zone in otoliths 
from Lebranche mullet (Mugil liza) caught at 4 sampling sites in southern Brazil, Chui (CHU), Patos Lagoon estuary 
(PLE), Tramandai (TRA), and Laguna (LAG), between May and October 2011. The elements are lithium (Li), magne- 
sium (Mg), copper (Cu), zinc (Zn), strontium (Sr), cadmium (Cd), barium (Ba), lanthanum (La), and lead (Pb). Signifi- 
cant differences in chemical signatures in otoliths between sites were examined by using Kruskal—Wallis and Dunn’s 
post hoc tests for each element with the significance level of 0.05. An asterisk (*) indicates that signatures were 
revealed in Kruskal—Wallis tests to be significantly different between sites. The results of Dunn’s test for pairwise 
comparisons are provided. Standard deviations (SDs) of mean ratios are given in parentheses. NS=not significant. 


Element 


Li 
Mg 
Cu 
Zn 
Cd 
Ba 
La 
Pb 
Sr 


CHU 


4.57 (1.62) 
3.20 (0.87) 
0.22 (0.20) 
0.43 (0.78) 
2.03 (1.75) 
0.07 (0.05) 
1.22 (1.29) 
0.21 (0.27) 
0.43 (0.10) 


Mean ratio (SD) 


PLE 


1.46 (1.16) 
3.47 (1.56) 
0.30 (0.44) 
0.60 (0.58) 
2.25 (1.81) 
0.08 (0.07) 
1.81 (1.62) 
0.30 (0.21) 
0.45 (0.12) 


TRA 


1.91 (1.06) 
3.17 (1.20) 
0.22 (0.14) 
0.67 (1.35) 
2.57 (1.85) 
0.06 (0.05) 
L-73A1,19) 
0.47 (0.40) 
0.46 (0.09) 


LAG 


2.01 (1.11) 
3.07 (1.29) 
0.20 (0.24) 
0.59 (0.53) 
2.22 (1.82) 
0.07 (0.05) 
1.44 (0.94) 
0.34 (0.33) 
0.48 (0.10) 


Dunn’s test 


CHU>PLE, TRA, LAG 
All NS 

PLE>LAG 

All NS 

AU NS 

All NS 

CHU<PLE, TRA 
TRA>CHU, PLE 
LAG>CHU, TRA 


1? 


caught at PLE and TRA. The ratios of concentrations of 
Pb to Ca were higher in otoliths from specimens caught 
at TRA than in the otoliths of those caught at CHU and 
PLE. The concentration ratios of Cu to Ca were signifi- 
cantly higher for otoliths from fish caught at PLE than 
for otoliths from those caught at LAG. Ratios of the con- 
centrations of Sr to Ca in otoliths was higher for speci- 
mens caught in LAG than for those from CHU and TRA, 
and ratios of the concentrations of Mg, Zn, Cd, and Ba to 
Ca in otoliths were similar between sites (Table 2). 

Results from Spearman correlation analysis (Table 3) 
indicate that some of the correlations between ratios of 
element concentrations were significant (P<0.05). The 
correlations between the concentration ratios of Cu, Zn, 
Cd, La, and Pb to Ca indicate a significantly positive rela- 
tionship. The only significant negative correlation was 
between the ratio of the concentration of Li to Ca and the 
ratio of the concentration of La to Ca. 


Chemical signatures in the nursery zone of otoliths 


Examination of the data set on the chemical composition 
of the nursery zone of otoliths included habitat-specific 
plots of elemental fingerprints with similarities and dif- 
ferences in proportions of element concentrations in oto- 
liths (Fig. 3). The nonmetric multidimensional scaling plot 
indicates differences in chemical signatures in otoliths for 
each sampling site. Element concentrations in otoliths of 
specimens from CHU had lower variability among indi- 
viduals than those in otoliths of fish from LAG, TRA and 
PLE, and they are clustered on the left side of the ordi- 
nation plot. Elemental compositions in the nursery zone 
of otoliths from specimens caught at PLE, LAG, and TRA 
were similar and clustered on the right side of the plot, 


indicating more similarity in otolith compositions among 
these sites than between these sites and CHU (Fig. 4). 

Results from the pairwise comparisons of multiple ele- 
ment concentrations for individuals within a single sam- 
pling site, made by using the EF', indicate that the otolith 
compositions for individuals collected at CHU are signifi- 
cantly more similar to those of each other (higher EFI 
value) and are significantly different from those of speci- 
mens caught at the other 3 sites (P<0.001) (Table 4). Com- 
parisons between sampling sites also identified similarities 
in element concentrations of otoliths from fish caught at 
PLE, TRA, and LAG (Table 4). The EFI values for compar- 
isons of chemical signatures in otoliths of fish caught at 
CHU to those of fish caught at each of the other 3 sites, 
LAG, PLE, and TRA, indicate that similarities in element 
concentrations between otoliths within these groups of fish 
were significantly lower (lower EFI value) than the similar- 
ities between otoliths from fish collected at CHU (P<0.001). 
The EFI values from comparisons between 2 sites at a time 
for all 4 nurseries are smaller than the EFI values from 
within-site comparisons, indicating that element concen- 
trations in otoliths of specimens are more similar within 
sites than between all 4 nurseries (Table 4). 


Discrimination of nursery grounds 


In using the random forest method to compare the 4 nurs- 
ery grounds sampled, we built classification trees by con- 
sidering multiple combinations of the 9 chemical elements 
tested. We observed that the addition of elements gener- 
ally improved the accuracy of the random forest classifier, 
increasing the average levels of accuracy in discriminat- 
ing fish nursery origin of Lebranche mullet (Fig. 5A). The 
maximum accuracy (55.1%) in assignment of nursery 
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Figure 2 


Mean ratios of the chemical concentrations of 9 elements to the concentrations of 
calcium in the nursery zone of otoliths of Lebranche mullet (Mugil liza) collected 
from May through October 2011 at 4 sampling sites in southern Brazil: Chui 
(CHU), Patos Lagoon estuary (PLE), Tramandai (TRA), and Laguna (LAG). The 
elements are cadmium (Cd), lanthanum (La), copper (Cu), barium (Ba), stron- 
tium (Sr), lead (Pb), lithium (Li), magnesium (Mg), and zinc (Zn). In each box 
plot, the thick black line is the median. The upper and lower parts of each box 
represent the first and third quartiles (the 25th and 75th percentiles). The whis- 
kers extend above and below the box no more than 1.5 times the interquartile 
range. Data points that appear beyond the end of the whiskers are outliers. 


origins to the 4 sampling sites was reached with the com- 
bination of the 4 elements Cu, Li, Pb, and Sr, specifically 
the ratios of their concentrations to concentrations of Ca. 
On the basis of evidence that element concentrations in 
otoliths of specimens from CHU were significantly differ- 
ent from those in otoliths of fish from the remaining sites 
sampled and that the nursery habitats for specimens from 


CHULAG PLE TRA 


PLE, TRA, and LAG could not be discrim- 
inated (Figs. 3 and 4, Table 4), we com- 
pared data for specimens from CHU to 
data for fish from the other 3 sampling 
sites combined (Fig. 5B). The maximum 
accuracy of 96% in assignment of nursery 
area was achieved with the combination 
of values for only 3 elements: the ratios 
of the concentrations of Li, Mg, and Sr to 
concentrations of Ca. 


Discussion 


Otolith chemical composition is not a 
direct representation of ambient water 
chemistry, and the relationship between 
environmental and physiological pro- 
cesses is still poorly understood (Thresher, 
1999). Nonetheless, it is widely accepted 
that the mechanisms controlling the 
uptake and incorporation of elements 
into otoliths are mediated by both exoge- 
nous (e.g., salinity) and endogenous (e.g., 
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growth) factors (Walther et al., 2010). 
Not all elements behave similarly in oto- 
liths (Vrdoljak, 2020). Many studies have 
attempted to determine the causes of vari- 
ations in otolith microchemistry (Elsdon 
and Gillanders, 2003; Sturrock et al., 2014, 
2015), and knowledge of the factors that 
influence otolith microchemistry is far 
from complete. However, several scientists 
have claimed that it is not necessary to 
understand all these mechanisms to be 
able to use otolith microchemistry as a tool 
in fisheries management (Thorrold et al., 
1998; Thresher, 1999; Elsdon et al., 2008). 
As expected, results from the analy- 
Sis in our study of element concentra- 
tions in the nursery zone of otoliths from 
Lebranche mullet indicate that the ratios 
of concentrations of Li, Sr, La, Pb, and Cu 
to concentrations of Ca were significantly 
different among the sampling sites. The 
concentrations of certain elements var- 
ied between nursery grounds in south- 
ern Brazil and the southernmost nursery 
areas in Uruguay and Argentina. 
Signatures of Sr and Ba can be used to 
distinguish residences in waters of dif- 
ferent salinities (Walther and Thorrold, 
2006; Elsdon et al., 2008) and have been used to success- 
fully reconstruct environmental histories and coast—-estuary 
movement for Lebranche mullet in South America (Callic6 
Fortunato et al., 2017; Mai et al., 2018, 2019). We found no 
significant difference between sites in the ratio of concentra- 
tions of Ba to Ca in the nursery zone of otoliths from spec- 
imens of Lebranche mullet. This finding indicates that all 
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Table 3 


Spearman correlation coefficients for the relationships between ratios of the concentra- 
tions of 9 elements, lithium (Li), magnesium (Mg), copper (Cu), zinc (Zn), strontium (Sr), 
cadmium (Cd), barium (Ba), lanthanum (La), and lead (Pb), to concentrations of calcium 
in the nursery zone of otoliths from Lebranche mullet (Wugil liza) collected from 4 sam- 
pling sites in southern Brazil from May through October 2011. An asterisk (*) indicates 
a significant correlation between concentrations of elements (P<0.05). 


Ratio of element to calcium 


Zn Cd Ba 


0.56* 

0.10 —0.04 

0.56* 0.47* 0.09 

0.67* 0.62* 0.07 0.55* 

0.24 0.13 -0.11 0.00 Grl7 


sampled individuals developed in estuaries. However, we 
found significantly higher ratios of concentrations of Sr to 
Ca in otoliths of Lebranche mullet from LAG, a trend that 
may be associated with differences in ambient chemistry 
and physical factors. The other sampled nursery areas are 
affected by large drainage basins and, therefore, are heav- 
ily influenced by freshwater discharges; in comparison, LAG 
is a saline wedge estuary located at 28°S and is under the 
influence of higher temperatures than those at other sites 
(D’Aquino et al., 2011). 

Fish absorb heavy metals (e.g., Pb and Cu) from their 
environments during their lives, and these elements are 
preserved within otolith structures (Reis-Santos et al., 2013; 
Sturrock et al., 2014). The accumulation of trace metals in 
fish otoliths depends on a number of factors, including the 
concentrations of elements in the water environment (Ran- 
aldi and Gagnon, 2010). We found higher concentration 
ratios of Pb and Cu to Ca in otoliths from Lebranche mullet 
caught in the estuary at TRA and in the PLE, respectively. 
This result may be an indication of anthropogenic effects in 
estuarine environments. In fact, the average concentrations 
of metals in estuarine water in Patos Lagoon correspond to 
the natural base levels described for other estuarine sys- 
tems, but sporadic increases in Cu and Pb concentrations 
have been related to the ingress of industrial effluents and 
activities in the limnic area of Patos Lagoon (Seeliger et al., 
1997). Although findings from some other research indicate 
associations of the highest Pb concentrations in otoliths 
with port areas (Cuevas et al., 2019), the effects of the time 
during which fish are exposed to Pb and of physiological 
processes should be considered (Vrdoljak, 2020). The nurs- 
eries sampled at both TRA and PLE are in highly popu- 
lated areas. Although concentrations of trace metals in the 
waters at these sites are low, they could be used to identify 
the contribution of polluted nurseries to the stock of Leb- 
ranche mullet in southern Brazil. 


We found significantly higher ratios of the concentra- 
tions of Li to Ca in the nursery zone of the otoliths of 
Lebranche mullet sampled from the southern migratory 
contingent (at CHU) during the time of northern migra- 
tion than in those of fish caught at other sampling sites, 
all estuaries in southern Brazil. In fact, Li was the most 
significant element in the use of otolith microchemistry 
to discriminate between sampling sites (Table 2, Fig. 2). 
According to Sturrock et al. (2012), among the various ele- 
ments found in otoliths, Li has the most positive relation- 
ship with the environment and is one of the elements that 
are reliable geographic markers in marine species. Lith- 
ium is highly bioavailable; however, little attention has 
been given to determining the environmental and physio- 
logical factors that affect incorporation of Li into otoliths 
(Brown, 2006). Furthermore, information about Li in the 
aquatic environment is limited (Tkatcheva et al., 2015). 
This element seems to be associated with estuaries in 
the southern regions of South America, and results from 
other studies already indicate the power of Li in discrim- 
inating habitats for species collected in Brazil, Uruguay, 
and Argentina: the Brazilian codling (Urophycis brasil- 
tensis) (Biolé et al., 2019) and the whitemouth croaker 
(Micropogonias furnieri) (Avigliano et al., 2021). 

Through none of our analyses were we able to use the 
nursery zone of otoliths to differentiate between groups of 
fish caught in each of the estuaries in southern Brazil 
(PLE, TRA, and LAG). However, otoliths of fish caught 
from CHU were different from those of the combined 
group of fish caught in PLE, TRA, and LAG. The similar- 
ity between the chemical signatures in the nursery zone 
of otoliths of Lebranche mullet from the estuaries of Bra- 
zil can be explained by the similar physicochemical condi- 
tions of the waters in those estuaries; after all, water 
chemistry influences the chemical composition of otoliths 
(Kerr and Campana, 2014). According to Adelir-Alves 
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All sites 


Figure 3 


Radar plots of the proportions of element concentrations in the nursery zone 
of otoliths from Lebranche mullet (Mugil liza) collected at 4 sampling sites in 
southern Brazil in May—October 2011: Chui (CHU, purple radii), Patos Lagoon 
estuary (PLE, red radii), Tramandai (TRA, green radii), and Laguna (LAG, 
blue radii). The proportions (shown as percentages) are for concentrations of 9 
elements, lithium (Li), magnesium (Mg), lead (Pb), copper (Cu), barium (Ba), 


zinc (Zn), strontium (Sr), cadmium (Cd), and lanthanum (La). 


et al. (2019), otolith microchemistry analyses provide 
information for a given time and on a fine spatial scale, 
but the usefulness of element concentrations in otoliths 
as natural biogeochemical markers depends on the exis- 
tence of measurable differences in otolith microchemistry 
on a geographic scale relevant to the life history of a spe- 
cies. In fact, the concentration of some elements in oto- 
liths differed significantly between the sampling sites in 


the estuaries of southern Brazil, but the 
general signature in otoliths of fish from 
each habitat had a relatively low level of 
accuracy (55%) in the analysis done 
with the random forest method (Fig. 5A). 
According to Mercier et al. (2011), an 
acceptable degree of accuracy should be 
greater than 75%. Only after grouping 
the specimens from the estuaries in 
Brazil and comparing the element con- 
centrations in their otoliths with those 
of specimens from CHU did the accu- 
racy of nursery assignment increase to 
96% (Fig. 5B). 

Temporal variability in otolith micro- 
chemistry must be considered because 
otolith elemental fingerprints are 
unlikely to remain consistent over 
time (Gillanders, 2002; Hamer et al., 
2003; Clarke et al., 2009; Walther and 
Thorrold, 2009). This temporal variability 
can confuse spatial differences, meaning 
differences in elemental composition of 
otoliths among fish from different nurser- 
ies. Therefore, it is important to consider 
interannual differences in otolith micro- 
chemistry (Gillanders, 2002; Cuveliers 
et al., 2010; Reis-Santos et al., 2012), and 
assignment of origin (to specific nurseries 
or sampling sites) should be done only for 
year-classes for which an otolith chemical 
signature has been identified for juve- 
niles. For this reason, it would be ideal if 
sampling and monitoring programs sys- 
tematically collect otoliths over a period 
of several years, compiling a library of 
reference elemental markers than can be 
later used to determine the nursery ori- 
gin of coastal adult fish (Gillanders, 2002; 
Reis-Santos et al., 2013). 

We assume that individuals collected 
at the southernmost sampling site, CHU, 
were from estuaries in Argentina and 
Uruguay, a notion that is reinforced by 
the distinction in otolith microchemistry 
for fish from this site and the low stan- 
dard deviations of ratios of element con- 
centrations. Several reports on previous 
scientific research have revealed that 
Lebranche mullet migrate from Argen- 
tina and Uruguay to waters in Brazil in 
order to reproduce (Gonzalez Castro et al., 2009; Lemos 
et al., 2014; Mai et al., 2014; Lemos et al., 2016). This 
common finding reflects the local ecological knowledge 
of Brazilian fishermen regarding the existence of various 
morphotypes of Lebranche mullet with different geograph- 
ical origins in the southwestern Atlantic Ocean (Herbst 
and Hanazaki, 2014). Distinct morphological characteris- 
tics are acquired in natural conditions of nurseries where 
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Lebranche mullet of these morphotypes 
sa grow for up to 4-5 years until they 
- 2D stress: 0.124 ™ migrate to reproduce in coastal zones 
(~26°S) (Garbin et al., 2014; Lemos et al., 
2014, 2016). Schroeder et al. (2023) rec- 
ognized 2 subpopulations of Lebranche 
mullet: one subpopulation is associated 
with the central stock in Brazil, and 
another subpopulation is associated 
with the southern stock that comprises 
a mix of individuals from Brazil and 
Argentina. Locations or regions where 
Lebranche mullet are known to spawn 
are more associated to oceanographic 
conditions than to specific geographi- 
cal locations (Lemos et al., 2014, 2016). 
The geographic breadth of the spawning 
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Coordinate 1 does the annual contribution of individ- 

Site @ CHU @ LAG @ PLE @ TRA uals to the different subpopulations or 
contingents. 


Figure 4 , 
The mullet fishery in southeastern 


Nonmetric multidimensional scaling plot of element concentrations in the : 
Are eA : Brazil targets a stock that has been 


nursery zone of 83 otoliths from Lebranche mullet (Mugil liza) caught from . 
May through October 2011 at 4 sampling sites in southern Brazil: Chui (CHU), reported to be overfished in recent decades 
Laguna (LAG), Patos Lagoon estuary (PLE), and Tramandai (TRA). The plot is (de Miranda et al., 2011; MPA and MMA, 
based on a Euclidean distance matrix. 2015). This exploitation is further aggra- 
vated by the fact that Lebranche mullet 
are fished more intensively during the 
reproductive period of this species, especially when shoals 


Table 4 are formed and gradually thickened as they migrate north 
from the nursery areas in Argentina, Uruguay, and south- 
Elemental fingerprinting indices (EFIs) for each sampling site ern Brazil (MPA and MMA, 2015; Lemos et al., 2016). 


(upper table) and between 2 sites (lower table), based on com- 


eee he 2 f 
parisons of chemical signatures in otoliths from Lebranche rovingithabthe 2.smbpopulahionsidesenibed:by chunedet 


mullet (Mugil liza) collected in May—October 2011 in southern at al. (2023) mix with each other during the reproductive 
Brazil. The EFI ranges from 0 to 1, with a value of 0 indicating period and determining in which proportions this mixing 
that the elemental compositions in the otoliths of 2 fish are occurs were the objectives of our study. Our work did not 
most different and a value of 1 indicating that the elemental incorporate analysis of individuals from the same cohort 
compositions in the otoliths of 2 specimens have the highest that were born in different locations; therefore, it was not 
similarity (Moll et al., 2019). An asterisk (*) indicates a signif- possible to estimate the actual contributions of specific 
icant difference between EFT values (P<0.05). Standard devi- nurseries to the adults of the southern population caught 
ations (SDs) of mean EFI values are provided in parentheses. by commercial fisheries during the reproductive period of 
this species. However, we were able to conclude that the 2 


2 ey ; j 
SS eae hone Beto BUTE ae subpopulations do mix during the reproductive period, and 


No. of comparisons Mean EFI (SD) we have established that it is possible to estimate the rel- 
ative contributions of different estuarine nurseries to the 

78 0.89 (0.04)* adult population of Lebranche mullet by using random for- 

276 0.86 (0.05) est classifiers generated with fingerprints of a combination 


231 0.86 (0.04) 


210 0.85 (0.05) of multiple elements in otoliths. 


Comparisons between 2 sites 


Conclusions 

Site No.ofcomparisons Mean EFI (SD) P 

PLE aY 0.81 (0.06)* By The results of our study indicate that the elemental com- 
TRA 208 0.82 (0.06)* <0.001 position of otoliths from Lebranche mullet differed signifi- 
LAG 195 0.83 (0.05)* <0.001 cantly between groups of nurseries and can potentially be 
TRA 297 0.86 (0.04) 0.960 used to estimate the relative contributions of nurseries to 
LAG 294 0.84 (0.06) 0.203 adults in the southern population. Using random forest 
LAG 252 0.87 (0.04) 0.074 classification algorithms with the chemical signatures of 3 
elements (Li, Mg, and Sr) combined, we achieved the 
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Maximum accuracy=96% 
Best combination::Li, Mg, and Sr 


Maximum accuracy=55% 
Best combination: Cu, Li, Pb, and Sr 
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Figure 5 


Accuracy in assigning nursery origins by using microchemistry in otoliths from Lebranche mullet (Mugil liza) caught 
between May and October 2011 at 4 sampling sites in southern Brazil, (A) when the origin was defined in a random 
forest algorithm as the site Chui, Patos Lagoon estuary, Tramandai, or Laguna and (B) when the origin was defined 
as either Chui or the 3 nursery sites combined (Patos Lagoon estuary, Tramandai, and Laguna). The classification 
trees of the algorithm were based on different combinations of up to 9 elements in otoliths. The polygon in each panel 
indicates the range of accuracy values. Error bars indicate standard deviations of mean accuracy levels. The dashed 


lines indicate maximum model accuracy. Cu=copper; Li=lithium; Mg=magnesium; Pb=lead; Sr=strontium. 


maximum accuracy of 96% in identification of nursery ori- 
gin for Lebranche mullet when assigning origin between 
CHU, where individuals of the southern migratory contin- 
gent were caught, and the nurseries in southern Brazil 
(PLE, TRA, and LAG combined). Given that analyses were 
done with individuals from the same population, our find- 
ings corroborate the existence of 2 subpopulations of Leb- 
ranche mullet off the coast of southern Brazil on the basis 
of their use of different nursery areas. Our results also 
confirm that signatures based on the concentrations of a 
combination of multiple elements in otoliths provide infor- 
mation valuable in establishing connectivity between 
nursery grounds and adult reproductive aggregation. 
Therefore, such analysis based on otolith microchemistry 
can be used to aid in the management of mullet fishing. 


Resumen 


El lisa Lebrancha (Mugil liza) es una especie estuarina de 
importancia econdmica que se encuentra a lo largo de la 
costa de Brasil. La poblacién meridional de esta especie se 
extiende desde la costa del estado de Sao Paulo (23°S) hasta 
la costa de Argentina (36°S). Migra anualmente entre estu- 
arios de Argentina, Uruguay y el sur de Brasil para repro- 
ducirse en sus zonas de desove marinas (~26°S). Evaluamos 
si las variaciones en la composicién quimica de los otolitos 
de los peces colectados en las areas de crianza pueden uti- 
lizarse para distinguir entre las 4 areas de crianza de la 
poblaci6n meridional. El andlisis de la microquimica de los 
otolitos incluy6é concentraciones de 9 elementos: litio (Li), 
magnesio (Mg), cobre, zinc, estroncio (Sr), cadmio, bario, 
lantano y plomo. Al utilizar algoritmos de clasificacién de 
bosque aleatorio, se alcanz6 una maxima precision del 96% 


en la asignacion del habitat de crianza de la lisa Lebrancha 
entre 2 grupos (peces capturados en criaderos de Brasil y 
peces migrando desde aguas de Uruguay y Argentina) con 
la combinacién de Li, Mg y Sr. Nuestros resultados indican 
que la composicién de elementos de los otolitos puede ser 
una herramienta importante para establecer la conectivi- 
dad entre las areas de crianza utilizadas por la lisa Lebran- 
cha. Discutimos las implicaciones de este resultado para la 
estructura de la poblacién y el manejo de la pesqueria de 
lisa en el sur de Brasil en relaci6n con las limitaciones de 
los métodos que empleamos. 
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Abstract—Dolphinfish (Coryphaena 
hippurus) are caught throughout the 
western Atlantic Ocean over varying 
spatial and temporal scales. Prior 
attempts to quantify the population 
dynamics of dolphinfish in this region 
have been inhibited by an inability to 
model the spatiotemporal dynamics 
of this stock. We fit a seasonal vector 
autoregressive spatiotemporal (VAST) 
model to quantify the spatiotemporal 
dynamics of western Atlantic dolphin- 
fish, to estimate standardized relative 
indices of abundance during 1986- 
2022 at regional scales, and to estimate 
changes in spatial distribution. The 
magnitude of abundance was greatest 
during spring and summer in north- 
ern spatial strata and was comparable 
over seasons in southern spatial strata. 
Abundance of dolphinfish appeared to 
be stable during 1986-2018 and then 
declined during 2019-2022. This trend 
occurred in all regions, except for in 
Atlantic waters from Cape Hatteras, 
North Carolina, to the southern border 
of Georgia, where abundance remained 
stable during 2019-2022. No shift in 
the distribution of the population was 
detected, but regional patterns of abun- 
dance provide insight into changes in 
the timing of availability. This study 
resulted in the first standardized index 
of relative abundance to capture the 
spatiotemporal dynamics of western 
Atlantic dolphinfish. These results have 
increased our understanding of the pop- 
ulation dynamics of this species in this 
region and should prove useful in future 
attempts to manage the population at 
different spatial and temporal scales. 
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The dolphinfish (Coryphaena hippu- 
rus) is a mid-level trophic predator 
and migratory, cosmopolitan species of 
pelagic fish that is found throughout 
tropical and subtropical regions (Palko 
et al., 1982; Luckhurst, 2017; Molt6 
et al., 2020). The stock structure of dol- 
phinfish caught in the western Atlan- 
tic Ocean, specifically in the southern 
extent of Food and Agriculture Orga- 
nization of the United Nations (FAO) 
major fishing area 21 (the northwest- 
ern Atlantic Ocean) and FAO major 
fishing area 31 (the western Central 
Atlantic Ocean), is not precisely known 
(Damiano, 2023). Using analyses of 
size composition in catches over space, 
life history traits, and allelic frequen- 
cies, Oxenford and Hunte (1986) found 
evidence of a northern and southern 
stock in the western Central Atlantic 
Ocean and hypothesized that there are 
2 unique seasonal migration patterns. 
The northern stock migration has been 
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defined by seasonal movements into 
waters of Puerto Rico during late fall 
and early winter, movement along the 
Greater Antilles and Bahamas and 
into Florida waters during the spring, 
then into waters of the eastern United 
States during late spring and summer, 
and finally outward toward Bermuda 
during late summer (Oxenford and 
Hunte, 1986; Oxenford, 1999). 

Tagging experiments conducted 
by Merten et al. (2014a, 2014b, 2016) 
revealed further complexity in the 
stock connectivity between dolphin- 
fish caught along the eastern United 
States, in the Caribbean Islands, and 
in the Caribbean Sea, expanding the 
northeastern circuit to include coastal 
waters along the northeastern United 
States, and dolphinfish caught along an 
additional pathway to waters in Florida 
from the Caribbean Sea and through 
the Straits of Florida (Damiano, 2023). 
The connectivity between populations 
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of dolphinfish in the western Atlantic Ocean and the Gulf 
of Mexico remains unclear (Damiano, 2023). Results from 
the DNA sequencing conducted by Merten et al. (2015) 
indicate shallow genetic differences among dolphinfish 
caught along the eastern United States and parts of the 
Caribbean Sea, and although the DNA sequencing results 
do not present conclusive evidence of a single stock, mul- 
tiple subpopulations have been identified at smaller geo- 
graphic scales in other ocean basins (Duarte-Neto et al., 
2008; Ochoa-Zavala et al., 2022). Strong genetic sepa- 
ration of populations is likely driven by oceanographic 
barriers (Maggio et al., 2019), of which there are none 
known for the regions composing the western Atlantic 
Ocean (Damiano, 2023). However, despite this uncer- 
tainty, sufficient evidence exists of stock mixing and 
movement linkages in dolphinfish caught in the eastern 
United States and parts of the Caribbean Sea (Merten 
et al., 2014a, 2015, 2016), and that evidence indicates the 
need to model the spatiotemporal dynamics of a western 
Atlantic stock (Damiano, 2023). 

According to a recent review of the condition of interna- 
tional dolphinfish fisheries in the western Central Atlantic 
Ocean that builds on the work of Mahon (1999), dolphin- 
fish continue to support numerous high-value commercial, 
recreational, and artisanal fisheries at different seasonal 
and spatial scales (Merten et al., 2022a). Management 
of dolphinfish in the western Atlantic Ocean falls under 
the purview of the International Commission for the Con- 
servation of Atlantic Tunas, but landings of this species 
are labeled as those of other fishes along with commer- 
cial bycatch taxa and therefore are not actively managed 
(Merten et al., 2022a). 

In the eastern United States, the dolphinfish is man- 
aged under the Dolphin and Wahoo Fishery Management 
Plan by the South Atlantic Fishery Management Council 
(SAFMC, 2003). The annual catch limit for this species in 
the eastern United States is determined by using an aver- 
age of past catches; this approach is a data-limited man- 
agement procedure employed when biological reference 
points cannot be reliably estimated from an assessment 
model (Prager'). At the state level, effort-reducing and 
size-based management measures, such as bag limits and 
minimum size limits, are implemented for recreational 
and for-hire fisheries, but not all states implement both 
types of measures (SAFMC, 2003; Damiano, 2023). 

The Caribbean Fishery Management Council, also a 
U.S. management body, and the Western Central Atlan- 
tic Fishery Commission, established by the FAO, do not 
have fishery management plans for dolphinfish, but 
they provide some monitoring and guidance for sus- 
tainable management of dolphinfish (Damiano, 2023). 
Practical management in Caribbean island waters and 
the Caribbean Sea falls to individual island nations and 


1 Prager, M. 2000. Exploratory assessment of dolphinfish, Cory- 
phaena hippurus, based on U.S. landings from the Atlantic 
Ocean and Gulf of Mexico, 18 p. Southeast Fish. Sci. Cent., Natl. 
Mar. Fish. Serv., NOAA, Beaufort, NC. [Available from website. | 


countries and territories of the Americas (Mahon, 1999; 
Merten et al., 2022a). 

Such a lack of spatially consistent management mea- 
sures 1S uncommon for transboundary stocks for which 
management and governance normally must adhere to an 
international policy (Levesque, 2008). Despite advances 
in spatiotemporal models and enhanced knowledge of the 
stock structure and movement connectivity of dolphinfish 
in the western Atlantic Ocean, no indices of relative abun- 
dance and no stock assessments are based on explicit mod- 
eling of the spatiotemporal dynamics of the population or 
subpopulations that make up this stock (Damiano, 2023). 
There are several other concerns regarding the western 
Atlantic stock. Lynch et al. (2018) reported a negative 
trend in estimated abundance of dolphinfish throughout 
the western Central Atlantic Ocean on the basis of catch 
per unit of effort (CPUE) data from the U.S. pelagic long- 
line (U.S. PLL) logbook for the period 1987—2013. In addi- 
tion, 23 nations still do not report landings of dolphinfish 
to the FAO (Merten et al., 2022a), and increases in sea- 
surface temperatures (SST's) due to climate change may 
affect stock productivity and distribution (Barange et al., 
2018; Damiano, 2023). 

As is the case for highly migratory species, such as the 
swordfish (Xiphias gladius) and tuna species, catch and 
effort information come primarily from fishery-dependent 
sources (Forrestal et al., 2019), and apart from tagging pro- 
grams, there are no fishery-independent data collection or 
monitoring programs for dolphinfish. Sources of fishery- 
dependent data, such as landings and CPUE data, also 
have mostly been disparate over both space and time for 
this species (Rose and Hassler, 1968; Mahon, 1999; Arocha 
et al., 2021; Merten et al., 2022a; Damiano, 2023). Conse- 
quently, past efforts to model the population dynamics of 
dolphinfish in the western Atlantic Ocean have generally 
been accomplished by using a surplus production model 
that has been fit to an index of relative abundance derived 
from fishery-dependent CPUE information for a discrete 
subsection of the stock’s distribution (Damiano, 2023). 
Each of these cases has involved fitting standardization 
models to CPUE data (Prager!; Parker et al.”; Kleisner, 
2008) in order to develop model-based indices that address 
the commonly violated assumption of constant catchabil- 
ity (Maunder and Punt, 2004), then fitting a surplus pro- 
duction model. 

In each of these efforts, CPUE data were not suffi- 
ciently informative because of short time series, “one-way 
trip” patterns, or other common features of fishery- 
dependent data (Maunder et al., 2006; Magnusson and 
Hilborn, 2007), and their use resulted in highly uncertain 
estimates of biological reference points that could not be 


2 Parker, C., K. M. Kleisner, and J. S. Nowlis. 2006. Preliminary 
assessment of the western Central Atlantic dolphinfish (Cory- 
phaena hippurus) stock: a Caribbean regional fisheries mecha- 
nism project, 13 p. Sustainable Fisheries Division Contribution 
SFD-2005-041. [Available from Sustain. Fish. Div., Southeast 
Fish. Sci. Cent., Natl. Mar. Fish. Serv., NOAA, 75 Virginia Beach 
Dr., Miami, FL 33149.] 
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used to inform fisheries management. Kleisner (2008) 
accurately concluded that for assessment of dolphinfish 
to be successful, estimating robust indices of relative 
abundance must be employed with new techniques that 
account for spatiotemporal structure in observed catch 
data (Damiano, 2023). 

The spatiotemporal dynamics of dolphinfish in the west- 
ern Atlantic Ocean are also likely linked to SST (Kleisner, 
2008), given that dolphinfish reproduce year-round within 
their optimal thermal habitat (Schwenke and Buckel, 
2008; Merten et al., 2015; Schlenker et al., 2021) and that 
high catch rates are associated with warm fronts with 
SSTs near 24—25C° (Farrell et al., 2014) and with mats of 
Sargassum macroalgae (Farrell et al., 2014; Rudershau- 
sen et al., 2019). Lynch et al. (2012) developed a 2-stage 
delta-generalized linear mixed model (GLMM) capable of 
modeling temperature regimes and habitat use of pelagic 
fish. They applied this method by fitting a delta-GLMM 
to CPUE data for highly migratory species, including dol- 
phinfish, from the U.S. PLL logbook for the period 1987— 
2013 and by treating year as a random effect (Lynch et al., 
2018). The indices derived by Lynch et al. (2018) are cur- 
rently the most robust indices of relative abundance for 
dolphinfish in the western Atlantic Ocean. Although the 
majority of dolphinfish landings that are reported within 
the western Atlantic Ocean come from recreational fish- 
eries (Merten et al., 2022a), the U.S. PLL logbook remains 
the most robust data set because of the long time series of 
observations of dolphinfish made at the fishing set level, 
which is a fine temporal resolution, throughout the west- 
ern Atlantic Ocean. 

Spatiotemporal autocorrelation is often present in com- 
mercial data because sampling (fishing) by fishing fleets 
is preferential and not random (Maunder and Punt, 2004). 
Spatiotemporal autocorrelation can lead to positively 
biased estimates of abundance (Xu et al., 2018), partic- 
ularly for undersampled regions (Pennino et al., 2019). 
Lynch et al. (2018) noted that future work should consider 
other treatments of spatiotemporal data when estimating 
indices of relative abundance from fishery-dependent data 
and made specific reference to geostatistical delta-GLMMs 
(Thorson et al., 2015). Geostatistical delta-GLMMs 
include the use of vector autoregressive spatiotemporal 
(VAST) models that are built on a 2-stage modeling philos- 
ophy similar to that of a delta-GLMM. In a VAST model, 
covariates can be assigned specifically to effects related 
to environmental conditions or habitat and to catchabil- 
ity, spatial and spatiotemporal variation can be estimated 
explicitly as random effects, and spatiotemporal autocor- 
relation can be estimated (Thorson et al., 2015; Thorson, 
2019) instead of implicitly assuming that sampling over 
space is independent (Conn et al., 2017). In the model in 
the R package VAST, the space over which spatial and spa- 
tiotemporal processes are measured can also be treated as 
being evenly distributed to address the unbalanced and 
preferential nature of sampling in fishery-dependent data 
(Thorson, 2019). 

Findings from several studies indicate that spa- 
tiotemporal models can be used to reliably estimate 


indices of relative abundance from fishery-dependent 
data (Maunder et al., 2020; Ducharme-Barth et al., 2022; 
Hansell et al., 2022). The use of VAST models fit to data 
has generally resulted in improved estimation of abun- 
dance (Thorson et al., 2020; O’Leary et al., 2022) and 
of spatial population dynamics (Cao et al., 2017), with 
success demonstrated when models have been fit to 
fishery-dependent longline data for other highly migra- 
tory pelagic fish species, such as tunas (Maunder et al., 
2020; Hansell et al., 2022). Our goals were to quantify 
the spatiotemporal dynamics of dolphinfish in the west- 
ern Atlantic Ocean and to develop a model-based index of 
relative abundance that is unaffected by variation over 
space and time. Additionally, we sought to answer the fol- 
lowing questions with our study: 7 


1. Can the VAST model be used to estimate the sea- 
sonal patterns in the spatiotemporal dynamics of the 
density of the dolphinfish population in the western 
Atlantic Ocean? 

2. How do the spatiotemporal dynamics of dolphinfish 
differ at regional scales? 

3. Has a range shift occurred at the population level for 
dolphinfish in the western Atlantic Ocean? 

4. Do estimates from use of the VAST package reveal 
spatiotemporal signals similar to those found in 
other data sources? 


Materials and methods 
Fishery data 


For our analysis, we used U.S. PLL logbook data from the 
National Marine Fisheries Service (NMFS). Although the 
fishing effort of the U.S. PLL fleet is not uniform over 
space and time (Lynch et al., 2018), the logbook data 
contain 37 years of detailed catch and effort information 
on dolphinfish to the longline set level and compose the 
only data set with long-term observations collected over 
nearly the full extent of the geographic range of west- 
ern Atlantic dolphinfish. Additionally, the dolphinfish, 
although targeted by the fleet, is not the primary target 
species of the U.S. PLL fleet, which primarily targets the 
swordfish and tunas (Luckhurst, 2017). Therefore, dol- 
phinfish sampling (i.e., CPUE) is less likely to be biased 
by preferential sampling that results from the targeting 
behavior of the fleet. Vessels in the U.S. PLL fleet report 
the date and location of catch to the set level, SST, the 
type of gear and bait used, the number of hooks per set, 
whether dolphinfish were targeted or not, the number 
of dolphinfish caught, and the numbers of dolphinfish 
discarded dead and discarded alive. We used all logbook 
data from reports of catch in the western Atlantic Ocean 
for which latitude, longitude, SST, and CPUE informa- 
tion were available. The data we used contain approx- 
imately 255,000 observations, of which approximately 
160,000 were observations in which zero encounters with 
dolphinfish were recorded. 
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Nominal catch calculation 


We treated the number of hooks used in a longline set 
as the unit of effort (Lynch et al., 2012). We calculated 
the nominal CPUE as the number of dolphinfish caught 
per 1000 hooks, including any dolphinfish caught and 
retained or caught and discarded, whether dead or alive. 
We explored the number of set hours as an additional 
measure of effort; however, the results of analysis of catch 
versus hourly effort indicate a gear saturation effect in 
which CPUE declines exponentially after 5 h of soak time. 
This decrease in CPUE violates the assumption that the 
probability of catch is the same over every hour (Peterson 
et al., 2017; Hansell et al., 2022). 


Spatial domain 


The spatial domain included what we have defined as the 
western Atlantic Ocean, which extends approximately 
from 0°N to 55°N and from 25°W to 86°W. We stratified 
the spatial domain by using the discrete management 
areas for the U.S. PLL fishery described in the U.S. PLL 
logbook data set. These areas are defined, east to west and 
north to south, as follows: the region closest to the coast 
from Cape Cod, Massachusetts, to Cape Hatteras, North 
Carolina, which in the logbook data set is referred to as 
the Mid-Atlantic Bight (MAB); northeast coastal waters 
(NEC), farther to the north and east from the MAB; north- 
east distant waters on the Grand Banks (NED); the region 
closest to the coast from Cape Hatteras to the southern 
border of Georgia, described in the logbook data set as the 
South Atlantic Bight (SAB); waters off the Florida east 
coast and around the Florida Keys (FEC); the eastern 
part of the Sargasso Sea (SAR); the region in the north- 
central Atlantic Ocean that forms the western part of 
the Sargasso Sea (NCA); the Caribbean Sea (CAR); and 
offshore waters (OFS) (Fig. 1) (Damiano, 2023). To map 
data, we used a custom extrapolation grid to encompass 
all strata (Damiano, 2023). 


Seasonality 


The VAST model was configured to estimate indices for 
4 seasons to capture the pattern of seasonal availability 
throughout the range of dolphinfish in the Western Atlantic 
Ocean. We assigned seasonal [Ds to U.S. PLL logbook data 
at the monthly level: winter included December of a given 
year and January and February of the following year; spring 
included March, April, and May of the same year; sum- 
mer included June, July, and August of the same year; and 
fall included September, October, and November of the same 
year (Damiano, 2023). The VAST model was configured with 
a total of 148 time steps over 4 seasons and 37 years. We also 
explored a monthly time step. 


Spatiotemporal model description 


The models implemented in the VAST package (vers. 3.5.0; 
Thorson and Barnett, 2017) in statistical software R (vers. 


4.0.3; R Core Team, 2020) were fit to nominal CPUE values 
to derive standardized indices of relative abundance. The 
model in the VAST package functions like a delta-GLMM, 
also called a hurdle or delta model, that separately models 
the probability of non-zero catches (i.e., positive encoun- 
ters) and positive catch rates from each encounter (i.e., 
positive catch rates) (Thorson et al., 2015; Lynch et al., 
2018). We used the binomial GLMM with a logit link for 
the first linear predictor and explored lognormal, gamma, 
and Tweedie distributions for the second linear predictor. 
We chose the delta-lognormal configuration so that we 
could compare the results with those of Lynch et al. (2018) 
and could report the optimal model configuration, which 
was identified on the basis of model diagnostics, by using 
the following equations: 


5 ae [be (t;)+ @, (s,) +e, (s,, 4] 


" YG (p) X(p,5,,4)}, and (1) 


Poi = Be (t;) + Wo (s,)+ Eo (s,, 4] 
+"? Y2(p)X(p,5:,4), (2) 


where P, ; = the probability of positive dolphinfish encoun- 
ter for observation 7; 

P, ; = the positive dolphinfish catch rates (i.e., num- 
ber of dolphinfish caught per 1000 hooks) 
given that they were encountered, for obser- 
vation 1; 

B(¢,) = the intercept for each time step ¢ for obser- 
vation 1; 

«(s;) = the spatial variation that is constant over time 
for knot s for observation 1; 

e(s., t,) = the time-varying spatiotemporal autocorrela- 
tion variation for knot s and time step ¢ for 
observation 1; 

yp) = the habitat effect p on the univariate response 
variable; 

n, = the number of habitat covariates; and 
X(p, s;, t;) = the value of covariate p in knot s in time step 
t for observation 7. 


Note that B(¢,) was treated as an independent and iden- 
tically distributed random effect. Because of the seasonal- 
ity in availability of dolphinfish to the U.S. PLL fishery 
over space, we included season as a spatially varying ran- 
dom effect of habitat; this treatment is also the standard 
for seasonal implementations of a VAST model. The spa- 
tial processes, @,(s;) and @o(s;), and spatiotemporal pro- 
cesses, €,(s;, t,) and €,(s;, t,), were modeled as Gaussian 
Markov random fields (random effects) (Lindgren et al., 
2011) with an isotropic Matérn covariance function (.e., 
correlation over a 2-dimensional spatial grid). We con- 
ducted correlation analyses by using the COR (vers. 0.0.1; 
Guo, 2020) and corrplot (vers. 0.92; Wei and Simko, 2021) 
packages in R with spatially stratified nominal CPUE for 
each season to determine the appropriate spatiotemporal 
structure for the second linear predictor. Although 


30 


40°N- 


United States 2m 


Fishery Bulletin 122(1-2) 


Freq. 
ad of 
trips 
High 
| & Medium 


Low 


Atlantic Ocean 


<3 


50°W A0°W 30°W 


Figure 1 


Map of the relative frequency of commercial trips in which dolphinfish (Coryphaena hippurus) were caught 
during 1986—2022, based on data from the U.S. pelagic longline logbook for the following management areas: 
the region closest to the coast from Cape Cod, Massachusetts, to Cape Hatteras, North Carolina, which in 
the logbook data set is referred to as the Mid Atlantic Bight (MAB); northeast coastal waters (NEC) farther 
to the north and east from the MAB; northeast distant waters on the Grand Banks (NED); the region closest 
to the coast from Cape Hatteras to the southern border of Georgia, referred to in the logbook data set as 
the South Atlantic Bight (SAB); waters off the Florida east coast and around the Florida Keys (FEC); the 
eastern part of the Sargasso Sea (SAR); the region in the north-central Atlantic Ocean that forms the west- 
ern part of the Sargasso Sea (NCA); the Caribbean Sea (CAR); and offshore waters (OFS). Colors indicate 
relative frequency of trips made to the geographic coordinates within each hexagon. 


correlation was high among some strata during specific 
seasons, there was no consistent pattern (Suppl. Fig. 1) 
(online only); therefore, we assumed that spatiotemporal 
random effects were independent and identically distrib- 
uted. We fit the VAST model by using a mesh approach 
with 500 knots, the unit of measure over which to measure 
the spatial and spatiotemporal processes. 

Parameters were determined through maximum like- 
lihood estimation, in which fixed and random effects 
were integrated, by using the R package Template Model 
Builder (vers. 1.9.6; Kristensen et al., 2016). We used esti- 
mated log-scale densities to calculate indices as the pre- 
dicted local density d for knots s and time step ¢ multiplied 


by area a, summed over the total number of knots (n,). In 
the VAST model, standardized indices of relative abun- 
dance are labeled as biomass in kilograms by default 
when the response is continuous CPUE (Thorson, 2019). 
However, we rescaled the index (J) estimates by using a 
constant so that the values would be unitless: 


1(¢)= (als) <a(s)) ® 


We rescaled the index for each time step ¢ as a product 
of local density and area a over knots s summed over the 
total number of knots. We calculated the change in spatial 
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distribution using the center of gravity, which is the esti- 
mated centroid of the population’s distribution: 


. (2 (s) x a(s) oS d(s,t)} 


Z(t) = be = Niee sane aaese (4) 


where Z(t) = the center of gravity over time ¢ that sums 
the center of gravity statistic z multiplied by 
area and density over the index J over the 
number of knots s. 


The changes in population distribution were measured in 
kilometers in eastings and northings (Thorson et al., 2015). 
Following Li et al. (2014), we also calculated the center of 
gravity for the nominal effort data (Z,), using the effort con- 
tained in the U.S. PLL logbook data set, as the sum of the 
product of seasonal effort (#,) and space (i.e., latitude and 
longitude) (Z,) over the sum of total effort to explore spatial 
changes in the distribution of U.S. PLL effort: 


7 Ely 
TE, (5) 


Model selection, sensitivities, and diagnostics 


We identified the optimal configuration of the model 
described in the previous section on the basis of the per- 
centage of deviance explained and Akaike information 
criterion (AIC) (Akaike, 1973). In the “Results” section, we 
report these values in relation to the null model, which 
was constant intercept-only and required for calculating 
the percentage of deviance explained. We also provide 
results for a time-varying intercept-only model. Model fit 
was evaluated by using the R package for residual diagnos- 
tics, DHARMa (vers. 0.4.6; Hartig, 2022); diagnostics for all 
models were acceptable, and the optimal model diagnostics 
are provided in Supplementary Figure 2 (online only). 


Correlation analysis 


The final research question was whether the VAST estimates 
could pick up spatiotemporal signals in other sources of data 
on dolphinfish. In order to address this question, we used the 
COR and corrplot packages to analyze correlations between 
indices of estimated abundance by spatial strata and season 
with independent catch time series. These data include rec- 
reational catch and effort information for the period 2002— 
2019 from the Large Pelagics Survey of private and for-hire 
recreational vessels conducted by the Marine Recreational 
Information Program (NMFS, website), commercial landings 
for fishing operations that were not longline in spring and 
summer in North Carolina and Florida during 1991-2021 
(M. Karnauskas, personal commun.), and recreational dol- 
phinfish tournament data for the period 2000—2019 from 
the Marine Recreational Fishing Program of the Puerto 
Rico Department of Natural and Environmental Resources 
(Rodriguez-Ferrer®). Analysis of correlations with other fish- 
eries data sets, especially those from recreational sources, 
will help to determine if an index fit to commercial data is 
consistent with noncommercial fishery trends (Damiano, 
2023). A description of the indices and catch time series that 
were analyzed is provided in Table 1. 


Results 
Optimal model selection 
We experimented with several model configurations to 


explore the effects of covariates and alternate distribu- 
tions. Excluding intercept-only models, only 3 models 


3 Rodriguez-Ferrer, Y. 2022. Personal commun. Mar. Recreat. 


Fish. Progr., Puerto Rico Dep. Nat. Environ. Resour., 1375 Ave. 
Ponce de Leén, San Juan, Puerto Rico 00926. 


Table 1 


The time series and seasons used to estimate the indices of abundance, of dolphinfish (Cory- 
phaena hippurus) in the western Atlantic Ocean by spatial strata, for which correlation 
analyses were conducted. The spatial strata include the following regions: the region closest 
to the coast from Cape Cod, Massachusetts, to Cape Hatteras, North Carolina (NC), which 
in the U.S. pelagic longline (U.S. PLL) logbook data set is referred to as the Mid-Atlantic 
Bight (MAB); the region closest to the coast from Cape Hatteras to the southern border of 
Georgia, referred to in the U.S. PLL logbook data as the South Atlantic Bight; waters off the 
Florida east coast and around the Florida Keys (FEC); and the Caribbean Sea (CAR). The 


time series include data from the Large Pelagics Survey (LPS) of the Marine Recreational 
Information Program, commercial landings in NC and Florida (FL) from the NOAA South- 
east Fisheries Science Center, and data from the Marine Recreational Fishing Program 
(MRFP) of the Puerto Rico Department of Natural and Environmental Resources. 


Region Catch time series 


MAB LPS recreational catch and effort 
Hook-and-line catch in FL and NC 
CAR MRF? recreational tournament 


SAB, FEC 


Period Seasons 

Summer and fall 
Spring and summer 
Winter, spring, and fall 


2002-2019 
1991-2021 
2000-2019 
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converged to a solution. We attempted to fit set-level 
observations of SST using a 3-parameter basis spline 
function because of the association of dolphinfish with 
a relatively narrow, dome-shaped band of warm SSTs 
that varies over space and likely leads to positive catch 
rates (Farrell et al., 2014; Molto et al., 2020; Schlenker 
et al., 2021). We also experimented with treating SST 
as a spatially varying random habitat effect. Given that 
the U.S. PLL targeted dolphinfish in only approximately 
8% of observations during the period examined, we also 
tested one catchability covariate to account for variation 
in fleet behavior: a 3-level factor describing targeting 
behavior; in the U.S. PLL logbook data, most observa- 
tions are assigned a Yes or No for targeting dolphinfish. 
For each observation with no targeting assignment, we 
added the letter U, for Unknown, in the blank cell. In this 
configuration, targeting behavior was treated as a spa- 
tially varying linear random effect. We considered bait 
as a potential catchability covariate, but bait type was 
nearly uniform throughout the data set. Therefore, it was 
not included. 

We also attempted to model an individual vessel effect 
to account for spatiotemporal changes in fishing power 
(O’Leary et al., 2022). However, it either became computa- 
tionally infeasible when used as a random effect, resulting 
in hundreds of thousands of new parameters, or resulted 
in computer memory allocation issues when treated as a 
fixed effect, likely because of approximately 700 individ- 
ual vessels being treated as individual levels of a factor. 
We conducted a sensitivity analysis in which a geomet- 
ric anisotropy spatial smoother was used in the optimal 
model; however, a solution could not be found with this 
model. We also experimented with gamma and Tweedie 
distributions for the observation model in the second lin- 
ear predictor, but a solution could not be obtained with 
either model. 

The optimal model (i.e., the basic seasonal spatiotem- 
poral model) explained 30.8% of deviance and had an AIC 
of 880,221. The time-varying intercept model, which was 
used to explore the effect of temporal variation, explained 
9.9% of deviance and had an AIC of 962,118. Adding SST 
as a spatially varying linear random effect resulted in a 
model with 31.6% of deviance explained and an AIC of 
876,181. Although the percentage of deviance explained 
is greater and the AIC is lower for this model, the rela- 
tionship between catch rates for dolphinfish and SST is 
not linear (Farrell et al., 2014). Furthermore, the relation- 
ship between catch rates and SST is likely confounded 
with seasonality, and the additional 0.8% of deviance 
explained may account for the variation in set-level SST 
observations. 

The model that included both SST and targeting 
behavior as spatially varying random linear covari- 
ates related to habitat and catchability, respectively, 
explained 32.4% of deviance and had an AIC of 871,608. 
Although this model technically achieved the best per- 
formance metrics, using it did not result in different 
patterns in estimates of abundance or center of gravity. 
Still, use of this model did result in substantially larger 


standards for those estimated quantities, indicating 
that the addition of spatially varying random linear 
covariates, which added tens of thousands of: additional 
random effects for estimation, did not alter the pattern 
and overall magnitude of trends obtained from the basic 
seasonal spatiotemporal model. This result further indi- 
cates that variation in estimates due to environmental 
effects or targeting behavior of fleets was confounded or 
“absorbed” by seasonal, spatial, and spatiotemporal pro- 
cess error. Therefore, we chose the basic spatiotemporal 
VAST model, which was the most parsimonious model, 
as the optimal configuration. 


Spatiotemporal patterns in abundance 


The total estimated seasonal index of relative abundance 
for the western Atlantic Ocean (Fig. 2) was defined by 
several peaks and troughs with no trend until 2018, after 
which the magnitude of subsequent peaks and troughs 
declined. Overall, the greatest levels of total abundance 
of dolphinfish occurred during 1986, 1995, and 2015, and 
abundance was generally highest in spring. The lowest 
abundance occurred in 2020. The trend of the total index 
among seasons was similar but stratified by magnitude; 
total relative abundance was greatest in spring, followed 
by summer and fall in some years, and then winter. Strat- 
ified by region, indices had considerable variability in 
trend and seasonality. 

Estimates from northern coastal strata, such as the 
NEC, MAB, and SAB, had the highest variability among 
seasons, with higher abundances during spring and sum- 
mer and lower abundances during fall and winter (Fig. 2); 
trends in abundance in spring and summer were similar 
through 2012 but began to deviate afterward. This pattern 
was most pronounced in the SAB, where abundance in 
summer declined to fall and winter levels after 2014 and 
abundance in spring remained relatively high and rela- 
tively stable during 2020-2022. Additionally, it is worth 
noting that the magnitude of abundance in fall increased 
relative to other seasons toward the north and the east 
(e.g., the magnitude of abundance in fall rose from the 
SAB to the MAB, from the MAB to the NEC, and from the 
NEC to the NED) (Fig. 2). 

Estimates from the more southern strata and off- 
shore strata (i.e., FEC, SAR, NCA, CAR, and OFS) indi- 
cate weaker or low variability in seasonality and mirror 
the index of estimated total abundance that has no trend 
until the decline during 2019-2022 (Fig. 2). Estimates 
were most certain (had the lowest standard errors) for 
abundance in the MAB and SAB, followed by those for 
abundance in the FEC and CAR, and were most uncer- 
tain (had the highest standard errors) for abundance in 
the offshore strata NED, SAR, and NCA (Fig. 2). These 
patterns in uncertainty were most likely due to higher 
and lower rates of sampling (i.e., fishing within those 
strata) respectively (Fig. 1). Furthermore, the highest 
estimates of abundance relative to the total index were 
found for the NCA, likely influencing the shape of the 
overall estimated trend in seasonal abundance. 
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Shift in center of gravity 


For the indices of estimated abundance, the eastings that 
approximate the spatial distribution from west to east 
from the bottom of the y-axis have no obvious trend 
(Fig. 3A). The northings that approximate the spatial dis- 
tribution from south to north from the bottom of the 
y-axis also have no trend (Fig. 3B). For the nominal 
CPUE, the center of gravity in eastings indicates a strong 
westward shift (~120 km) during 1986-1994, then a 
gradual eastward shift during 1995-2022 (~90 km) that 
was consistent among seasons (Fig. 3C). The center of 
gravity in northings for nominal CPUE indicates an 
overall southern shift during 1993-2022 (~40 km), with 
some seasonal stratification: trends were similar between 
summer and winter and between spring and fall 
(Fig. 3D). We used box plots to visualize the average sea- 
sonal pattern in centers of gravity for the index of esti- 
mated abundance and for nominal CPUE (Fig. 4, A—D). 
The average seasonal pattern in the center of gravity for 
the estimated abundance index was defined by abun- 
dance shifting north and west as winter transitions to 
spring, slightly south and east during the transition to 
summer, and then south and east during fall (Fig. 4, A 
and B). For the nominal CPUE, the average seasonal pat- 
tern in the box plots more closely mirrored the time 
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series shown in Figure 3, C and D, with a net southeast- 
ern movement throughout the year (Fig. 4, C and D). 


Correlation analysis 


Correlations among the VAST-estimated index and indices 
based on data sources independent of the U.S. PLL log- 
book, by spatial strata, were consistent with the timing of 
large migrations of dolphinfish and the overall seasonal 
pattern in the center of gravity (Fig. 4, A and B). The 
strong positive correlations between the Marine Recre- 
ational Fishing Program tournament catches in fall and 
the abundance estimates for winter and spring in the CAR 
(Fig. 5) are consistent with the timing of the movement of 
dolphinfish through the Antilles Current, which flows 
along the archipelago that separates the Caribbean Sea 
from the Atlantic Ocean, and with lower seasonal varia- 
tion in abundance in this region compared with the varia- 
tion in other regions (Fig. 3). Similarly, hook-and-line 
catches in North Carolina and Florida had the strongest 
correlations with estimates of abundance in the FEC and 
SAB during spring and summer months, and the CPUE 
data from the Large Pelagics Survey was most correlated 
with the estimate of abundance for the MAB in the fall 
(Fig. 5). Weak or negative correlations generally occurred 
between the estimated abundance index and the indices 
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Figure 3 


Seasonal estimates of the centroid or center of gravity of the distribution of dolphinfish (Coryphaena hippurus), by (A) eastings 
and (B) northings, in the western Atlantic Ocean during 1986-2022, based on abundance indices created with a vector autore- 
gressive spatiotemporal model. Also shown are the seasonal centers of gravity calculated for nominal catch per unit of effort 
from the U.S. pelagic longline logbook for the same period, by (C) eastings and (D) northings. In the top panels, a gray area for 
each season outlined with colored dashed lines indicates the 95% confidence interval for that season. 
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Figure 4 


Box plots of seasonal estimates of the centroid or center of gravity of the distribution of dolphin- 
fish (Coryphaena hippurus), by (A) eastings and (B) northings, in the western Atlantic Ocean 
during 1986—2022, based on abundance indices created with a vector autoregressive spatiotem- 
poral model. Also shown are the seasonal centers of gravity calculated for the nominal catch per 
unit of effort from the U.S. pelagic longline logbook for the same time period, by (C) eastings and 
(D) northings. The thick line within the box indicates the median, the upper and lower parts of the 
box represent the first and third quartiles (the 25th and 75th percentiles), the whiskers extending 
above and below the box correspond to 1.5 times the interquartile range, and points represent 


values outside this range. 


based on other data sources when season or region did not 
overlap (Fig. 5). Additionally, results of the correlation 
analysis provide some validation for estimated spatiotem- 
poral patterns in abundance and indicate that using a 
VAST model to estimate spatial and spatiotemporal ran- 
dom effects can reveal broader spatiotemporal signals in 
the population dynamics of dolphinfish. 


Discussion 


Our model results reveal several important features of the 
population dynamics for western Atlantic dolphinfish. The 
overall abundance of dolphinfish in the western Atlantic 
Ocean was relatively stable during 1986-2018, followed 
by a decline during 2019-2022; seasonal patterns in abun- 
dance vary at different spatial scales. The stable trend in 
estimated abundance for 1986-2018 was somewhat sur- 
prising given the more pronounced negative trends found 
by Kleisner (2008) and Lynch et al. (2018). However, as of 


2023, the U.S. PLL logbook has become a much larger data 
set than the one that was available during either of these 
previous studies, and in light of the results of our study 
(Fig. 2), the negative trends identified previously may have 
been driven by natural downward fluctuations in abun- 
dance that happened to occur during the terminal year of 
each analysis. Apart from the overall negative trend, our 
results are consistent with the general pattern in estimated 
abundance of dolphinfish that has been reported for the 
periods 1987—2005 (Kleisner, 2008) and 1987-2013 (Lynch 
et al., 2018). Furthermore, because we modeled spatial and 
spatiotemporal process errors, our results likely explain 
variation in numerous albeit aggregated latent processes 
(e.g., variation due to the behavior of the Gulf Stream or 
changes in fleet composition and behavior) that could not 
be explicitly modeled in previous research efforts. 
Movement of dolphinfish generates seasonally and spa- 
tially explicit availability of abundance to many fisheries 
(Oxenford and Hunte, 1986; Mahon, 1999). Our results 
(Fig. 4, A and B) matched expectations for both the 
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Figure 5 


Correlations between indices of estimated abundance of dolphinfish (Coryphaena hippurus) in the western 
Atlantic Ocean and independent catch time series, by spatial strata and season. The correlations on the left 
side are between the indices of estimated abundance in summer and fall in the region closest to the coast 
from Cape Cod, Massachusetts, to Cape Hatteras, North Carolina, which in the U.S. pelagic longline logbook 
data set is referred to as the Mid-Atlantic Bight (MAB), and catch and effort data for the period 2002-2019 
from the Marine Recreational Information Program’s Large Pelagics Survey (LPS). The correlations in the 
middle are between estimated abundance indices for spring and summer in the region closest to the coast 
from Cape Hatteras to the southern border of Georgia, referred to in the logbook data set as the South Atlan- 
tic Bight (SAB), and in waters on the east coast of Florida and around the Florida Keys (FEC) and commer- 
cial catches (comm.) in spring and summer in North Carolina (NC) and Florida (FL) during 1991-2021. The 
correlations on the right side are between indices of estimated abundance in the fall, winter, and spring in 
the Caribbean Sea (CAR) and recreational tournament catch data for the same seasons during 2000-2019 
from the Marine Recreational Fishing Program (MRFP) of the Puerto Rico Department of Natural and 


Environmental Resources. 


seasonal distribution of abundance of dolphinfish and the 
spatiotemporal pattern of movement outlined by Oxen- 
ford and Hunte (1986) for the northern stock and 
expanded on by Merten et al. (2014a, 2014b, 2016). The 
estimates of abundance during spring and summer 
within strata along the East Coast of the United States 
were consistent with documented patterns of seasonal 
abundance (Merten et al., 2014a, 2015, 2016). Addition- 
ally, movement rates of dolphinfish slow in the SAB 
because of the Charleston Gyre (Merten et al., 2014a), an 
eddy off the coast of Charleston, South Carolina, that 
may be responsible for increased availability for extended 
pulses of fishing effort in spring and summer (Damiano, 
2023). In the northern regions, the magnitude of abun- 
dance in spring and summer is substantially greater 
than in fall and winter, but this seasonal difference in 
magnitude (i.e., inter-seasonal variability) is lower in the 
southern regions. 

The timing of the movement of the northern stock through 
the Caribbean Sea generates broad seasonal availability 
of dolphinfish to fisheries (Oxenford and Hunte, 1986), as 
reflected in estimates of abundance for southern spatial 
strata (Damiano, 2023). For example, although estimates 
of abundance in spring in the CAR were still the greatest 
among seasons, the magnitude of abundance in fall and 
winter were comparable with that of abundance in spring. 
These peaks in abundance in fall and winter may coincide 


with the movement of dolphinfish traveling to and from 
the Strait of Florida (Merten et al., 2014b, 2015) and the 
appearance of large dolphinfish (>90 cm in total length) in 
the waters of Caribbean islands (such large dolphinfish are 
caught in recreational tournaments; W. Merten, personal 
commun.), but these peaks may be more indicative of pat- 
terns in abundance of dolphinfish in certain U.S. territo- 
ries in the Caribbean region than patterns in the broader 
CAR stratum (Damiano, 2023). 

The increases in abundance in spring in the MAB and 
SAB during 2010-2018 and the associated decrease in 
abundance in summer in the SAB are also consistent with 
perceptions of commercial, recreational, and for-hire fish- 
ermen regarding local availability of dolphinfish in waters 
of North Carolina, including those north of the Outer 
Banks; they have observed either an increase in or stable 
catches of dolphinfish during the period approximately 
from 2015 through 2018 or 2019 (McPherson et al., 2022). 
The linear negative trend in abundance in the FEC region 
during 2015-2022 is similarly consistent with observa- 
tions by stakeholders from eastern Florida and the Flor- 
ida Keys of lower catches of dolphinfish in recent years 
leading up to a workshop held in 2021 (McPherson et al., 
2022). Some local perceptions of abundance in Puerto 
Rico also have matched the recent decrease in estimated 
abundance in the CAR (Merten et al., 2022a; W. Merten, 
personal commun.). 
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These stakeholder observations provide some validation 
for the regional estimates of spatiotemporal indices of rel- 
ative abundance (Damiano, 2023). Additionally, it is worth 
noting that the estimated trends in total abundance, and 
in abundance in the MAB, SAB, and FEC, strongly resem- 
ble those of the catch time series for dolphinfish from 
the Marine Recreational Information Program (NMFS, 
website). Further, the results of our correlation analysis 
indicate that the modeled spatial and spatiotemporal pro- 
cesses are identifying signals observed in data for other 
fisheries and in the Large Pelagics Survey that are inde- 
pendent of the U.S. PLL fleet. However, additional infor- 
mation, specifically from other Caribbean island nations 
and countries bordering the Caribbean Sea, would be 
required to determine if the index for the CAR is suffi- 
ciently representative of spatiotemporal dynamics of dol- 
phinfish in the western extent of the Caribbean Sea where 
the U.S. PLL fleet fishes less than it does in U.S. territories 
(Fig. 1) (Damiano, 2023). 

We suggest that there are 2 potential explanations (which 
are not mutually exclusive) for the recent decrease in abun- 
dance of dolphinfish during 2019-2022. The recent decline 
may be a lagged effect of a protracted period of high fish- 
ing pressure by commercial, recreational, and artisanal 
sectors throughout the western Atlantic Ocean (Damiano, 
2023). Private recreational fishing in the coastal United 
States has been the dominant source of reported remov- 
als of dolphinfish in the western Atlantic Ocean (Merten 
et al., 2022a), and recreational effort linearly increased 
during 1986-2010 and then plateaued (NMFS, website). 
In the northeastern states located within the NEC and 
MAB, recreational catches have remained high (NMFS, 
website) despite the estimated decreases in abundance 
(Fig. 2). Consequently, these increases in recreational 
exploitation rates in specific areas may explain the recent 
reductions in abundance (Fig. 2, Suppl. Fig. 3 [online only]). 

With respect to commercial effort, although the num- 
ber of trips of the U.S. PLL fleet has decreased since 1994 
and effort has linearly declined (senior author, personal 
observ.), commercial fishing throughout the northwest 
Atlantic and western Central Atlantic Oceans in numer- 
ous countries in aggregate has linearly increased (FAO 
FishStat, website); this trend does not include those coun- 
tries where fishing occurs but landings of dolphinfish are 
not reported (Merten et al., 2022a). In Barbados, landings 
of dolphinfish have slowly decreased while CPUE has 
remained relatively stable during 1994-2018 (Oxenford 
et al.*), with higher sustained effort coupled with reduced 
landings (i.e., availability of fish), indicating the possibil- 
ity of an effect on stock productivity. 

Concerning the effect of artisanal fisheries, artisanal 
fishing in Venezuela also is reportedly a large source of 
removals of dolphinfish in the Caribbean Sea (Arocha 


* Oxenford, H. A., D. Johnson, S.-A. Cox, and J. Franks. 2019. 
Report on the relationships between Sargassum events, oceanic 
variables, and dolphinfish and flyingfish fisheries, 25 p. Cent. 
Resour. Manage. Environ. Stud., Univ. West Indies, Cave Hill 
Campus. Bridgetown, Barbados. [Available from website.] 


et al., 2021), and other Caribbean island countries are 
reportedly developing artisanal fleets that target pelagic 
fish, including dolphinfish, for which data are not avail- 
able (W. Merten, personal commun.). Cumulatively, a sce- 
nario in which overexploitation is occurring throughout 
the western Atlantic Ocean is plausible, although this 
notion is subject to a high degree of uncertainty because of 
the lack of reporting. 

Alternatively, or in concert with fishing activity, envi- 
ronmental effects may explain the overall trend and 
recent decline in estimated abundance. Mean SST is rising 
throughout the world (Garcia-Soto, 2021) and in the large 
marine ecosystems that compose the western Atlantic 
Ocean (Kessler et al., 2022). Within the western Atlantic 
Ocean, the most pronounced rates of increased warming 
are occurring in the North and northwest Atlantic Oceans 
(Saba et al., 2016; Kessler et al., 2022) due to a combina- 
tion of natural climate variability and anthropogenic cli- 
mate change (Chen et al., 2020). Consequently, northward 
shifts in the distribution of abundance have been observed 
for many species along the continental shelf of the Atlantic 
coast of the United States (Kleisner et al., 2017; Morley 
et al., 2018). Additionally, the heat transported northward 
by the Gulf Stream that had normally been quite variable 
has steadily declined since 2008 (Bryden et al., 2020), with 
a pronounced decrease occurring in 2020 that has since 
been followed by a period of uncharacteristically low heat 
transport (NOAA Atlantic Oceanographic and Meteorolog- 
ical Laboratory, website). 

Dolphinfish appear to prefer warm water within a rela- 
tively narrow band (Molt6 et al., 2020), and like individuals 
of other marine fishes, they will seek to remain within their 
optimal thermal habitat (Schlenker et al., 2021). However, 
temperature preferences of populations of dolphinfish are 
not geographically uniform; for example, dolphinfish prefer 
waters with temperatures of 22—24°C off the coast of Peru 
(Torrején-Magallanes et al., 2019), but dolphinfish have 
been caught in waters as warm as 36—37°C in the Gulf of 
Oman, an arm of the Arabian Sea (W. Merten, personal com- 
mun.). A recent lab experiment conducted by Heuer et al. 
(2021) off eastern Florida confirmed that a dolphinfish’s 
ability to physically move was optimal within the well- 
documented preference of western Atlantic dolphinfish for 
waters with temperatures of 24—29°C (Farrell et al., 2014; 
Merten et al., 2014b; Schlenker et al., 2021). However, expo- 
sure to temperatures beyond 31°C has resulted in reduced 
mobility or death for this species (Heuer et al., 2021). 
These results may explain the reduction in abundance of 
dolphinfish in spring in the FEC during 2015-2022 and 
subsequent divergence in abundance in spring and sum- 
mer in the SAB (Fig. 2); the timing of fish migrating north 
from the Antilles may be occurring earlier during the year 
so that dolphinfish can avoid inhospitable temperatures. 
Salvadeo et al. (2020) estimated a northward shift in the 
population of dolphinfish in the eastern Pacific Ocean at 
similar latitudes. These distribution changes across differ- 
ent ocean basins are consistent with the broadscale effects 
of climate change on marine species globally (Poloczanska 
et al., 2013; Damiano, 2023). 
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A shift in population distribution is also consistent 
with stakeholder observations in the FEC region: fisher- 
men from the Florida Keys could not find dolphinfish in 
waters warmer that than the temperatures preferred by 
dolphinfish (>29°C), and fishermen from more northern 
regions (e.g., NEC, MAB, and SAB) have perceived more 
stable abundance and availability and done so earlier in 
the year (McPherson et al., 2022). Although it is highly 
plausible that dolphinfish exposed to near-lethal tempera- 
tures (>31°C) would move, the overall pattern in estimated 
abundance closely mirrors that of an SST “fingerprint” 
that indicates the recent weakening of the Atlantic Meridi- 
onal Overturning Circulation and consists of cooling in the 
subpolar Atlantic Ocean and warming in the Gulf Stream 
region (Caesar et al., 2018; Ditlevsen and Ditlevsen, 2023). 
The behavior of these oceanographic features may indicate 
a more complex relationship between the Gulf Stream and 
spatiotemporal dynamics of dolphinfish. 

Changes in fishing pressure and in the environment are 
likely contributing to the observed patterns in abundance, 
and data gaps make it challenging to discern the drivers of 
shifts with more accuracy. Additionally, some of the changes 
in fishing pressure and in the environment are aligned 
spatially and temporally, making it more challenging to 
disentangle the forces. For example, increases in landings 
from the Eastern Caribbean in 2014 coincided with signif- 
icant changes in temperature, upwelling, and productivity 
in the region from Cape Hatteras through the Florida Keys 
(Craig et al., 2021). Sudden warming of the Florida Straits 
(Volkov et al., 2019), which connects the Gulf of Mexico and 
Atlantic Ocean and is a major migration pathway for dol- 
phinfish, may also have contributed to changes in presence 
of dolphinfish. Most of the variation in these overlapping 
processes is likely absorbed by the estimation of spatial 
and spatiotemporal process errors with small contribu- 
tions to deviance (i.e., variation) explained by modeling 
explicit habitat and catchability covariates. Differentiating 
between the drivers related to habitat and catchability will 
require more intensive analysis of patterns in ocean condi- 
tions and fishing exploitation in relation to patterns in the 
abundance discussed herein. 

We recommend further investigation into the relation- 
ship between SST anomalies and the behavior of the 
Atlantic Meridional Overturning Circulation and Gulf 
Stream and the abundance of western Atlantic dolphin- 
fish by using fine-scale fishery-independent oceanographic 
and movement data. Such an analysis would allow further 
differentiation between changes in fishing pressure and 
changes in the environment as drivers of distribution pat- 
terns and would provide some insight into whether these 
patterns will continue in the future given the shifts in tem- 
perature and current systems that are predicted to occur 
with climate change (Alexander et al., 2020; Damiano, 
2023). Spatiotemporal models may not be appropriate 
for these analyses because of the tendency to aggregate 
estimated spatial and spatiotemporal process errors. Fur- 
thermore, although the seasonal patterns in the estimated 
center of gravity (Fig. 4, A and B) indicate that a general- 
ized movement circuit was captured with the VAST model, 


we acknowledge that movement of dolphinfish is likely 
much more complex than the patterns described herein, 
potentially with patterns markedly different from those 
reported in the literature (Peterson”) and may include con- 
nectivity with regions not included in this analysis (e.g., 
the Gulf of Mexico and West Africa). 

The configuration of the VAST model used in this study 
and described herein does not include any age or size 
structure in the population. Therefore, any spatiotempo- 
ral differences in ontogeny (e.g., documented differences 
in feeding habits) (Pincay-Espinoza and Varela, 2022), in 
spawning and larval distribution (Park et al., 2017), in 
changes in growth rate with size and age, or in vertical 
habitat use (1.e., association with fish aggregating devices 
or Sargassum) (Ortega-Garcia et al., 2022) were also 
likely to have been absorbed by the estimation of spatial 
and spatiotemporal random effects. No relationship was 
detected between influxes of Sargassum and CPUE for 
dolphinfish in the Eastern Caribbean (Oxenford et al.°), 
but large quantities of Sargassum may make fishing 
more difficult (.e., fishing methods are negatively affected 
and fish are harder to detect and, therefore, less likely to 
be encountered). This notion warrants some investiga- 
tion into the precise nature of the effects of Sargassum 
on catchability. We recommend that future research pri- 
oritize collection of age- or length-structure information 
sufficient for estimating age or size bins by using VAST 
or other spatiotemporal modeling techniques (Cao et al., 
2020; Hansell et al., 2022). 

The U.S. PLL logbook data have substantial spatial lim- 
itations and, therefore, uncertainty in eastern areas out- 
side of the U.S. exclusive economic zone and the western 
extent of the Caribbean Sea (Fig. 1). This uncertainty is 
evidenced by stronger extrapolation to the CAR and OFS 
regions and by larger standard errors and confidence inter- 
vals for regionally stratified indices from those regions 
(Fig. 3). Furthermore, although the average pattern of 
the seasonal center of gravity (Fig. 4, A and B) appears to 
have captured the movement circuit of the northern stock, 
the dynamics of a purported southern stock found in the 
CAR and OFS strata are likely confounded (Oxenford and 
Hunte, 1986). Despite this caveat, fitting a VAST model to 
these data may help inform the estimation of abundance 
in neighboring cells in the model’s extrapolation grid if 
both stocks occur. Future efforts to quantify the spatio- 
temporal dynamics of western Atlantic dolphinfish may 
benefit from calibrating the CPUE data from the U.S. PLL 
logbook with foreign longline fleet data sets from areas 
with low U.S. PLL fleet activity (O’Leary et al., 2022). Such 
an exercise will also help to determine the accuracy of this 
study’s model extrapolations to undersampled regions 
and consistency with signals appearing in other data sets. 
Lastly, we acknowledge that the optimal VAST model 
was sensitive to different configurations, smoothers, and 
assumptions concerning the second linear predictor’s error 


° Peterson, C. 2023. Personal commun. Southeast Fish. Sci. Cent., 
Natl. Mar. Fish. Serv., NOAA, 101 Pivers Island Rd., Beaufort, 
NC 28516-9722. 


Damiano et al.: Spatiotemporal dynamics of Coryphaena hippurus in the western Atlantic Ocean Sy 


structure. Future efforts dedicated to teasing apart fish- 
ery and environmental drivers of the population dynamics 
of western Atlantic dolphinfish by using models that are 
not spatiotemporal may help elucidate reasons for these 
sensitivities. 

The development of an index of relative abundance for 
western Atlantic dolphinfish through modeling with spa- 
tiotemporal variation is a vital step toward understanding 
the population dynamics of western Atlantic dolphinfish 
and improving the robustness of management measures 
(Damiano, 2023). Ideally, a spatially explicit stock assess- 
ment would be the next step toward determining stock 
status. However, despite advances in modeling software 
and in technology available for data collection and despite 
a growing citizen science data collection program (Merten 
et al., 2022b), the quality and temporal resolution of most 
data sources for western Atlantic dolphinfish have not 
changed radically since previous assessments were con- 
ducted in 2000 and 2006. Therefore, assessments would 
likely continue to be confined to surplus production mod- 
els. Kleisner (2008) found that fitting surplus production 
models to geostatistical indices did not improve the abil- 
ity to obtain reference points. Although McDonald et al. 
(2021) demonstrated that a spatially explicit biomass 
dynamics model can capture time-varying spatial pat- 
terns in both population processes and fishing effort, most 
removal data are available only at an annual time step, 
precluding modeling of the seasonal dynamics reported in 
this study. 

If improved data collection or monitoring are not possi- 
ble or stock status cannot be determined from an assess- 
ment, an alternative and perhaps more realistic use of 
the spatiotemporal index of relative abundance could be 
to guide an index-based management approach to setting 
catch advice (Apostolaki and Hillary, 2009; Little et al., 
2011; Fischer et al., 2020; SAanchez-Marono et al., 2021; 
Damiano, 2023). Indeed, empirical index-based man- 
agement procedures that have been simulation-tested 
through a management strategy evaluation have been 
shown to have promise for data-limited stocks (Geromont 
and Butterworth®; Sagarese et al., 2019). For western 
Atlantic dolphinfish, correlations at various time lags 
should be analyzed to identify appropriate lead times, 
and rigorous evaluation through simulations or a man- 
agement strategy evaluation should be pursued to test 
the robustness of such management procedures to alter- 
native spatial and temporal scopes (Bosley et al., 2019; 
Cadrin et al., 2020) and to the substantial uncertainty 
found in this study and described herein. Additionally, 
a management strategy evaluation would be useful in 
determining if an empirical management procedure will 
be more or less responsive than, or as responsive as, the 
current catch-averaging strategy for quantifying the spa- 


6 Geromont, H., and D. Butterworth. 2015. A review of assess- 
ment methods and the development of management procedures 
for data-poor fisheries, 198 p. University of Cape Town, Cape 
Town, South Africa. [Report produced for the FAO.] [Available 
from website.] 


tiotemporal dynamics of western Atlantic dolphinfish 
(Walter et al., 2023). 


Conclusions 


The dolphinfish possesses a unique life history character- 
ized by fast growth and early maturation (Oxenford, 1999; 
Perrichon et al., 2019) that hypothetically makes individ- 
uals of this species capable of withstanding high levels of 
fishing pressure (Schwenke and Buckel, 2008; Damiano, 
2023). The population abundance and distribution of dol- 
phinfish in the western Atlantic Ocean appear to have been 
relatively stable during 1986-2018, but that period was 
followed by one of declining abundance during 2019-2022. 
High fishing pressure in concert with changing ocean con- 
ditions may be responsible for the recent decline in and 
changes to abundance at finer spatial scales (Damiano, 
2023). Having used VAST models in our study, we provide 
the first spatiotemporal index of relative abundance for 
dolphinfish caught in the western Atlantic Ocean at sea- 
sonally and spatially explicit scales. The general consis- 
tency of our results with those of past studies and with 
stakeholder observations and the relatively strong spatio- 
temporal correlations with recreational catch time series 
indicate that our index is robust and reliably captures 
the spatiotemporal dynamics of western Atlantic dolphin- 
fish (Damiano, 2023). This research addressed the need 
for a reliable index of relative abundance and improved 
upon the work of Kleisner (2008) and Lynch et al. (2018) 
by accounting for spatial and spatiotemporal variation 
in abundance of dolphinfish. Spatiotemporal modeling 
techniques, although complex and data-intensive, provide 
more information than nonspatial configurations (Berger 
et al., 2017) and are especially useful for standardizing 
fishery-dependent data (Maunder et al., 2020; Izquierdo 
et al., 2022) when improved monitoring programs are too 
costly (Damiano, 2023). 


Resumen 


El dorado (Coryphaena hippurus) se captura en todo el 
océano Atlantico occidental a distintas escalas espacia- 
les y temporales. Los intentos previos de cuantificar la 
dindmica poblacional del dorado en la regién se han visto 
inhibidos por la incapacidad de modelar la dinamica espa- 
ciotemporal de esta poblacién. Ajustamos un modelo vec- 
torial espaciotemporal, estacional y autorregresivo (VAST) 
para cuantificar la dinamica espaciotemporal del dorado 
del Atlantico occidental, para estimar los indices abun- 
dancia relativos estandarizados de durante 1986—2022 a 
escalas regionales y estimar los cambios en la distribu- 
cidn espacial. La magnitud de la abundancia fue mayor 
durante la primavera y el verano en los estratos espacia- 
les del norte y fue comparable a lo largo de las estaciones 
en los estratos espaciales del sur. La abundancia del 
dorado parecié estabilizarse durante 1986—2018y luego 
disminuy6 durante 2019-2022. Esta tendencia ocurri6 
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en todas las regiones, excepto en las aguas del Atlantico 
desde el cabo Hatteras, Carolina del Norte, hasta la fron- 
tera sur de Georgia, donde la abundancia se mantuvo 
estable durante 2019-2022. No se detect6 ningtin cambio 
en la distribucién de la poblacion, pero los patrones regio- 
nales de abundancia proporcionan informacion sobre los 
cambios en el momento de la disponibilidad. En el pre- 
sente estudio se reporta el primer indice de abundancia 
relativa estandarizado que capta la dinamica espaciotem- 
poral del dorado del Atlantico occidental. Estos resultados 
incrementaron nuestra comprensién sobre la dinémica 
poblacional de esta especie en la regién y deberian resul- 
tar utiles en futuros intentos de manejar la poblacién a 
diferentes escalas espaciales y temporales. 
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Abstract—The white hake (Urophy- 
cis tenuis) is a groundfish distributed 
throughout the Gulf of Maine. Catch 
advice is based on stock assessments 
done with age-based population dynam- 
ics models; however, otolith aging is 
challenging because of unclear growth 
increments. To address this concern, we 
compared the consistency of aging with 
counts of visual annuli to that of aging 
with cycles of elemental concentra- 
tions measured by using laser ablation 
inductively coupled plasma mass spec- 
trometry. We tested the hypothesis that 
oscillations in both environmental con- 
ditions and internal physiology through 
time influence uptake of elements 
during otolith mineralization. Concen- 
trations of manganese, in comparison 
with those of the other investigated 
trace elements (magnesium, strontium, 
and barium), had the most promising 
correlation with visual growth incre- 
ments (~100% age agreement, +1 year), 
offering an additional tool to enhance 
increment identification. In our exam- 
ination of 550 otoliths collected during 
2007-2021, we found that white hake 
lived a maximum of 10.3 years and 
exhibited sexual dimorphism in maxi- 
mum length and age. By using gener- 
ated von Bertalanffy growth functions, 
L@140Geen ee) for males and 
L@)=140(1-e 18") for females 
(where L(t) is length at time 12), size 
and age at maturity were calculated 
for males (37.4 cm in total length (TL), 
3.3 years) and females (47.4 cm TL, 
4.2 years). These results demonstrate 
that otolith geochemistry can be used 
to improve the accuracy and precision 
of the estimation of fish age and matu- 
rity, even for challenging species. 
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The use of biological hard struc- 
tures, such as otoliths, is the most 
common method of aging bony fish 
(Pannella, 1971; Gauldie and Nelson, 
1990; Stevenson and Campana, 1992), 
dating back to the late 18th century 
(Hederstrém, 1959; Jones, 1992). Oto- 
liths are made of the biomineral calcium 
carbonate and form part of the mechano- 
receptor system in bony fish (Campana, 
1999). These biomineral structures are 
well established aging tools because 
they grow in proportion to the individ- 
ual with continuous accretion (Pannella, 
1971; Stevenson and Campana, 1992; 
Vigliola and Meekan, 2009). 

Like the fish itself, the rate at which 
the otolith grows is not always constant, 
changing because of a range of internal 
(age, sex, metabolism, and reproduc- 
tive status) and external (environmen- 
tal conditions and prey availability) 
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factors (Campana and Neilson, 1985; 
Neilson and Geen, 1985; Wright, 1991; 
Yamamoto et al., 1997; Wright et al., 
2001; Hiissy, 2008). Such variability in 
growth rate by an individual is tran- 
scribed into otolith formation through 
physiological processes, linking otolith 
growth to time (Pannella, 1971). For 
species that inhabit temperate climates, 
seasonal oscillations in environmental 
conditions translate to variations in 
growth, which are observed as distinct 
translucent (winter) and opaque (sum- 
mer) growth increments (visible under 
reflected light) that, as a pair, repre- 
sent a single year of growth (Victor 
and Brothers, 1982; Campana, 1999; 
Wright et al., 2002a). Although well 
studied, the drivers of seasonal growth 
increments are still debated and have 
been shown to be highly species spe- 
cific (Ding et al., 2020). Species-specific 
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studies that validate that one annulus spans one year of 
growth are required for accurate age estimation through 
otolith annuli counting. 

Although this process of counting growth increments 
in otoliths is well established, it remains challeng- 
ing in application for certain species like the white 
hake (Urophycis tenuis), whose age has been validated 
in waters of Canada (Hunt, 1982) but not in the Gulf 
of Maine (GOM). In past research, white hake were 
observed to have a combination of false and weak growth 
increments consistent among sexes, sizes, and popula- 
tions (Clay and Clay, 1991). Unlike true growth incre- 
ments, false growth increments do not form a full ring 
around the otolith core (Katayama, 2018). The biologi- 
cal mechanism for false growth increment deposition 
is still debated, with research results indicating that a 
combination of physiological (feeding and lifestyle tran- 
sitions) or environmental (temperature, salinity, and 
oxygen) conditions are typically responsible (Watson, 
1967; Berghahn, 2000; (Cappo et al., 2000; Wright et al., 
2002b; Katayama and Isshiki, 2007; Katayama, 2018). 
Weak growth increments are true growth increments 
that are indistinct between the opaque and translucent 
zones, likely because of decreased metabolism between 
seasons with age, and are ultimately difficult to identify. 
The occurrence of either false or weak growth increments 
raises the potential for both overestimating and under- 
estimating ages, resulting in bias in growth models. 

For white hake, the appearance of the first annulus 
has been characterized in detail by Lang et al. (1996), 
and identification of the first several annuli has been cor- 
roborated by the modes in length—frequency plots (Hunt, 
1982). Although these 2 studies serve as a foundation for 
the validation of the frequency of annulus formation in 
white hake, continued research efforts are needed to pro- 
vide evidence that assigned ages can be validated through 
aging with the entire otolith section. 

The white hake is an ecologically and economically 
important groundfish species found within the GOM 
(Ames, 2012; Ames and Lichter, 2013) and is managed by 
the New England Fishery Management Council. With a 
full range from Florida to Iceland, this species has been 
found along the eastern coast of North America and pri- 
marily inhabits muddy substrates along the continental 
shelf and slope (Musick, 1974; Ames, 2012). Limited age 
and growth research has been conducted for the neighbor- 
ing northern stock of white hake in the Gulf of St. Law- 
rence (Clay and Clay, 1991), but the GOM stock has not 
been fully explored because of the difficulty in aging this 
species (Penttila, and Dery, 1988), limiting age-based 
fishery management efforts. For a comparable species 
for which ages have been difficult to assign, the Atlantic 
cod (Gadus morhua), analysis of otolith geochemistry has 
been applied recently as a new tool (Heimbrand et al., 
2020; Hiissy et al., 2021). 

Otolith core-to-edge geochemistry has been used to 
determine life history information for a range of species. 
Otoliths are metabolically inert; therefore, incorporated 
trace elements are retained permanently during growth 


throughout the life of a fish (Campana and Neilson, 
1985). Such chemical histories have been shown to cor- 
respond to both environmental and physiological condi- 
tions experienced by individuals during life (Campana, 
1999). Heimbrand et al. (2020) demonstrated that the 
uptake of certain trace elements, magnesium (Mg) and 
phosphorus (P), in otoliths of Atlantic cod had seasonal 
oscillations, resulting in distinct elemental peaks that 
were correlated with age. Htissy et al. (2021) found simi- 
lar results for the uptake of Mg, P, and manganese (Mn) 
and proposed counting of enrichment peaks for select 
elements as a means of age determination. Also, age esti- 
mates from the use of 2-dimensional (2D) geochemical 
maps of the entire surface of otoliths have been shown 
to correlate with those from the use of visual annuli 
(Ulrich et al., 2009), providing confidence in the use of 
elemental peak counting techniques for age determi- 
nation (Heimbrand et al., 2020). Methods for counting 
peaks in uptake of trace elements in otoliths show prom- 
ise as alternative aging techniques for species for which 
the traditional method of counting growth increments 
is prone to error because of false or weak annuli, as has 
been the case for the white hake. 

The purpose of this study was to develop models for esti- 
mation of age and growth of male and female white hake 
by using traditional growth increment counting informed 
by core-to-edge geochemistry and by using 2D geochemical 
mapping. We hypothesized that growth trends for white 
hake would be sex specific, with females growing larger than 
males as a reproductive benefit, on the basis of the results of 
previous studies in waters of Canada (Hunt, 1982; Markle 
et al., 1982). From the use of our generated growth models, 
we also report estimates for both length and age at maturity 
for both sexes. Overall, our results can be used to inform 
decisions concerning techniques for aging white hake in the 
future, with regard to factors including time, operation scale, 
and the costs associated with laser ablation inductively cou- 
pled plasma mass spectrometry (LA-ICP-MS). 


Materials and methods 
Sample collection 


White hake are sampled by the Maine Department of 
Marine Resources (DMR) during the biannual Maine— 
New Hampshire Inshore Trawl Survey. Conducted since 
the spring of 2000, this fishery-independent, multispe- 
cies trawl survey involves sampling inshore waters out 
12 nautical miles from the coast from Seabrook, New 
Hampshire, to the Canada—United States border. Data 
from the surveys in inshore waters of New Hampshire 
and Maine are used to generate species abundance indi- 
ces and enhance management practices across fisheries, 
and specimens caught during the surveys are used in bio- 
logical studies. Survey gear consists of a modified shrimp 
net, selected to minimize habitat disruption and sam- 
ple multiple species of varying sizes. Fishermen in this 
region collaborated with the DMR on the trawl survey 
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design and selection of gear type. The aim for tows of the 
trawl net is a period of 20 min (mean: 18.8 min [stan- 
dard deviation (SD) 2.2]), with distances of tows ranging 
from 0.29 to 1.51 nautical miles. Although the survey is 
conducted biannually, sagittal otoliths are collected from 
sampled white hake only during the fall, when white hake 
are primarily (~87%) caught. The sampled white hake are 
measured for total length (TL) to the nearest 1 cm and for 
weight to the nearest whole gram, their sex is determined, 
and their maturity status is assessed (see the “Maturity” 
section) prior to extraction of otoliths. 


Otolith preparation 


Otoliths were mounted in Buehler EpoThin 2’ epoxy resin 
(Buehler Ltd., Lake Bluff, IL), then transversely cut, with 
a target thickness of 0.5-0.75 mm, by using an IsoMet 
Low Speed Saw (Buehler Ltd.) with a diamond blade. 
After cross-sectioning, otolith samples were polished with 
a range of low-grit sandpaper (300—60,000 grit), to smooth 
the surface and increase optical clarity for counting of 
growth increments and to prepare otoliths for trace ele- 
mental analysis. 

All otoliths were assigned an age by the primary 
reader (see the “Aging” section) prior to elemental analy- 
sis. Otolith samples were mounted to glass petrographic 
slides with Crystalbond thermoplastic cement (Aremco 
Products Inc., Valley Cottage, NY), then sent off for trace 
elemental analysis. In this study, we used 3 different 
approaches for trace elemental analysis: line transects 
(number of otolith samples [n]=43), partial 2D maps 
(n=30), and complete 2D maps (n=7). This analysis was 
conducted at 3 separate geochemistry facilities: Univer- 
sity of Texas at Austin, University of Maine, and Oak 
Ridge National Laboratory. Line transects consisted of 
a single line scan run from the otolith core to the edge, 
along the primary aging axis. Partial 2D maps delin- 
eated a region consisting of the otolith’s core and dor- 
sal surface, and complete 2D maps delineated the entire 
otolith’s surface. The high-resolution 2D elemental maps 
provide spatial detail to better understand the regional 
compositional context of the elemental time series 
obtained with the single core-to-edge line scans. These 3 
approaches for trace elemental analysis varied in effort 
(time) and were selected as a means to accurately iden- 
tify chemical increments. We chose multiple geochemis- 
try facilities because of logistical concerns about funding 
and equipment availability. 


Aging 


The primary reader for this study (senior author) assigned 
ages for all otoliths (n=550). Samples were selected across 
a range of years, from 2007 to 2021 (Suppl. Table 1) (online 
only), in which the DMR fall inshore trawl survey was 
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conducted, with the assumption that variation in growth 
patterns between years sampled was negligible. Aging of 
white hake was conducted by using a Leica S9 D micro- 
scope and LAS X imaging software (Leica Microsystems, 
Wetzlar, Germany), following techniques outlined by Hunt 
(1982). Samples were submerged, with the sulcus groove 
facing down and the dorsal side to the right (Fig. 1), in a 
small, blacked-out, water-filled petri dish that sat on the 
microscope stage during the aging process. Otoliths were 
illuminated by using reflected light, then imaged with a 
gray scale to increase growth increment contrast. Addi- 
tional band-clarifying methods, such as baking, burning, 
and staining, have not been proven to enhance contrast 
between growth increments (Hunt, 1982) and were not 
used in this study. 

The most consistent counts of growth increments were 
done along the primary aging axis, which starts at the 
otolith core and runs along the dorsal axis (Fig. 1) (Clay 
and Clay, 1991). Attempts to assign ages by counting in 
other core-to-edge directions, such as the secondary aging 
axis that runs parallel to the sulcus groove (Fig. 1), were 
necessary when poor growth incrementing on the dorsal 
plane occurred. Following previous research on white hake 
(Beacham and Nepszy, 1980; Hunt, 1982; Clay and Clay, 
1991) and other groundfish species (Gauldie and Nelson, 
1990); Stevenson and Campana, 1992; Ferri, 2023), we 
assumed that a pair of translucent and opaque growth 
increments represents 1 year of growth. 

The secondary reader (J. Price) assigned ages for a ran- 
dom subset (n=130) of otoliths. The precision of agreements 
between reader-assigned ages were evaluated with 2 indi- 
ces, percent agreement and Chang’s coefficient of varia- 
tion (CV) (Chang, 1982), and bias of disagreements were 
evaluated with 3 tests of symmetry, Evans and Hoenig’s 
(1998), Bowker’s (1948), and McNemar’s (1947) tests, fol- 
lowing the recommendations of McBride (2015). Compar- 
isons between ages assigned by the primary reader and 
ages estimated by analyzing otolith geochemistry (see the 
“Klemental data analysis” section) used the same 2 preci- 
sion indices and 3 bias tests. All such analyses and data 
plotting were done by using the FSA package (vers. 0.9.1; 
Ogle et al., 2021) in RStudio (vers. 2022.07.2; RStudio 
Team, 2022) (Ogle, 2016). 

The final age assignments, which involved adjust- 
ments for age differences between sampling months, 
were determined after accounting for between-reader 
discrepancies. Although spawning seasonality of white 
hake in the GOM has not been observed, juveniles of 
this species are thought to hatch in the late spring or 
early summer months in the GOM (Ames, 2012); there- 
fore, birthdates were assigned as 1 June and their date 
of capture was factored into the age estimate by using 
the following equation: 


Age = Age assignment 
+ ((Month of capture — 6)/12). (1) 


By using these final age estimates, von Bertalanffy 
growth curves for male and female white hake were fit. 
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Figure 1 


Image of a transversely sectioned sagittal otolith from a female white hake (Urophycis 
tenuis) sampled on 26 September 2012 in inshore waters of southern Maine, with proper 
orientation for aging. This fish had a total length of 42 cm, a weight of 0.5 kg, and an 
estimated age of 3.3 years. The primary aging axis is the axis where growth increments 
were primarily sampled, and the secondary aging axis is the axis where growth incre- 
ments were sampled if the primary axis was difficult to interpret or to check the count 


of increments on the primary axis. 


These growth equations were modeled using the following 
formula (von Bertalanffy, 1934): 


15 (pes (1 & ery (2) 


where L(¢) = the length (in centimeters) of an individual at 
a certain time ¢ (in years); 
L, = the hypothetical average maximum length; 
k = the rate of growth to L,,; and 
ty = the age at which length is zero. 


The Gompertz growth curve was also explored for these 
age estimates but produced a weaker model fit (AIC: 1460 
for males and 1726 for females) compared with that of the 
von Bertalanffy growth curve (AIC: 1457 for males and 
1579 for females), and it was not investigated further. A 
full combination of fixed variables from growth equations 
for both males and females was explored with the von 
Bertalanffy growth function, as well (Table 1). 


Trace element analysis 


For trace elemental analysis, a random subset of 73 oto- 
lith samples was chosen from the 550 samples for which 
traditional aging was done in this study. This subset was 
from fish that ranged in sex (40 females and 33 males), 
maturity stage, age, sampling year, and age readability 
of their otoliths. Across laboratories, analytes included 
isotopes of Mg, calcium (Ca), Mn, strontium (Sr), and 
barium (Ba): **Mg, 43Cq, °°Mn, °8Sr, and /°°Ba. Our anal- 
yses entailed otolith core-to-edge line scans along the 


Table 1 


Comparison of the various configurations of the von 
Bertalanffy growth function used in this study. Models 
differed in which, if any, parameters were fixed. The 
parameters include the hypothetical average maximum 
length (L,,), the rate of growth (k), and the age at which 
length is zero (t,). The asterisk (*) indicates the model 
with the lowest Akaike information criterion (AIC) and, 
therefore, the selected model. Models were fit to data 
from analysis of otoliths from white hake (Urophycis 
tenuis) captured in New Hampshire and Maine during 
fall from 2007 through 2021. 
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primary aging axis (Fig. 1), as well as partial (from the 
core to the dorsal surface only) and complete 2D maps. 
For complete 2D mapping, contiguous parallel line scans 
were run over the entire sample surface, allowing for 
observation of elemental patterns, such as oscillations in 
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concentrations of elements consistent with growth incre- 
ments (Heimbrand et al., 2020). The 2D maps clearly 


show high enrichment at the otolith core across elements - 


(Fig. 2); therefore, line scans were trimmed from the first 
growth increment to the edge, in an effort to reduce the 


Figure 2 


2D maps of the geochemistry of the otolith of a female 
white hake (Urophycis tenuis) sampled on 24 October 2012 
in inshore waters of northern Maine. Peaks in concentra- 
tions of isotopes of magnesium (Mg), calcium (Ca), manga- 
nese (Mn), strontium (Sr), and barium (Ba)—**Mg, °°Mn, 
88Sr and }°°Ba—in parts per million (ppm) were analyzed. 
This fish had a total length of 85 cm, a weight of 5.1 kg, 
and an estimated age of 7.3 years. 


Fishery Bulletin 122(1-2) 


potential for misreading the highly enriched core as a 
growth increment. 


University of Texas at Austin Elemental variations along 
a line transect on each otolith (n=43) of white hake were 
measured through LA-ICP-MS, by using an ESI NWR193 
excimer LA system (193-nm wavelength and 4-ns pulse 
width; Elemental Scientific Inc., Omaha, NE) coupled 
with an Agilent 7500ce ICP-MS (Agilent Technologies, 
Santa Clara, CA). The LA-ICP-MS system is equipped 
with a large format, 2-volume sample cell with fast wash- 
out (<1 s) that accommodated all samples and standards 
(ECRM-752-NP, USGS MACS-3, and NIST 612) in a single 
cell loading. The system was optimized daily for sensitivity 
across the atomic mass unit range and for low oxide pro- 
duction (mean ratio of thorium oxide to thorium [ThO/Th]: 
0.40 [SD 0.01]) by tuning on a reference material (NIST 
612), and these parameters were checked with trial tran- 
sects on representative samples. 

Following pre-ablation (75-um spot, 50-ym/s scan rate, 
2-J/em energy density [fluence]) to remove shallow sur- 
face contaminants, scans along transect lines on sec- 
tioned otoliths were run from the core to the marginal 
edge, by using a 25-by-50-um rectangular aperture, 
with a long axis perpendicular to the radial growth axis, 
a 5-um/s scan rate, a fluence of 2.73 J/em (SD 0.05), a 
20-Hz repetition rate, and a carrier gas flow of 0.80 L/min 
for both argon and helium. The quadrupole duty cycle of 
0.4152 s corresponded to 94% of measurement time, with 
a corresponding linear sampling rate of 2.08 m/s, equiv- 
alent to 12 measurements within the 25-ym-high aper- 
ture footprint. 

Measured intensities were converted to elemental con- 
centrations (in parts per million) by using iolite 4 software 
(Paton et al., 2011); **Ca was used as the internal stan- 
dard element with the assumption of a Ca index weight 
percentage of 38.3% (Serre et al., 2018; Hiissy et al., 2021). 
Carbonate pressed pellet USGS MACS-3 was used as the 
primary calibration material, and NIST 612 and ECRM- 
752-NP were measured as external reference materials. 
Analyte recoveries were typically within 7% and 9% of the 
reference values for NIST 612 (n=33) and ECRM-752-NP 
(n=32), respectively. The grand average of secondary 
standard (n=27) recovery fractions for all elements was 
typically within 7% of preferred values in the GeoReM 
database (Jochum et al., 2005) for NIST 612 and within 
14% of values for USGS MACS-3. Typical (mean) analyte 
concentrations over otolith transects were tens to hun- 
dreds of times higher than the limits of detection for all 
elements. 


University of Maine Elemental variations along multiple 
line transects covering the core and dorsal surface of oto- 
liths of white hake for partial 2D maps (n=27) and cov- 
ering the entire surface of otoliths for complete 2D maps 
(n=4) were measured through LA-ICP-MS at the Micro- 
Analytical Geochemistry and Isotope Characterization 
Laboratory, by using an ESI NWR193UC excimer LA sys- 
tem (Elemental Scientific Inc.) coupled to an Agilent 8900 
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triple quadrupole ICP-MS (Agilent Technologies) and fol- 
lowing the methods described in Cruz-Uribe et al. (2021) 
and Walters et al. (2022). The fast washout mapping setup 
in the MicroAnalytical Geochemistry and Isotope Char- 
acterization Laboratory consisted of 1-mm ID tubing con- 
nected directly from the analytical cup of a TwoVol3 LA 
chamber to a dual concentric injector ICP-MS torch (DCI1) 
from Elemental Scientific Inc.. A 10-by-10-yym square spot 
was rastered in parallel lines across the sample surface at 
a repetition rate of 200 Hz and a scan speed of 500 pm/s 
with a fluence of 2 J/cm. 

Two lines of NIST SRM 610 glass were ablated at the 
beginning and end of map acquisition. The USGS GSD-1G 
basalt glass and a speleothem nano-powder pressed pellet 
(KCSp-1NP) from the International Association of Geoana- 
lysts (IAG) were used as quality control materials. Element 
recoveries were within 10% of preferred values for GSD-1G 
in the GeoReM database. Measured trace element values 
for KCSp-1NP were within the uncertainty (SD of 2 on rep- 
licate analyses) of the preferred values in the IAG report 
for this material. Time-resolved signals were processed by 
using the Trace Elements data reduction scheme in iolite 4, 
and trace element maps were produced by using the Cell- 
Space module in iolite 4 (Woodhead et al., 2007; Paton et al., 
2011; Paul et al., 2012; Petrus et al., 2017). Trace element 
mass fractions were determined relative to the NIST SRM 
610 reference glass with the assumption of an internal 
standard of a mass fraction of Ca of 40.04%. 


Oak Ridge National Laboratory Elemental variations on 
the entire surface of otoliths of white hake for complete 
2D maps (n=3) were measured through LA-ICP-MS, by 
using an ESI imageGEO193 laser (Elemental Scientific 
Inc.) equipped with a TV3 ablation cell coupled to a 
Thermo Scientific iC AP TQ ICP-MS (Thermo Fisher Sci- 
entific, Waltham, MA). A 10-by-10-um square spot was 
rastered in parallel lines across the sample surface at a 
repetition rate of 200 Hz and a scan speed of 125 pm/s 
with a fluence of 5 J/em. The USGS GSD-2G basalt glass 
and a speleothem nano-powder pressed pellet (KCSp- 
INP) from the IAG were used as quality control mate- 
rials. Measured trace element values for GSD-2G and 
KCSp-1NP were within the uncertainty (SD of 2 on 
replicate analyses) of the preferred values in the IAG 
reports for these materials. 


Elemental data analysis 


Generated 2D maps (partial and complete) were ana- 
lyzed by using iolite 4 to investigate spatial variations 
in elemental concentrations. In these 2D maps, only Mn 
had patterns of growth incrementing across individuals 
(Fig. 3); therefore, only Mn was investigated in the line 
scan analysis. 

Line scan elemental data were analyzed by using the 
geplot2 package (vers. 3.3.6; Wickham, 2016) in RStudio. 
All derived elemental time series were smoothed by con- 
secutive moving median and average filters with a 25- 
point boxcar width (equivalent to a distance of 175 pm), 


resulting in smooth, locally weighted signals free from 
high-frequency outliers. Signals were converted to dis- 
tances (in micrometers) from the core on the basis of the 
scan rate and duty cycle. Concentrations of Mn were nor- 
malized across individuals, to apply the same parameters 
for counting peaks in elemental concentrations across indi- 
viduals. Peaks in concentrations of elements were counted 
by using the stat_peak and stat_valley functions in RStu- 
dio, with a peak width (span) of 101, peak height (thresh- 
old) of 0.15, valley span of 151, and valley threshold of 0.80 
(see the R script in Supplementary Material [online only]), 
as has been done similarly in previous research (Htissy 
et al., 2021). Peaks in concentrations of elements between 
2 valleys were counted as 1 growth increment, with mul- 
tiple peaks between 2 valleys considered a single growth 
increment (Fig. 3). As was done previously for agreements 
in age assignments between the primary and secondary 
readers, tests of both precision and symmetry were con- 
ducted between age estimates based on assignments of the 
primary reader and counts of peaks in Mn concentrations, 
generating an age bias plot for both males and females. 
The methods used were the same as those used in the pre- 
vious precision and symmetry analysis. 


Maturity 


Classification of the maturity of white hake followed 
the protocol of Burnett et al. (1989), with a final desig- 
nation of either immature or mature male or female. The 
maturity of a few fish (n=4) could not be classified. The 
maturity for the remaining aged white hake (n=546) was 
analyzed. The median length (L;,) and age (A;,) at matu- 
rity were estimated by fitting binomial logistic models to 
sex-specific data with the sizeMAT package (vers. 1.1.2; 
Torrejon-Magallanes, 2020) in RStudio, and 95% confi- 
dence intervals were estimated by using bootstrapping 
(Torrejon-Magallanes, 2020). Sexual dimorphism, which 
had not been reported for this species in the GOM prior 
to our study, was investigated for both the y-intercept and 
slope by using a chi-square distribution summarized in an 
analysis of variance table. 


Results 


From 2007 through 2021, 550 of the white hake caught 
during the fall as part of the Maine-New Hampshire Inshore 
Trawl Survey were sampled for this study (Suppl. Fig. 1) 
(online only). Sampled white hake ranged in age between 0.3 
and 10.3 years for females (n=291) and between 0.3 and 9.3 
years for males (n=259) (Suppl. Table 2) (online only). The 
average female in this study was 41.6 cm TL, 0.9 kg, and 
3.5 years old, and males were on average 32.3 cm TL, 0.4 kg, 
and 2.7 years old. The primary and secondary readers did 
not discern growth variation among years sampled, lending 
confidence to our initial assumption that growth variation 
between year classes is negligible. 

The precision of paired age comparisons between age 
readers is acceptable but indicates the difficulty of aging 
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Figure 3 


Normalized concentrations of the manganese isotope °’Mn in the otoliths of (A) a male and 
(B) a female white hake (Urophycis tenuis) sampled on 4 October 2021 and 28 September 
2012, respectively, in inshore waters of northern Maine. Concentrations are based on data 
from laser ablation line scans. Green (peaks) and red (valleys) dots are annotated on the 
transect of the line scan. The red line on the otolith identifies the line scan path, and the 
white and yellow dots are counted winter and false winter growth increments identified 
by the primary reader. The male had a total length of 33 cm, a weight of 0.3 kg, and an 
estimated age of 3.3 years, and the female had a total length of 102 cm, a weight of 9.3 kg, 
and an estimated age of 10.3 years. 


both sexes of white hake (e.g., Chang’s CVs between 5 and Of all trace elements for which concentrations were 
10) (Table 2). Statistically significant bias was evident quantified, only Mn had concentration patterns that corre- 
between ages assigned by readers (e.g., P<0.01 from sponded to presumed summer and winter growth incre- 
Evans and Hoenig’s test of symmetry), but the average menting (Figs. 2 and 3). The precision of paired age 
differences by age class were small and without trend comparisons between counts of peaks in concentrations of 


(Fig. 4, A and B). Mn and counts of growth increments by the primary 
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Table 2 


Results from evaluation of precision and symmetry of age estimates for white hake (Urophycis 
tenuis) captured in inshore waters of New Hampshire and Maine from 2007 through 2021 in the 
fall. Age assignments made by the primary reader of growth increments on otoliths were com- 
pared either to age assignments made by the secondary reader (Secondary) or to age estimates 
based on peaks in concentrations of manganese in otoliths (Mn peaks). Precision was examined 
with 2 indices: percent agreement (PA), which is reported for both the exact year of the estimate 
and +1 year of the estimate, and Chang’s coefficient of variation (CV). Three symmetry tests 
were used. n=number of otoliths. 


PA (%) 
Method PA (%) (+1 year) CV Symmetry test 


Secondary 71.01 98.55 9.38 McNemar’s 
Evans and Hoenig’s 
Bowker’s 
Mn peaks McNemar’s 
Evans and Hoenig’s 
Bowker’s 
Female Secondary ; McNemar’s 
Evans and Hoenig’s 
Bowker’s 
Mn peaks McNemar’s 
Evans and Hoenig’s 
Bowker’s 


Primary (years) Primary (years) 


Figure 4 


Age bias plots of the differences in ages between assignments by the primary and secondary readers of otoliths from 
(A) male (number of otoliths [n]=69) and (B) female (n=61) white hake (Urophycis tenuis) and between assignments 
by the primary reader and estimates based on counts of peaks in the concentration of the manganese istotope °°Mn for 
otoliths of (C) male (n=33) and (D) female (n=40) white hake. Fish were captured in inshore waters of New Hampshire 
and Maine from 2007 through 2021 in the fall. 
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reader is acceptable but also indicates the difficulty of 
aging both sexes of white hake (e.g., Chang’s CVs between 
9 and 10) (Table 2): No statistically significant bias is evi- 
dent (e.g., P>0.05 from Evans and Hoenig’s test of symme- 
try), with minor disagreement in age assignments across 
lengths for both sexes (Fig. 4, C and D). 

The estimated growth curves differed by sex. The 
growth curve was L(t)=110(1-e°13*°*) for males and 
L(t)=140(1-e7°: 113-93) for females. The von Bertalanffy 
growth function is best described with a fixed value for k 
(0.113), a smaller ¢,) value for males (—0.45 year [SD 0.16]) 
than for females (0.30 year [SD 0.13]), and a larger L,, 
value for females (140 cm TL [SD 29.0]) than for males 
(110 cm TL [SD 19.0]) (Table 1, Fig. 5). 

Median length and age at maturity also differed by sex. 
Female L;, was 10 cm longer than male Lz, (47.4 cm TL 
[SD 2.3] versus 37.4 cm TL [SD 2.4]). Female A;) was 0.9 
year older than male Az, (4.2 years [SD 0.2] versus 3.3 
years [SD 0.3]) (Fig. 6). Model slopes were not significantly 
different between sexes (P>0.05), although y-intercepts 
were (P<0.01). 


Discussion 
This study is the first to demonstrate that sex-specific growth 


models are most suitable for white hake in U.S. waters of 
the GOM and the first to formally test the hypothesis of 


Females (n=291) 
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sexually dimorphic growth with an AIC approach. In pre- 
vious work in the Gulf of St. Lawrence, Hunt (1982) found 
that females attained a larger size at age than males, with 
estimates for the parameter L,, very similar to estimates 
from our study (Suppl. Table 2) (online only). 

The variation in L,, between females (140 cm TL) and 
males (110 cm TL) aligns with previous reports of the lifes- 
pan of females exceeding 20 years and that of males never 
exceeding 8 years (Burnett et al., 1984; Langton et al., 
1994). Although in our study females reached a maximum 
age of only 10.3 years and maximum sizes of only 102 cm 
TL and 9.3 kg, larger and older individuals in the popula- 
tion likely exist. The age of a male white hake sampled in 
our study exceeds the previous maximum age observed for 
this species (8 years) (Burnett et al., 1984; Langton et al., 
1994), with an estimate of 9.3 years and measured sizes of 
85 cm TL and 5.0 kg. This individual was ~2 years older and 
17 cm TL longer than the next largest male in our study, 
indicating that this specimen was anomalous in the inshore 
environment. The DMR, as part of the inshore trawl sur- 
vey, does not sample beyond 12 nautical miles from shore, 
into areas possibly inhabited by white hake that are larger 
and older than those found in shallow waters (Chang et al., 
1999; Ames, 2012). Overall, our results support the notion 
of sexually dimorphic growth that has been previously sug- 
gested for this species (Burnett et al., 1984; Langton et al., 
1994; Haedrich, 2003) but not previously modeled for the 
GOM (Suppl. Table 2) (online only). 


Males (n=259) 
EO (emyp= 110 = e50 eesevert=):+ 0.45)) 


Age (years) 


Figure 5 


von Bertalanffy growth curves for both male (number of fish [n]=259) and female (n=291) white 
hake (Urophycis tenuis) caught in New Hampshire and Maine during 2007-2021 in the fall, based 
on age estimates from counting of visual growth increments in otoliths. Sex-specific von Berta- 
lanffy growth functions are also provided. L(t)=length at time ¢. 
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Figure 6 
Maturity at length for (A) female (number of fish [n]=289) and (B) male (n=257) white hake (Urophycis tenuis) and 
maturity at age for (C) female and (D) male white hake captured in New Hampshire and Maine from 2007 through 
2021 in the fall. The curved dashed lines indicate 95% confidence intervals (CIs). The vertical and horizontal dashed 
lines indicate the median length or (Z;)) or age (As,) at maturity. 


We investigated the possible seasonal and annual pat- 
terns of concentrations of Mg, Mn, Sr, and Ba in otoliths of 
white hake by using 2D maps. These elements were chosen 
because of their prominence in previous otolith age and 
growth studies. We hypothesized that distinguishable vari- 
ation in the concentration of one or more of these elements 
over time may be used as an indicator of age. Our hypothe- 
Sis is supported by the finding that oscillations in the con- 
centration of Mn, among those of other elements, most 
clearly had patterns in 2D maps similar to those of growth 
increments in otoliths, and we inferred that the patterns 
correspond to changes in Mn concentration in otoliths 
between seasonal growth increments (summer and winter). 
Variation in concentrations of Ba, Mg, and Sr each had 
varying patterns similar to those of growth increments, but 
the patterns were not consistent across individuals and, 
therefore, could not be attributed to growth (Suppl. Fig. 2, 
Suppl. Table 3) (online only). We speculate that variations in 
the concentration of these elements across individuals may 
be due to a combination of environmental (e.g., temperature 
and water chemistry) and physiological (e.g., diet and repro- 
duction) factors. Research into these variations is ongoing, 
and we selected only Mn for line scan and peak counting as 
a means of estimating ages that could be later compared to 
age assignments by the primary reader of growth incre- 
ments on otoliths. We made these comparisons by using 


symmetry and precision analyses, finding that age estimates 
from analysis of Mn concentration were ~100% in agree- 
ment, +1 year, with ages assigned by readers. Again, we 
highlight the lack of age validation for this species and note 
that validation would allow increased confidence in the 
notion that Mn oscillations align with age. 

Uptake of Mn into otoliths occurs as a substitute for 
Ca during biomineralization (Thomas et al., 2017; Hiissy 
et al., 2021), with reports from previous research also 
describing correlations between Mn and matrix proteins 
(Thomas et al., 2017). The element Mn has been shown 
to be a cofactor of extracellular serine and threonine pro- 
tein kinase (Fam20C) (Tagliabracci et al., 2012), a biomin- 
eralization matrix protein. Growth increment patterns 
observed in 2D maps led us to hypothesize that Mn is a 
tracer of historical protein content, with seasonal oscilla- 
tions in protein incorporation driven by growth, account- 
ing for consistent age estimations across size classes. 

In previous studies, elevated levels of Mn in otoliths 
cores have decreased through time, possibly in relation to 
maternal offloading and metabolism (Brophy et al., 2004; 
Ruttenberg et al., 2005; Miller, 2009; Friedrich and Halden, 
2010; Clarke et al., 2011; Limburg et al., 2015; Hughes 
et al., 2016). Our 2D maps and data from line scans indi- 
cate comparable spatial variations on otoliths, with oto- 
lith cores enriched with Mn and the concentration of this 
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element decreasing toward the edge of otoliths over suc- 
cessive inferred summer and winter growth increments. 

In previous trace element research, Mn concentration 
has been shown to correlate with both hypoxic redox 
conditions that promote solubility of aqueous Mn** and 
somatic growth (Limburg et al., 2011; Mohan and Walther, 
2015, 2016). Hypoxia is uncommon in the GOM, as its 
waters are cold and well mixed and its substrates are 
rocky (Hale and Heltshe, 2008; Irish et al., 2008); there- 
fore, we attribute patterns of Mn concentration in oto- 
liths of white hake to somatic growth. Heimbrand et al. 
(2020) made similar conclusions, ranking signals of Mn 
in otoliths with a high readability score for Atlantic cod 
and attributing distinct oscillations in concentration to 
growth. Variations in readability of Mn signals has been 
theorized to be a result of major influx in oxygen in the 
Baltic Sea, making Mn mineralization less reliable as a 
recorder of age compared with uptake of Mg and P, as 
substitution for Ca becomes altered (Heimbrand et al., 
2020; Hiissy et al., 2021). Unlike in the Baltic Sea, hypoxic 
redox conditions do not regularly occur in the GOM, 
which has stable oxygen availability (Hale and Heltshe, 
2008; Irish et al., 2008); therefore, somatic growth, which 
influences the seasonal incorporation of Mn in otoliths, is 
theorized to be less disrupted in white hake. In our inves- 
tigation of concentration of Mg in otoliths of white hake, 
we found no discernable patterns indicative of growth, a 
result that varies from those of Heimbrand (2020) and 
Hiissy (2021). We interpret this result as a potential spe- 
cies (white hake versus Atlantic cod) or ecosystem (GOM 
versus Baltic Sea) effect. 

Our analyses did not include P because of a lack of 
access to a standard to calibrate LA, and we strongly advo- 
cate for continued investigations into the concentration 
of P, and other trace elements, in otoliths and into envi- 
ronmental drivers of incorporation not examined in this 
study, for white hake in the GOM and for other species 
and ecosystems. 

Similar to results from our examination of age esti- 
mates and growth, maturity schedules differed by sex. 
Previously, white hake had been reported to reach Lz, 
at 32.7 cm TL for males and at 35.1 cm TL for females, 
with an estimated A;, of 1.4 years for both sexes (O’Brien 
et al., 1993). In comparison to those previously described 
differences, our observations indicate a high degree of sex- 
ual dimorphism in both length (L;, differed by 10 cm TL) 
and age (A,, differed by 1 year). Larger sizes of females, 
in comparison to the sizes of males, may indicate that 
female white hake increase fecundity with size, as noted 
by Beacham and Nepszy (1980). Delayed length and age at 
maturity has been reported previously as a sign of stock 
recovery, declining fishing pressure, or increased repro- 
ductive success for northern pike (Esox lucius) or chum 
salmon (Oncorhynchus keta) (Carlson et al., 2007; Edeline 
et al., 2007; Fukuwaka and Morita, 2008). The current 
understanding of the reproduction of white hake is lim- 
ited; therefore, continued research is necessary to be able 
to pinpoint population trajectories. Although white hake 
have been previously described as highly fecund (Han and 
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Kulka, 2008), significant spawning in the GOM has yet to 
be observed (Fahay and Able, 1989; O’Brien et al., 1993). 


‘Continued investigations into the maturity patterns of 


white hake by using historical data sets may help describe 
the future trajectory of stocks. 

In this study, we compared techniques of aging otoliths 
from white hake, a traditional method of counting visual 
growth increments and a more recent method of analyzing 
geochemical markers. We found no variation in cost, effort, 
and time related to otolith preparation for each method. 
The traditional method of aging was more accessible (had 
lower costs and was faster) at a high volume of samples 
than the other investigated method. Aging through the 
use of otolith geochemistry was done quickly once ele- 
mental data were received from laboratories, but data 
collection and processing requires specialized instruments 
(had higher costs and was slower). The cost of LA-ICP-MS 
varies greatly between laboratories and equipment used, 
but results from using it have been promising in aging 
difficult species, like white hake. On the basis of our find- 
ings, we suggest that geochemical aging be considered 
for species for which aging through traditional reading of 
growth increments has proven difficult. Our approach of 
selecting white hake in a range of varying sizes provided 
us the ability to discern which elements had oscillations 
in their concentrations in otoliths that correlated to age, 
and the age estimates from analysis of these oscillations 
can be used moving forward as a reference for aging white 
hake by using traditional methods. With advancements in 
LA-ICP-MS technology, lowering of operation costs, and 
continued research, we see the potential for this approach 
to become an increasingly popular alternative method of 


aging. 


Conclusions 


Results from this study of white hake indicate the useful- 
ness of analyzing otolith geochemistry as an alternative 
method of assigning age. Herein, we report accepted mea- 
sures of precision and bias from a traditional compari- 
son of between-reader agreement in age estimates and 
from comparisons of annulus counts of one reader to age 
estimates inferred from peaks in concentrations of Mn 
discerned through 2D mapping of otolith geochemistry. 
This expanded approach is likely to become more com- 
mon (e.g., Heimbrand et al., 2020; Hiissy et al., 2021), 
although we cannot predict which elements will produce 
reliable annual patterns for different species, especially 
among different ecosystems. We emphasize that anal- 
ysis of otolith geochemistry is worth consideration as 
a means to gain confidence in age assignments for dif- 
ficult to age species; of course, growth validation, cost, 
and effort requirements of this method must be part of 
such consideration. The ability to increase certainty in 
age estimates for a species allows fishery managers to 
confidently determine growth rates, longevity, and repro- 
ductive output, all critical information for accurate stock 
assessment. 
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Resumen 


La merluza blanca (Urophycis tenuis) es un pez de fondo 
distribuido por todo el Golfo de Maine. Las capturas se 
recomiendan con base en la evaluacion de stock utilizando 
modelos de dinamica de poblaciones estructurados por edad; 
sin embargo, datar los otolitos es un reto debido a la falta de 
claridad de los incrementos de crecimiento. Para abordar esta 
preocupacion, comparamos la consistencia entre las edades 
de recuentos visuales de anillos con las edades asignadas 
con ciclos de las concentraciones de elementos medidos utili- 
zando ablacién laser mediante espectrometria de masas con 
plasma acoplado inductivamente. Probamos la hipétesis de 
que las oscilaciones tanto en las condiciones ambientales 
como en la fisiologia interna a lo largo del tiempo influyen en 
la absorcién de elementos durante la mineralizaci6én del oto- 
lito. Las concentraciones de manganeso, en comparacion con 
las de otros elementos traza investigados (magnesio, estron- 
cio y bario), tuvieron la correlacién mas prometedora con los 
incrementos visuales de crecimiento (~100% de acuerdo con 
la edad, +1 afio), ofreciendo una herramienta adicional para 
mejorar la identificacién de los incrementos. Examinamos 
550 otolitos colectados durante 2007-2021, encontrando que 
la merluza blanca vivia un maximo de 10.3 anos con dimorf- 
ismo sexual en longitud y edad maxima. Utilizando funciones 
de crecimiento de von Bertalanffy, L(t)=110(1-e 013°”) 
para machos y L(t)=140(1-e 01°) para las hembras 
(donde Lit) es la longitud en el tiempo t), se calcularon la 
talla y la edad de madurez de los machos (37.4 cm de longi- 
tud total [LT], 3.3 anos) y las hembras (47.4 cm LT, 4.2 afios). 
Estos resultados demuestran que la geoquimica de los otoli- 
tos puede utilizarse para mejorar la exactitud y precision de 
la estimacion de la edad y madurez de los peces, incluso en 
especies dificiles. 
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Guidelines for authors 


Contributions published in Fishery Bulletin describe 
original research in marine fishery science, fishery engi- 
neering and economics, and marine environmental and 
ecological sciences (including modeling). Preference will 
be given to manuscripts that examine processes and 
underlying patterns. Descriptive reports, surveys, and 
observational papers may occasionally be published but 
should appeal to an audience outside the locale in which 
the study was conducted. 

Although all contributions are subject to peer review, 
responsibility for the contents of papers rests upon the 
authors and not on the editor or publisher. Submission 
of an article implies that the article is original and is not 
being considered for publication elsewhere. 

All submissions are subject to a double-blind review 
to remove the identity of author and reviewer during the 
review process. 

Plagiarism and double publication are considered 
serious breaches of publication ethics. To verify the orig- 
inality of the research in papers and to identify possible 
previous publication, manuscripts may be screened with 
plagiarism-detection software. 

Manuscripts must be written in English; authors 
whose native language is not English are strongly advised 
to have their manuscripts checked by English-speaking 
colleagues before submission. 

Artificial intelligence (AI), machine learning, or similar 
algorithmic tools should not be used to write or generate 
scientific content, including text, tables, and figures, for 
manuscripts submitted to Fishery Bulletin. Although some 
authors may use such tools to improve the readability or 
syntax of their manuscript, it must be done with extreme 
caution, and authors must review and edit the result. If 
such a tool is used, authors must declare the name of the 
tool, where it was used, and the purpose for which it was 
used in an “Artificial intelligence declaration” section after 
the “Acknowledgments” section. Authors remain account- 
able for all elements of their manuscript and should state 
that they take full responsibility for the content of their 
manuscript in the AI declaration. 

Once a paper has been accepted for publication, online 
publication takes approximately 2 weeks. 

There is no cost for publication in Fishery Bulletin. 


Types of manuscripts accepted by the journal 


Articles generally range from 20 to 380 double-spaced 
typed pages (12-point font) and describe an original 
contribution to fisheries science, engineering, or econom- 
ics. Tables and figures are not included in this page count, 
but the number of figures should not exceed 1 figure for 


every 4 pages of text. Articles contain the following divi- 
sions: abstract, introduction, methods, results, discussion, 
and conclusions. 


Notes are generally less than 10 double-spaced typed 
pages (12-point font), including the Literature cited sec- 
tion. Like articles, notes describe an original contribu- 
tion to fisheries science. They follow a format similar to 
that for articles: abstract, introduction, methods, results, 
and discussion, but the results and discussion sections 
may be combined and a conclusions section should not 
be included. They should include no more than 2 figures 
or tables (2 of each would be too many). They are distin- 
guished from full articles in that they report a noteworthy 
new observation or discovery—such as the first report of 
a new species, a unique finding, condition, or event that 
expands our knowledge of fisheries science, engineering, 
or economics—and do not require a lengthy discussion. 
Manuscripts on range extensions will not be considered. 


Companion articles should be submitted together and 
are published together as a scientific contribution. Both 
articles should address a closely related topic and may be 
articles that result from a workshop or conference. 


Review articles of exceptional quality may be considered 
occasionally for publication on a case-by-case basis. They 
address a timely topic that is relevant to aspects of fisheries 
science. They should include an abstract, but the format of 
the article, per se, will be up to the author. Please contact 
the scientific editor to discuss your ideas regarding a poten- 
tial review article before embarking on such a project. 


Perspective essays generally are less than 5 double- 
spaced typed pages (12-point font) and provide clarifica- 
tion of misused terms or misunderstood topics in fisheries 
science. Essays are accepted at the discretion of the 
scientific editor. Opinion pieces on policy will not be con- 
sidered. They should include an abstract and a list of the 
literature cited, but no other sections should divide the 
article. Please contact the scientific editor to discuss your 
idea before making a submission. 


Preparation of manuscript 


Title page should include authors’ full names, affiliations, 
mailing addresses, and the senior author’s email address. 


Abstract should be limited to 200 words (one-half typed 
page), state the main scope of the research, and empha- 
size the authors’ conclusions and relevant findings. Do 
not review the methods of the study or list the contents of 
the paper. Because abstracts are circulated by abstract- 
ing agencies, it is important that they represent the 
research clearly and concisely. 


General text must be typed in 12-point Times New Roman 
font throughout. A brief introduction should convey the 
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broad significance of the paper; the remainder of the paper 
should be divided into the following sections: Materials and 
methods, Results, Discussion, Conclusions, and Acknowl- 
edgments. Headings within each section must be short, 
reflect a logical sequence, and follow the rules of subdivision 
(i.e., there can be no subdivision without at least 2 subhead- 
ings). The entire text should be intelligible to interdisciplin- 
ary readers; therefore, all acronyms, abbreviations, and 
technical terms should be written out in full and defined 
the first time they are mentioned. Abbreviations should be 
used sparingly because they are not carried over to index- 
ing databases and slow readability for those readers out- 
side a discipline. They should never be used for the main 
subject (e.g., species or method) of a paper. 

For general style, follow the U.S. Government Publish- 
ing Office Style Manual (2016, available at website) and 
Scientific Style and Format: the CSE Manual for Authors, 
Editors, and Publishers (2014, 8th ed.) published by the 
Council of Science Editors. For scientific nomenclature, 
use the current edition of the American Fisheries Soci- 
ety’s (AFS) Common and Scientific Names of Fishes from 
the United States, Canada, and Mexico and its companion 
volumes (Crustaceans, Mollusks, Cnidaria and Ctenophora, 
and World Fishes Important to North Americans). For spe- 
cies not found in the previously mentioned AFS publica- 
tions and for more recent changes in nomenclature, use the 
Integrated Taxonomic Information System (ITIS, available 
at website), or, secondarily, the California Academy of Sci- 
ences Catalog of Fishes (available at website) for species 
names not included in ITIS. Common (vernacular) names of 
species should be lowercase. Citations must be given for the 
identification of specimens. For example, “Fish species were 
identified according to Collette and Klein-MacPhee (2002); 
sponges were identified according to Stone et al. (2011).” 

Dates should be written as follows: 11 November 2018. 
Measurements should be expressed in metric units, for 
example, “58 metric tons (t);” if other units of measurement 
are used, please make this fact explicit to the reader. Use 
numerals, not words, to express whole and decimal num- 
bers in the general text, tables, and figure captions (except 
at the beginning of a sentence). For example, “We consid- 
ered 3 hypotheses. We collected 7 samples in this location.” 
Use American spelling. Refrain from using the shorthand 
slash (/), an ambiguous symbol, in the general text. 

Cite all software, special equipment, and chemical solu- 
tions used in the study within parentheses in the general 
text, including the version number, company name, and the 
city and state (or nation) of the company headquarters, for 
example, “SAS, vers. 9.4 (SAS Institute Inc., Cary, NC).” 


Word usage and grammar that may be useful are the 
following: 


e Aging 
For our journal, the word aging is used to mean both 
age determination and the aging process (senescence). 
Authors should make clear which meaning is intended 
where ambiguity may arise. 


e Fish and fishes 
The plural of the word fish (a collective noun that 
implies individuals without regard to species) is fish. 
Example: The fish were collected by trawl net. 
Example: The numbers of fish collected that season 
were less than the numbers from previous years. 


The plural for fish species is fishes (a contrived plural 
used by taxonomists to mean several or more fish spe- 
cies) or one can use fish species (which is preferred in 
this journal for clarity across disciplines). 
Example: The fishes of Puget Sound [biodiversity is 
implied] or 
Example: The fish species of Puget Sound [preferred 
plural for clarity across disciplines]. 


e Crab and crabs, squid and squids, etc. 
The plural of the word crab (i.e., many individuals 
without regard to species) is crab. 
Example: The crab were sorted by weight. 
Example: Many red king crab were dying [Many 
individuals of one species of crab.] 


The plural of crab species is crabs (a word used by tax- 
onomists) or crab species (the latter is preferred in this 
journal for clarity). 

Example: These crabs were selected for treatment. 
[Different crab species are implied.]| 

Example: These crab species were selected for 
treatment. [Preferred word choice for clarity.] 

Example: Snow crabs are found throughout the 
North Pacific Ocean and Bering Sea. [There are 
2 species of snow crab; therefore the word crabs can 
be used here.] 

Example: Two species of snow crab are found through- 
out the North Pacific Ocean and Bering Sea. [Pre- 
ferred usage for clarity.] 

Example: Three crabs were selected for treatment. 
[3 species of crab are implied.] 

Example: Three crab species were selected for treat- 
ment. [Preferred word choice for clarity.] 


e We use fisherman and fishermen, not fisher and fishers, 
in this journal. One can always use crew member, vessel 
operator, and angler (the latter for recreational fishing). 


e The definite article with common names of species 
When the singular common name of a species rep- 
resents the entire class or group to which it belongs, 
use the definite article. 

Example: Only one species of the genus Salmo is 
found in the Atlantic Ocean—the Atlantic salmon 
(Salmo salar). 

Example: The sonic emissions of the bottlenose dol- 
phin are complex. 


For plural common names, this rule does not apply. 
Example: Chinook salmon are found throughout the 
Pacific Ocean. 
Example: Bottlenose dolphins are found in temper- 
ate and tropical waters. 
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e Sex 
For the meaning of male and female, use the word sex, 
not gender. Do not write, “fish were sexed.” Write, “sex 
was determined.” 


e Participles 
As adjectives, participles must modify a specific noun 
or pronoun. 

Example: Using mark-recapture methods, these sci- 
entists determined the size of the population. [Correct. 
The participle using modifies the word scientists.] 

Example: These scientists, based on the collected 
data, concluded that the mortality rate of these fish 
had increased. [Incorrect. The scientists were not 
based on the collected data.] 

Example: These scientists concluded, on the basis of 
collected data, that the mortality rate of these fish 
had increased. [Correct. The offending participle has 
been eliminated and an adverbial phrase modifies 
the verb concluded.| 


Equations and mathematical symbols should be set from a 
standard mathematical program (MathType or Equation 
Editor). Equations formatted in LaTex are not accept- 
able. For mathematical symbols in the general text (a, y”, 
mM, +, etc.), use the symbols provided by the MS Word pro- 
gram and italicize all variables, except those variables 
represented by Greek letters and the superscript and 
subscript parts of variables and expressions. Do not use 
photo mode when creating these symbols in the general 
text, and do not cut and paste equations, letters, or sym- 
bols from a different software program. 

Number equations (if there are more than one) for 
future reference by scientists; place the number within 
parentheses at the end of the first line of the equation. 


Literature cited section comprises published works and 
those accepted for publication (in press) in peer-reviewed 
journals. Follow the name and year system for citation for- 
mat in this section (i.e., citations should be listed alpha- 
betically by the authors’ last names, and then by year if 
there is more than one citation by the same author. A list 
of abbreviations for citing journal titles can be found on 
our website. 

Authors are responsible for the accuracy and com- 
pleteness of all citations. Avoid the use of multiple 
citations when a single citation sufficiently supports a 
statement; cite the work that first reported the informa- 
tion that supports a statement, not all of the subsequent 
works. Literature citation format: Authors (last name, 
followed by initials for first name and, if given, mid- 
dle name of first author; then list names of additional 
authors with initials before last names). Year. Title of 
article. Abbreviated title of the journal in which it was 
published. Always include either the range of page num- 
bers (for a journal article) or a total number of pages (for 
a book or other type of publication). List a sequence of 
citations in the general text chronologically, for example, 
“(Smith, 1932; Green, 1998; Smith and Jones, 2015).” 


Acknowledgments should be no more than 6 lines of 
text. Only those who have contributed in an outstanding 
way should be acknowledged by name. For recognition of 
other persons or groups, use a general term, such as crew, 
observers, or research coordinators, and do not include 
names with these terms. 


Digital object identifier (doi) code ensures that a publica- 
tion has a permanent location online. A doi link (which 
may include a doi code) should be included at the end of 
citations of published literature. Authors are responsi- 
ble for submitting accurate doi links. Faulty links will be 
deleted at the page-proof stage. 


Footnotes are used for all documents that have not been 
formally peer reviewed and for observations and personal 
communications, but these types of references should be 
cited sparingly in manuscripts submitted to the journal. 

All reference documents, administrative reports, inter- 
nal reports, progress reports, project reports, contract 
reports, personal observations, personal communications, 
unpublished data, manuscripts in review, and council meet- 
ing notes are footnoted in 10-point font and placed at the 
bottom of the page on which they are first cited. Footnote 
format is the same as that for formal literature citations. A 
link to the online source (e.g., [Available from http://www... , 
accessed July 2017.]), or the mailing address of the agency 
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